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EXECUTIVE  SUMMARY 


This  technical  reDort  presents  the  detailed  methodology  of  various  modeling 
techniques  developed  by  the  Office  of  the  Assistant  for  Study  Support  (OAS)  during 
the  course  of  analytical  studies  on  life  cycle  costs  and  test  and  evaluation. 

These  studies  were  initiated  at  the  request  of  DCS/Development  Plans,  Headquarters 
AFSC  and  were  designed  to  investigate  new  and  innovative  methods  of  reliability 
management  and  to  develop  models  and  methodology  applicable  to  life  cycle  cost  and 
test  and  evaluation  analyses.  The  principal  impetus  in  developing  the  models  came 
from  one  of  the  major  findings  in  the  initial  study  performed  by  OAS  on  test  and 
evaluation  analysis.  This  finding  concerned  the  generally  poor  field  reliability 
experienced  by  the  A-7D  avionics  systems.  Further  investigations  showed  that 
Air  Force  avionics  systems  generally  experience  poor  reliability  in  the  operational 
environment.  There  are  a  multitude  of  reasons  for  this;  one  reason  being  that 
current  avionics  systems  are  extremely  complex  and  low  reliability  is  therefore 
to  be  expected.  This  does  not  mean,  however,  that  low  reliability  of  avionics 
systems  should  be  accepted  as  a  way  of  life.  In  fact,  every  attempt  siio^d  be 
made  to  achieve  the  highest  reliability  possible  commensurate  with  the  complexity 
of  the  system  and  cost  constraints. 

In  probing  deeper  into  reliability  problems  associated  with  Air  Force  avionics 
systems,  one  must  examine  the  regulations,  tests,  and  demonstrations  that  form  the 
basis  for  Air  Force  acceptance  and  procurement  of  avionics  systems.  Air  Force 
Regulation  80-5  states  that  "realistic  and  meaningful  R&M  characteristics  and 
levels  will  be  determined  by  cost  effectiveness  analysis,  reflecting  both  system 
effectiveness  and  life  cycle  costs."  AFR  80-5  states  further  that  "the  establish¬ 
ment  of  minimum  acceptable  R&M  levels  must  be  determined  on  the  basis  of  realistic 
operational  needs."  Unfortunately,  there  has  not  been  an  effective  technique  for 
selecting  levels  of  subsystem  R&M  which  would  optimize  the  total  cost  effectiveness 
of  a  system.  Therefore,  the  desired  system  effectiveness  is  not  being  achieved, 
and  life  cycle  costs  are  greatly  exceeding  desired  levels. 

A  corollary  problem  has  been  the  total  lack  of  similarity  between  contractual 
R&M  requirements  and  the  values  actually  achieved  in  the  field.  Much  of  this 
problem  is  due  to  deficiencies  in  the  test  procedures  used  (normally  from 
MIL-STD-781B  for  electronic  equipment),  but  it  is  also  affected  by  the  lack  of  a 
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rational  approach  to  the  establishment  of  requirements  or  for  evaluating  alter¬ 
natives  when  testing  indicates  that  the  initial  requirements  are  not  being 
achieved.  The  result  has  been  a  spiral  of  higher  and  higher  user  requirements 
for  which  the  developers  have  been  expending  more  and  more  resources  in  a  futile 
attempt  to  achieve.  The  models  and  methodology  developed  herein  represent  a  first 
step  in  attempting  to  bridge  this  gap  between  initial  requirements  and  achievable 
operational  capabilities.  In  addition,  the  models  and  methodology  can  provide 
other  information  of  interest  to  decision-makers  concerned  with  either  development, 
acquisition,  testing,  logistic  support,  or  life  cycle  costs.  The  models  are  quite 
general  and  can  be  applied  at  various  stages  of  system  development.  The  models 
were  developed  in  the  context  of  a  total  system  consisting  of  a  number  of  subsys¬ 
tems.  However,  the  models  can  be  applied  at  the  subsystem  level  by  considering 
the  total  system  to  be  the  subsystem  of  interest  consisting  of  its  components. 

The  initial  model  is  a  Mission  Completion  Success  Probability  (MCSP)  model. 

MCSP  models  are  applied  to  show  the  dependence  of  mission  success  upon  the  aggre¬ 
gate  of  subsystems.  MCSP  models  have  not  been  used  extensively  in  Air  Force 
programs.  In  some  instances  where  MCSP  modeling  techniques  were  employed,  they 
involved  complicated  simulation  methods.  Generally,  an  MCSP  model  involving 
simulation  does  not  readily  lend  itself  to  identification  of  critical  subsystems, 
or  to  evaluation  of  critical  subsystem  improvement.  The  OAS  MCSP  model  developed 
during  this  study  is  a  generalized,  probabilistic  model.  Using  A-7D  data,  the 
utility  of  the  model  has  been  demonstrated  by  ranking  subsystems  according  to 
abort  causing  failures  and  also  in  determining  the  MCSP  enhancement  due  to 
improvements  in  individual  subsystem  reliability.  The  next  step  in  developing 
the  overall  methodology  is  to  consider  reliability  optimization,  i.e.,  the  tradeoff 
between  levels  of  reliability  and  lifetime  support  cost  to  decrease  system  life 
cycle  cost.  MCSP  models  alone  are  inadequate  for  this  task  since  they  do  not 
measure  the  impact  of  subsystem  reliability  levels  on  system  life  cycle  cost. 
ComDining  reliability  optimization  with  MCSP  considerations  leads  to  the  develop¬ 
ment  of  the  Designing  to  System  Performance/Cost  (DSPC)  model. 

The  DSPC  methodology  represents  a  new  and  innovative  approach  to  system  acqui¬ 
sition,  and  preliminary  results  indicate  that  this  technique  will  provide  very 
valuable  information  to  the  decision-maker.  This  methodology  systematically  iden¬ 
tifies  those  subsystem  options  which  provide  the  highest  system  performance  at  any 
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prescribed  level  of  cost  (either  acquisition  cost  or  acquisition  plus  logistic 
support  cost).  The  DSPC  model  is  compatible  with  designing  to  system  cost,  or 
performance,  or  both.  Once  total  system  reliability  specifications  are  established, 
eacn  individual  subsystem  has  a  corresponding  installed  reliability  and  cost  goal 
which  allows  realistic  and  continuous  evaluation  and  adjustment  as  the  subsystem 
is  developed  to  maturity. 

Along  with  the  DSPC  methodology  appropriate  measures  of  effectiveness  must  be 
tailored  to  the  particular  mission  of  interest  and  related  to  system  performance 
parameters.  In  this  way  the  methodology  can  provide  some  of  the  many  inputs  the 
decision-maker  requires.  In  this  report  two  measures  of  effectiveness  for  fighter 
aircraft  are  presented.  In  the  case  of  air-to-ground  fighters,  it  is  shown  that 
an  evaluation  of  the  effectiveness  must  account  for  the  interaction  of  availability, 
abort  probability,  kill  potential,  and  survivability;  and  survivability  is  often 
the  most  dominant  factor.  For  air-to-air  fighters,  the  exchange  ratio  (Red  Aircraft 
destroyed  per  Blue  Aircraft  destroyed)  is  an  important  measure  of  worth,  and  it  can 
be  expressed  as  a  function  of  weapon  effectiveness,  maneuver  capability,  and  first 
shot  probability  with  first  shot  probability  being  the  most  important  parameter. 

As  mentioned  previously,  the  models  and  methodology  can  be  applied  at  various 
stages  of  system  development  and  were  developed  to  augment  established  Air  Force 
procedures.  One  of  the  more  important  applications  of  the  models  would  be  in 
providing  information-  for  the  establishment  of  meaningful  reliability  requirements 
during  the  conceptual  and  validation  phases.  Another  important  application  would 
be  in  employing  the  DSPC  model  during  reliability  validation  tests.  AFR  80-5 
makes  provisions  for  reliability  evaluation  tests,  i.e.,  tests  to  determine  reli¬ 
ability  deficiencies  rather  than  to  demonstrate  achievement  of  specified  values. 
After  identifying  the  reliability  deficiencies  in  a  given  subsystem,  there  are 
various  options  available  for  taking  corrective  action  such  as  redesign,  use  of 
higher  quality  components,  redundancy,  environmental  protection,  etc.  Each  of 
these  options  will  have  associated  with  it  a  certain  reliability  improvement  along 
with  the  cost  of  achieving  his  improvement.  The  DPSC  model  applied  to  this 
subsystem  would  identify  those  corrective  action  options  which  would  provide  the 
highest  performance  at  a  prescribed  cost. 

In  conclusion,  it  is  useful  to  review  briefly  the  stepwise  procedures  and 
inputs  required  for  implementing  the  OAS  analytical  models.  These  procedures  and 
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inputs  are  as  follows: 


Specify  the  mission  profile  by  phases  and  the  subsystem  operating  time 
during  each  phase. 

Identify  the  mission  critical  subsystems  and  specify  their  MTBFs. 

From  failure  modes  effects  analysis  or  other  data  determine  the  condi¬ 
tional  probability  of  abort  given  failure. 

With  the  above  data,  the  mission  completion  success  probability  can  be  calculated 
and  the  subsystems  ranked  according  to  their  probability  of  causing  a  mission 
abort.  In  addition,  a  sensitivity  analysis  can  be  performed  to  determine  the 
increase  in  MCSP  due  to  increasing  the  MTBF  of  any  selected  subsystem.  Even 
without  cost  data  the  above  information  is  useful  for  the  planner  early  in  the 
program  in  identifying  the  most  troublesome  subsystems  and  indicating  those  sub¬ 
systems  for  which  additional  options  are  desired. 

When  options  for  the  various  subsystems  are  available  and  the  acquisition  cost 
of  each  subsystem  option  is  estimated,  the  OAS  model  can  optimize  system  perfor¬ 
mance  over  acquisition  cost,  i.e.,  for  any  level  of  system  acquisition  cost  those 
options  are  identified  which  will  yield  maximum  performance. 

The  next  step  is  to  obtain  the  average  cost  per  repair  for  each  subsystem 
option.  Then  the  model  can  optimize  over  total  system  cost  (acquisition  plus 
logistic  support  cost). 

Finally,  an  appropriate  measure  of  system  effectiveness  must  be  established 
and  the  relationship  between  the  subsystems  and  system  effectiveness  must  be 
determined.  Then  the  model  will  identify  those  subsystem  options  which  yield 
maximum  system  effectiveness  for  any  level  of  total  system  cost. 

The  models  and  methodology  presented  herein  are  just  one  approach  to  providing 
the  decision-maker  with  important  information.  These  models  can  be  extended  if 
more  detailed  analysis  is  required,  and  it  is  hoped  that  this  methodology  will 
provide  some  guidelines  for  other  workers  in  developing  and  formulating  models 
for  their  own  particular  applications. 
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SECTION  I 
INTRODUCTION 


1 .  GENERAL 

During  the  past  few  years,  the  Office  of  the  Assistant  for  Study  Support  (OAS) 
has  been  engaged  in  a  variety  of  analyses  concerning  test  and  evaluation  and  life 
cycle  costs  (References  1  and  2),  and  in  the  course  of  these  analyses,  various 
mathematical  models  were  developed.  An  overview  of  this  approach  to  life  cycle 
cost  and  test  and  evaluation  analysis  is  presented  in  Reference  2. 

For  convenience  of  reference  and  In  the  belief  that  the  models  might  be  of 
use  to  other  workers  in  the  areas  of  test  and  evaluation  and  life  cycle  costs, 
the  models  are  presented  here  along  with  detailed  methodology  and  examples. 

It  should  be  noted  that  a  complete  life  cycle  cost  model  was  not  developed 
during  the  course  of  this  study,  but  rather  such  things  as  subsystem  reliability 
levels  and  logistic  support  cost  and  their  Impact  on  life  cycle  cost  were 
analyzed.  In  addition,  a  generalized  approach  for  relating  system  effectiveness 
to  system  life  cycle  cost  is  developed. 

2.  BACKGROUND 

a.  Operational  Reliability  Deficiencies.  The  principal  Impetus  in  developing 
the  models  came  from  one  of  the  major  findings  in  the  initial  study  performed  by 
OAS  on  test  and  evaluation  analysis.  This  finding  concerned  the  generally  poor 
field  reliability  experienced  by  the  A-7D  avionics  systems.  Low  operational  reli¬ 
ability  of  sophisticated  avionics  equipment  Is  not  in  itself  an  unexpected  revela¬ 
tion,  but  the  low  reliabilities  In  conjunction  with  wide  discrepancies  between  the 
established  reliability  requirements  for  the  A-7D  avionics  systems  and  their 
respective  operational  reliability  levels  does  appear  to  be  significant.  Further 
investigation  showed  that  such  discrepancies  are  not  unique  to  the  A-7D  program 
but  are  also  prevalent  in  other  Air  Force  weapon  systems. 
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b .  Impact  of  Reliability  Deficiencies  on  System  Effectiveness  and  Life  Cycle 
Costs.  In  probing  deeper  into  reliability  problems  associated  with  Air  Force 
avionics  systems,  one  must  examine  the  regulations,  tests,  and  demonstrations 
that  form  the  basis  for  Air  Force  acceptance  and  procurement  of  avionics  systems. 

Air  Force  Regulation  80-5  states  that  "realistic  and  meaningful  R&M  characteris¬ 
tics  and  levels  will  be  determined  by  cost  effectiveness  analysis,  reflecting 
both  system  effectiveness  and  life  cycle  costs."  AFR  80-5  states  further  that 
"the  establishment  of  minimum  acceptable  R&M  levels  must  be  determined  on  the 
basis  of  realistic  operational  needs."  Unfortunately,  there  has  not  been  an 
effective  technique  for  selecting  levels  of  subsystem  R&M  which  would  optimize 
the  total  cost  effectiveness  of  a  system.  Therefore,  the  desired  system  effec¬ 
tiveness  is  not  being  achieved  (Figure  1),  ^n^^fr'Vr^tr-r-iirtr  (f,^yrr  2)  are 
greatly  exceeding  desired  levels. 

In  Figure  1,  the  A-7D  mission  completion  success  probability  (MCSP)  is 
shown  as  a  function  of  mean  time  between  failure  (MTBF)  for  the  A-7D  forward 
looking  radar  (FLR).  MCSP  is  a  measure  of  overall  system  reliability  from  the 
mission  success  standpoint,  and  in  the  results  depicted  in  Figure  1,  the  MTBFs  of 
all  other  subsystems  are  held  constant  at  their  operational  values  while  the  FLR 
MTBF  is  varied  as  shown.  As  shown  in  Figure  1,  there  is  a  wide  discrepancy 
between  the  MCSP  values  corresponding  to  the  operational  MTBF  and  the  MTBF 
requirement  demonstrated  by  MIL-STD-781B  reliability  qualification  testing.  An 
exact  correspondence  is  not  to  be  expected  between  the  operational  MTBF  and  the 
laboratory  demonstration  because  of  differing  environments  and  various  other  factors. 
However,  one  of  the  purposes  of  MIL-STD-781B  testing  is  "facilitating  the  deter¬ 
mination  of  more  realistic  correlation  factors  between  test  reliability  and 
operational  reliability."  The  MCSP  value  corresponding  to  the  operational  MTBF 
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Figure  2 


Logistic  Support  Costs  as  a  Function  of  MTBF  for  the  A-7D 
Forward  Looking  Radar. 
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predicted  on  the  basis  of  reliability  testing  is  also  shown  in  Figure  1.  Not 
only  do  the  operational  and  predicted  MTBFs  differ  significantly,  but  in  addi¬ 
tion,  predicted  system  performance  as  measured  by  MCSP  is  not  being  achieved. 

In  Figure  2,  the  10-year  logistic  support  costs  for  the  A-7D  FLR  are  shown 
as  a  function  of  FLR  MTBF.  Although  logistic  support  costs  are  just  one  part  of 
total  life  cycle  costs,  for  low  reliability  systems  the  logistic  support  cost  can 
be  a  very  significant  part  of  total  life  cycle  costs.  Figure  2  shows  the  logistic 
support  costs  associated  with  the  various  MTBFs  of  interest.  Two  important  points 
should  be  noted  in  Figure  2.  Firstly,  low  MTBF  values  result  in  inordinate 
support  costs,  and  secondly,  failure  to  correlate  reliability  test  results  to 
operational  levels  can  cause  logistic  support  costs  to  be  underestimated  by 
millions  of  dollars. 

3.  OVERVIEW 

In  the  subsequent  sections,  the  details  of  the  models  are  presented  along  with 
examples  of  their  application.  It  is  anticipated  that  these  models  will  have  a 
variety  of  uses  as  management  tools  in  the  systems  acquisition  process. 

The  first  model  to  be  considered  is  the  MCSP  model.  An  MCSP  model  determines 
the  probability  that  the  system  completes  its  mission  without  experiencing  an 
abort  causing  failure.  (With  the  proper  data  input  and  interpretation  of  results 
the  model  can  also  determine  the  probability  that  the  system  completes  its  mission 
with  degraded  effectiveness,  i.e.,  the  mission  is  not  aborted  by  the  failure  but 
full  system  capability  is  not  available.)  MCSP  models  are  quite  useful  during 
test  programs  since  they  provide  a  continuous,  single,  easily  comprehensible 
measure  of  reliability  growth  for  the  total  system.  They  also  highlight  any 
problem  areas  early  in  the  program  so  that  appropriate  action  can  be  taken.  In 
addition,  the  results  of  MCSP  modeling  techniques  can  provide  the  potential  user 
with  early  insight  into  the  operational  suitability  of  the  system  from  the  reli¬ 
ability  standpoint. 
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MCSP  models  by  themselves  are  inadequate  for  life  cycle  cost  analyses  since 
they  do  not  consider  the  cost  of  reliability  development/ improvement,  nor  do  they 
consider  logistic  support  costs.  For  example,  it  is  of  little  value  to  determine 
that  improving  the  reliability  of  a  critical  subsystem  leads  to  dramatic  enhance¬ 
ment  of  the  MCSP  if  the  cost  ramifications  associated  with  the  improvement  are  not 
carefully  considered.  It  could  happen  that  the  cost  of  reliability  improvement 
is  exorbitant  and  exceeds  any  expected  savings  in  logistic  support  costs.  On 
the  other  hand,  by  selecting  subsystems  for  reliability  improvement  based  on 
MCSP,  the  cost  of  reliability  improvement,  and  logistic  support  cost  considerations, 
the  reliability  of  the  total  system  can  be  improved  in  an  optimum  manner.  The 
next  step  in  developing  the  methodology  is  to  consider  reliability  management, 
i.e.,  the  tradeoff  between  reliability  development/improvement  costs  and  logistic 
support  cost  savings.  Optimum  reliability  levels  can  be  selected  in  this  way. 

Combining  reliability  optimization  with  MCSP  considerations  leads  to  the 
development  of  the  Designing  to  System  Performance/Cost  (DSPC)  model.  The  DSPC 
methodology  represents  a  new  and  innovative  approach  to  system  acquisition,  and 
preliminary  results  indicate  that  this  technique  will  provide  very  valuable 
information  to  the  decision-maker.  This  methodology  systematically  identifies 
those  subsystem  options  which  provide  the  highest  system  performance  at  any 
prescribed  level  of  cost  (either  acquisition  cost  or  acquisition  plus  logistic 
support  cost).  The  DSPC  model  is  compatible  with  designing  to  system  cost,  or 
performance,  or  both.  Once  total  system  reliability  specifications  are  established 
each  individual  subsystem  has  a  corresponding  installed  reliability  and  cost  goal 
which  allows  realistic  and  continuous  evaluation  and  adjustment  as  the  subsystem 
is  developed  to  maturity. 

Two  important  applications  of  the  DSPC  model  are  in  establishing  reliability 
requirements  and  reliability  testing.  The  model  would  determine  the  most 
realistic  reliability  levels  for  the  available  funding,  and  would  also  measure 
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the  cost  consequences  and  impact  on  system  performance  if  higher  reliability 
levels  are  desired.  When  applied  at  the  subsystem  level  during  reliability 
testing,  the  model  would  determine  the  most  cost  effective  technique  for  correc¬ 
ting  reliability  deficiencies. 

Finally,  a  generalized  approach  for  combining  system  effectiveness  with  the 
results  of  the  DSPC  methodology  is  presented.  '  The  input  data  required  for  this 
step  is  a  valid  measure  of  effectiveness  for  the  system  under  consideration.  As 
examples,  two  measures  of  effectiveness  for  fighter  aircraft  are  developed.  In 
the  case  of  air-to-ground  fighters,  it  is  shown  that  an  evaluation  of  the  effec¬ 
tiveness  must  account  for  the  interaction  of  availability,  abort  probability, 
kill  potential,  and  survivability;  and  survivability  is  often  the  most  dominant 
factor.  For  air-to-air  fighters,  the  exchange  ratio  (Red  aircraft  destroyed  per 
Blue  aircraft  destroyed)  is  an  important  measure  of  worth,  and  it  can  be  expressed 
as  a  function  of  weapon  effectiveness,  maneuver  capability,  and  first  shot 
probability  with  first  shot  probability  being  the  most  Important  parameter. 

The  Appendix  Section  contains  descriptions  and  listings  for  the  computer 
programs  developed  in  the  study. 
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SECTION  II 


GENERALIZED  MISSION  COMPLETION  SUCCESS  PROBABILITY  MODEL 

1.  INTRODUCTION 

This  section  presents  the  development  of  a  generalized  MCSP  model.  Such 
a  model  can  be  applied  to  any  system  which  can  be  divided  into  mission 
critical  subsystems  for  which  mean  time  between  failure  (MTBF)  data  either 
exists  or  can  be  estimated,  and  for  which  a  mission  profile  can  be  defined. 

OAS  experience  to  date  is  only  with  aircraft  systems.  Therefore,  the  examples 
and  terminology  presented  in  this  report  are  aircraft  oriented.  A  digital 
computer  program  listing  for  the  model  is  presented  in  Appendix  A. 

2.  MISSION  PROFILE 

The  mission  profile  should  be  typical  or  representative  for  the  given 
system.  In  addition,  the  profile  should  be  divided  into  phases  and  the 
subsystems  critical  to  each  phase  should  be  identified.  Figure  3  is  an 
example  of  a  close  air  support  mission  profile  for  the  A-7D. 

3.  MATHEMATICAL  DEVELOPMENT 

a.  Basic  MCSP  Model.  The  MCSP  model  is  based  on  subsystem  failures  which 
follow  the  exponential  distribution.  This  distribution  is  characterized  by  a 
constant  failure  rate  which  is  usually  a  valid  assumption  for  most  of  the 
subsystems  of  interest.  Physically,  a  constant  failure  rate  indicates  the 
subsystem  has  gone  through  a  burn-in  period  so  that  failures  due  to  design 
deficiencies  are  negligible,  and  also  subsystem  components  are  repaired  or 
replaced  on  a  regular  basis  so  that  physical  wearout  does  not  cause  the 
failure  rate  to  increase  with  time.  For  subsystems  for  which  the  exponential 
distribution  is  not  applicable,  the  mathematical  formulation  presented  below 
remains  the  same  with  the  appropriate  distribution  being  utilized,  and  the 
equations  changed  accordingly. 
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Figure  3.  A-7D  Mission  Profile. 
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It  is  first  assumed  that  no  subsystem  is  redundant.  The  modifications 

« 

required  to  account  for  redundancy  are  described  later.  For  an  exponentially 
distributed  failure  pattern,  the  probability,  P-c,  that  the  i-th  subsystem 
completes  its  function  without  a  failure  is  given  by 

np 

E* 

j  =  i 

where  is  the  mean  operating  time  between  failures  for  the  i-th  subsystem, 
t^  is  the  time  the  i-th  subsystem  is  used  in  the  j-th  phase,  and  Np  is  the 
number  of  phases  in  the  mission. 

Before  proceeding  further  with  the  mathematical  development  of  the 
MCSP  model,  some  discussion  of  failure  types  is  warranted.  This  is  important 
for  the  proper  understanding  of  the  very  important  concept  of  probability  of 
abort  causing  failure.  The  importance  of  this  concept  lies  in  the  fact  that 
not  all  subsystem  failures  cause  aborts,  and  furthermore,  not  all  failures 
which  would  normally  cause  an  abort  (depending  during  which  phase  they  occurred) 
reduce  mission  effectiveness  to  zero. 

Failure  types  are  defined  according  to  their  effect  on  completing  the 
mission  and  the  mission  effectiveness.  Failure  types  are  further  categorized 
by  the  mission  phase  during  which  they  occur.  Table  I  lists  an  example  of  the 
various  failure  types  and  their  effects  on  mission  completion  and  mission 
effectiveness. 

After  determing  the  effect  on  the  mission  of  various  failure  types,  the 

I 

next  step  is  to  determine  the  probability  that  the  i-th  subsystem  completes  its  , 
function  without  an  abort  causing  failure.  The  most  convenient  way  to  treat 

i 
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Table  I 


FAILURE  TYPES  AND  THEIR  EFFECTS  ON 
MISSION  COMPLETION  AND  MISSION  EFFECTIVENESS 


TYPE 

EFFulT  ON  MISSION 

MISSION  EFFECTIVENESS 

Before 

Reduction  in  Safety 

Aborted 

None 

Mission 

Objective 

Phase 

Extreme  Reduction 
in  Capability 

Aborted 

None 

Reduction  in 
Capability 

Completed 

Reduced 

Minor  Malfunctions 

Completed 

Full  Capability 

During 

Reduction  In  Safety 

Aborted 

Reduced 

Mission 

Objective 

Phase 

Extreme  Reduction 
in  Capability 

Completed 

Reduced 

Reduction  in 
Capability 

Completed 

Reduced 

Minor  Malfunctions 

Completed 

Full  Capability 

After 

Reduction  in  Safety 

Completed 

Full  Capability 

Mission 

Objective 

Phase 

Extreme  Reduction 
in  Capability 

Completed 

Full  Capability 

Reduction  in 
Capability 

Completed 

Full  Capability 

Minor  Malfunctions 

Completed 

Full  Capability 
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this  is  to  introduce  the  concept  of  mean  operating  time  between  abort  causing 
fa i lures.  For  this  discussion  let  r  denote  the  mean  operating  time  between 
failures  for  a  certain  subsystem.  This  means  that  the  various  failure  modes 
have  been  defined  for  that  subsystem.  Suppose  there  are  n  different  failure 
modes  possible  for  the  subsystem,  i.e.,  in  the  determination  of  t  each  failure 
had  to  be  classified  as  one  of  these  n  modes.  Let  the  n  failure  modes  be 
denoted  by 


frf2 . fn  •  ( 1 1-2) 

Suppose  the  first  k  failure  modes  are  abort  causing  failures.  Given  that  a 
failure  occurs,  the  relative  frequency  of  occurrence  of  the  first  k  failure 
modes  is  some  number  P  .  This  value  P,  is  the  probability  of  abort  given  a 

a  a 

failure  of  the  subsystem.  If  the  subsystem  operates  for  time  T  then  the 
expected  number  of  failures  is 


7  •  (II-3) 

Since  P  is  the  fraction  of  failures  causing  an  abort,  the  expected  number  of 

a 

abort  causing  failures  is 

pa  t  '  (II-4) 

Letting  t,  denote  the  mean  operating  time  between  abort  causing  failures,  then 

a 

it  follows  from  the  definition  of  t,  that  the  expected  number  of  abort  causing 

a 

failures  is  given  by 

—  =  Pa  -  .  ( 1 1 -5 ) 

tax 

a 
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From  equation  ( 1 1 -5 )  the  relation  between  ra  and  t  is  found  to  be 

d 

T  =  jh  •  ( 1 1-6 ) 

a  pa 

Clearly,  the  value  of  Pa  is  dependent  upon  the  definition  of  failure  used  in 

a 

the  determination  of  x.  It  also  depends  on  how  critical  the  subsystem  is 
during  the  j-th  phase. 

To  clarify  these  ideas,  suppose  a  subsystem  has  6  failure  modes  where 
the  consequence  and  effect  on  the  mission  for  each  failure  mode  is  given  in 
Table  II. 

Table  II 

FAILURE  MODES  AND  THEIR  EFFECT  ON  MISSION 


FAILURE 

MODES 

CONSEQUENCE 

EFFECT  ON 
MISSION 

fi 

Extreme  Reduction  in  Effectiveness 

Abort 

f2 

Reduction  in  Safety 

Abort 

f3 

Reduction  in  Safety 

Abort 

f4 

Reduced  Effectiveness 

Continue  with 
Reduced  Eff. 

f5 

Reduced  Effectiveness 

Continue  with 
Reduced  Eff. 

f6 

Minor  Repairs  Required 

None 

Since  only  failure  modes  f -j ,  fg,  and  fj  cause  an  abort,  it  follows  that 

,  (II-7a) 


,  _  P{  fl»  f2»  f3  } 

a  P  |  failure  j 
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where  P  { f -j ,  f g »  f3|  denotes  the  probability  that  failure  mode  or  or 
f3  occurs.  The  probability  of  abort  due  to  safety  factors  (given  that  the 
subsystem  fails)  is 


P 


as 


r\ 

|  f2’  f3  } 

P1 

I  failure 

1 

(II-7b) 


The  probability  of  reduced  (or  zero)  effectiveness  (given  a  failure)  is 


P 


E 


f2»  f3» 
f  failure 


(II-7c) 


The  point  to  be  made  is  that  such  factors  as  Pgs  and  PE  can  be  used  in  the  same 
manner  as  Pa  to  calculate  other  measures,  for  instance,  the  probability  of 

d 

completing  the  mission  without  a  safety  abort  or  the  probability  of  completing 
the  mission  with  maximum  effectiveness. 

The  failure  modes,  the  associated  failure  rates,  and  the  impact  on 
mission  performance  can  be  estimated  for  a  new  subsystem  design  by  component 
analysis,  initial  testing,  or  from  Air  Force  Logistics  Command  data  for  similar 
systems.  As  the  development  of  the  subsystem  progresses  these  estimates  can 
be  updated. 

Using  the  concept  of  mean  operating  time  between  abort  causing  failures, 
the  probability  that  the  i-th  subsystem  does  not  cause  an  abort  (given  that  the 
mission  was  not  aborted  due  to  other  causes)  is  given  by 


fc1  jPai j 


j  =  1 
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( 1 1-8) 


where  P^j  is  the  probability  of  mission  abort  given  that  the  i-th  subsystem 

fails  during  the  j-th  phase.  The  abort  probability  P  . .  depends  upon  how 

a  i  j 

mission  critical  the  1-th  subsystem  is  during  the  j-th  phase.  Since  it  is 
the  relative  frequency  of  failures  which  are  abort  causing  failures  (abort 
type  failures),  P  ..  is  also  dependent  upon  the  definition  of  a  failure.  In 

al  J 

most  cases  a  failure  is  considered  an  abort  type  failure  for  reasons  of  safety 
or  reduced  effectiveness. 

By  calculating  P.  for  each  subsystem,  the  subsystems  can  be  ranked 
according  to  their  likelihood  of  aborting  the  mission.  An  example  of  this 
aspect  of  the  MCSP  model  is  presented  below. 

The  next  item  of  interest  is  the  probability  P  of  completing  the 

L  JC 

Jt-th  phase  without  an  abort  causing  failure.  (In  order  to  reach  the  Ji-th 
phase  all  previous  phases  must  have  been  completed  without  an  abort  causing 
failure.)  This  probability  is  given  by 


ci 


Ns 

n  exp< 
i  =  1 


.  .P  . 
■J  aiJ 


j  =  1 


(II-9) 


where  N$  is  the  total  number  of  subsystems  and  i  is  the  mission  phase  of 
interest.  The  case  l  =  Np  yields  the  mission  completion  success  probability 


MCSP  =  n  Pi  , 
i  =  1 


(H-10) 


where  P.  is  given  by  equation  ( I 1-8) . 

MCSP  by  cumulative  phases  is  of  interest  because  it  makes  it  possible 
to  examine  the  mission  up  to  and  including  any  phase.  For  example,  in  the 
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case  of  a  single  mission,  abort  causing  failures  occurring  after  the  target 
phase  do  not  affect  mission  effectiveness.  However,  in  the  more  interesting 
cases  involving  repeated  sorties,  failures  occurring  during  all  phases  are 
important  since  they  affect  maintenance  requirements  between  sorties.  An 
important  measure  of  maintenance  requirements  is  the  probability  of  completing 
the  mission  without  any  subsystem  failures.  This  measure  is  obtained  by 
setting  all  abort  probabilities,  pa^j»  equal  to  unity  and  using  equations 
( 1 1-8)  and  (11-10). 

Two  other  items  of  Interest  regarding  MCSP  are  the  probability,  P.., 

■  J 

that  the  i-th  subsystem  causes  an  abort  In  phase  j  given  no  abort  before 

phase  j;  and  the  probability,  P  .,  of  abort  in  phase  j  given  no  abort  before 

apj 

phase  j.  These  probabilities  are  given  respectively  by 


Examples  of  applying  the  methodology  are  presented  in  Table  III  and 
Figure  4.  Table  III  shows  the  critical  subsystem  identification  for  the  A-7D. 
On  the  left  hand  side,  the  eight  A-7D  subsystems  with  the  highest  failure  rates 
during  Category  II  testing  are  shown.  On  the  right  hand  side,  the  ranking  of 
the  eight  A-7D  subsystems  causing  the  greatest  number  of  aborts  during 
Category  II  testing  are  shown.  (The  number  in  parenthesis,  1  -  P^ ,  is  the 
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probability  that  the  subsystem  will  cause  an  abort  during  the  mission.)  An 
examination  of  the  Table  reveals  that  only  the  Forward  Looking  Radar  (FLR) 
preserves  the  same  ranking.  Some  subsystems  appearing  in  the  MTBF  ranking  do 
not  appear  in  the  abort  ranking  and  vice  versa.  This  illustrates  the  fact 
that  MTBF  alone  is  not  a  good  indicator  of  the  effect  a  subsystem  will  have 
on  mission  success. 


Table  III 

CRITICAL  SUBSYSTEM  IDENTIFICATION 
MTBF  RANKING  ABORT  RANKING 


MTBF 

1  -pi 

1. 

Forward  Looking  Radar 

(12  hr) 

1. 

Forward  Looking  Radar 

(.117) 

2. 

Inertial  Measurement 

(31  hr) 

2. 

Navigation  Weapon 

System 

Delivery  Computer 

(.064) 

3. 

Lighting 

(34  hr) 

3. 

Inertial  Measurement 
System 

(.056) 

4. 

Navigation  Weapon 
Delivery  Computer 

(35  hr) 

4. 

M61  Gun 

(-047) 

5. 

M61  Gun 

(38  hr) 

5. 

Tactical  Air 

Navigation 

(.032) 

6. 

Tactical  Air 

Navigation 

(44  hr) 

6. 

Radar  Altimeter 

(.032) 

7. 

Radar  Altimeter 

(44  hr) 

7. 

Head  Up  Display 

(.030) 

8. 

Landing  Gear 

(64  hr) 

8. 

Weapons  Release 

(.027) 

The  abort  ranking  is  dependent  upon  the  length  of  time  the  subsystem 
is  used  during  the  mission,  the  MTBF  of  the  subsystem,  and  the  conditional 
probability  that  the  mission  will  be  aborted  given  that  the  subsystem  fails. 
Thus,  in  general,  the  abort  ranking  does  not  correspond  to  the  MTBF  ranking. 
This  example  illustrates  the  way  the  MCSP  model  can  be  utilized  to  identify 
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those  subsystems  whose  reliability  improvement  most  enhances  probability  of 
mission  completion.  The  next  example  illustrates  the  evaluation  of  those 
subsystems  so  identified.  t 

Figure  4  shows  the  results  of  the  type  of  sensitivity  analysis  that 
can  be  conducted  using  the  MCSP  model.  Starting  with  the  baseline  system, 
the  effect  of  improving  the  reliability  of  single  subsystems  or  combination 
of  subsystems  can  be  analyzed.  The  abort  ranking  in  the  previous  Table 
identified  the  Forward  Looking  Radar  (FLR)  and  the  Navigation  Weapon  Delivery 
Computer  (NWDC)  as  the  two  A-7D  subsystems  having  the  most  impact  on  mission 
success.  Increasing  the  MTBF's  of  the  NWDC  and  the  FLR  results  in  dramatic 
improvements  in  MCSP,  while  increasing  the  MTBF  of  relatively  high  reliability 
subsystems  such  as  the  engine  has  essentially  no  effect  on  MCSP.  However,  it 
does  not  follow  that  the  reliability  of  the  engine  should  not  be  improved 

*4 

since  it  is  possible  that  the  cost  of  improvement  could  be  more  than  compensated 
for  by  the  resultant  savings  in  logistic  support  cost.  These  important  considera 
tions  will  be  discussed  later. 

The  methodology  presented  so  far  can  be  used  to  analyze  a  large  number 

of  systems.  In  the  sections  below,  extensions  of  the  basic  methodology  which 

may  be  of  interest  in  other  applications  are  presented. 

b.  Redundant  Subsystems.  To  achieve  an  increase  in  system  reliability  it 

may  be  necessary  to  introduce  redundant  subsystems  provided,  of  course,  certain 

constraints  such  as  weight  and  volume  can  be  met.  If  a  subsystem  has  redundant 

units  then  in  the  expression  (11-10)  for  MCSP,  the  probability  that  the 

subsystem's  function  is  performed  successfully  must  be  adjusted  to  account 

for  redundancy.  The  purpose  of  this  section  is  to  derive  the  expressions  for 

the  successful  performance  of  a  redundant  subsystem's  function  and  also  the 

associated  logistic  support  cost  resulting  from  redundancy.  Two  types  of 
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redundancy  will  be  considered.  The  first  type  is  Operative  Redundancy  which 
will  here  mean  that  all  redundant  units  operate  simultaneously  and  an  abort 
occurs  only  if  all  units  have  an  abort  causing  failure.  The  second  type  of 
redundancy  is  Standby  Redundancy  meaning  that  only  one  unit  is  operating  and 
a  redundant  unit  will  be  switched  on  only  when  the  operating  unit  has  an  abort 
causing  failure.  It  is  assumed  that  there  are  no  switching  failures. 

Let  n  denote  the  total  number  of  redundant  units  for  a  subsystem,  t  • 
the  MTBF  of  the  i-th  redundant  unit,  and  p  .  the  probability  that  a  failure 
of  the  i-th  redundant  unit  is  an  abort  causing  type. 

Operative  Redundancy  will  be  considered  first.  If  the  subsystem 
function  must  be  performed  for  time  T  during  a  successful  mission  then  the 
probability  that  the  subsystem  will  not  cause  an  abort  is 


Pc(n,  T)  = 


(11-12) 


If  all  redundant  units  have  the  same  characteristics  (i.e.,  p  ■  =  p  , 

a  I  a 

Tri  =  Tr,  i  =  1 ,  2,  . . . ,  n)  then  (11-12)  becomes 


Pc(n,  T)  = 


(11-13) 


In  implementing  the  DSPC  model  described  in  Section  IV  it  is  important  to 
know  the  total  cost  (acquisition  plus  logistic  support  costs)  of  n  redundant 
units.  If  all  units  are  the  same  then  the  acquisition  cost  for  n-th  order 
redundancy  is 

nCa  ,  (11-14) 


where  Ca  is  the  unit  acquisition  cost.  Since  each  unit  operates  the  same  amount 

a 
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of  time  the  logistic  support  cost  for  Operative  Redundancy  is 


(11-15) 


where  Cg  is  the  logistic  support  cost  for  a  single  unit.  Thus,  increasing  the 
reliability  by  going  from  a  single  subsystem  to  n  units  with  Operative 
Redundancy  increases  total  cost  by  a  factor  of  n.  It  will  be  shown  that  this 
is  not  true  for  Standby  Redundancy. 

It  is  now  assumed  that  a  subsystem  is  Standby  Redundant  and  that  all 
units  have  the  same  characteristics  (i.e.,  each  unit  characterized  by  p,  and 

a 

ip).  Let  Tgr  denote  the  mean  time  between  abort  type  failures,  i.e.. 


T 


ar 


(11-16) 


With  Standby  Redundancy  only  one  unit  is  operating  and  a  redundant  unit  is 
switched  on  only  when  the  operating  unit  has  an  abort  type  failure.  The 
probability  that  the  subsystem  will  cause  an  abort  will  first  be  determined 
for  the  case  of  two  units  (primary  and  backup  unit). 

The  probability  that  the  primary  unit  fails  during  the  small  time 
interval  (t,  t  +  At)  is  the  product  of  the  probability  that  the  unit  operates 
successfully  for  time  t  multiplied  by  the  probability  it  fails  during  the  next 
time  interval  At,  i.e., 


If  the  primary  unit  fails  at  time  t  the  backup  unit  is  switched  on  to  operate 
the  remaining  time  T-t.  The  probability  that  the  backup  unit  fails  during 
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time  T-t  is 


1  -  exp 


(11-18) 


Therefore,  the  probability  that  both  units  fail  is  obtained  by  summing 
(integrating)  the  product  of  (11-17)  and  (11-18)  over  all  possible  failure 
times  of  the  primary  unit.  Thus, 


1  -  P 


:(2.  T)  ex, p(-  J-)j, 


dJL 

ar 


(11-19) 


Integration  of  (11-19)  yields 


P  (2,  T)  =  exp(-  -4(1  +  -4 

ar  ar 


(11-20) 


In  general 


1  -  Pc(j,  T)  =  f  exp(-  -S-jjl  -  Pfj  -  1,  T  -  t) 
J  Tar  c 


dt 

Tar 


(II-21) 


Repeated  application  of  (11-21)  yields 


n  -  1 


Pc(n.  T)  = 


i  =  0 


(11-22) 


In  other  words,  Pc(n,  T)  is  equal  to  Pc ( 1 ,  T)  multiplied  by  the  first  n  terms 
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in  the  expansion  of  exp 


Observe, 


Lim  Pc(n,  T)  =  exp/-^  W(~-)  -  1. 
n  -*■  "  \  ar  /  \  a r/ 


(11-23) 


Probably  the  most  common  redundancy  is  with  two  units,  i.e.,  n  =  2.  To  calcu¬ 
late  the  logistic  support  cost  of  the  backup  unit  it  is  necessary  to  know  its 
operating  time.  During  a  successful  mission  the  subsystem  operates  for  time 
T;  the  average  operating  time  of  the  second  unit  is 


Dividing  by  T,  the  average  fraction  of  the  to|al  time  T  the  second  unit  is 
operating  is  then 


T 


(11-25) 


Thus,  the  logistic  support  cost  for  the  second  unit  can  be  approximated  by 


T 


(11-26) 


This  means  that  the  change  from  one  unit  to  two  units  increases  the  subsystem 
reliability  from 

PC(1>T)  =  exp  (-7-)  (II-27a) 
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to 


( 1 1 -27b ) 


pc(2.  r) 


Furthermore,  the  corresponding  increase  in  total  unit  cost  is 


AC 


(11-28) 


The  increase  in  cost  for  Standby  Redundancy  can  be  considerably  less  than  that 
with  Operative  Redundancy.  In  Section  IV  some  of  the  above  redundancy  equations 
will  be  utilized  in  an  example  to  show  how  the  DSPC  model  will  identify  condi¬ 
tions  under  which  redundancy  is  the  optimal  choice. 

Next  to  be  considered  is  the  case  of  two  units  with  Standby  Redundancy 
where  the  MTBF's  of  the  primary  and  backup  units  are  different.  Letting  t  , 
and  Tfl2  denote  the  mean  system  operating  time  between  abort  causing  failures 
of  the  primary  and  backup  units,  respectively,  the  expression  for  P  (2,  T) 

V 

becomes 


exp 


For  n  units  with  Standby  Redundancy,  all  MTBF's  being  different,  it  can  be 
shown  that 
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where  r  .  denotes  the  mean  operating  time  between  an  abort  causing  failure  of 

a  i 

the  i-th  redundant  subsystem. 

It  is  instructive  to  Illuminate  some  of  the  above  ideas  by  means  of 
a  s.mple  example  of  a  subsystem  which  operates  3  hours  during  a  mission, 

Pa  =  0.8,  xr  =  10  hrs,  =  $1 M ,  and  C$  =  $3M.  Table  IV  shows  the  increase  in 
performance  and  cost  resulting  from  redundancy.  It  also  shows  the  advantage 
of  Standby  Redundancy,  i.e.,  Operative  Redundancy  is  equivalent  to  increasing 
the  subsystem  MTBF  from  10  hours  to  55  hours  with  total  cost  increasing  from 
$4M  to  $8M,  whereas  Standby  Redundancy  is  equivalent  to  an  MTBF  of  96  hours 
at  a  cost  of  $5.4M. 


Table  IV 

REDUNDANCY  EXAMPLE 


NUMBER 

OF 

UNITS 

TYPE 

REDUNDANCY 

MCSP 

MTBF 

EQUIVALENT 

(hrs) 

ACQ. 

COST 

($M) 

15-YEAR 

LSC 

($M) 

n 

1 

None 

.7927 

10 

1 

3 

4 

2 

Operative 

.9570 

.55 

2 

6 

8 

2 

Standby 

.9754 

96 

2 

3.4 

5.4 

4.  SUMMARY 

In  this  section  the  basic  MCSP  methodology  has  been  presented  along  with 
clarifying  examples  and  extensions  of  the  basic  model  to  include  redundant 
subsystems. 

The  MCSP  model  can  be  used  to  assess  the  reliability  of  the  total  system 
based  on  the  reliability  of  the  individual  subsystems;  rank  the  subsystems  in 
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terms  of  the  probability  of  abort  causing  failures;  and  determine  the  MCSP 
enhancement  due  to  improvements  in  individual  subsystem  reliability. 

In  the  next  section  reliability  management  techniques  are  discussed. 
These  techniques  include  the  cost  considerations  that  must  be  combined 
with  the  MCSP  results  in  order  to  extend  the  methodology  for  applications 
to  life  cycle  cost  analyses. 
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SECTION  III 


RELIABILITY  MANAGEMENT 


1.  INTRODUCTION 

MCSP  models  are  quite  useful  for  identifying  critical  subsystems  and  also 
for  determining  the  enhanced  performance  to  be  gained  by  improving  the  reli¬ 
ability  of  the  critical  subsystems.  However,  as  mentioned  previously,  MCSP 
models  by  themselves  are  inadequate  for  life  cycle  cost  analyses  since  they 
do  not  take  into  account  the  cost  of  reliability  development/improvement  or 
logistic  support  costs.  For  example,  Figure  4  in  the  previous  section  showed 
that  improving  the  MTBF  of  the  engine  had  essentially  no  effect  01  MCSP. 
However,  improvements  in  engine  MTBF  could  significantly  decrease  logistic 
support  costs.  On  the  other  hand,  reliability  improvement  of  the  Navigation 
Weapon  Delivery  Computer  and  the  Forward  Looking  Radar  significantly  improves 
MCSP,  but  the  cost  might  be  so  prohibitive  as  to  preclude  reliability 
improvement  for  these  subsystems. 

2.  METHODOLOGY 

a.  Logistic  Support  Cost.  Logistic  support  costs  can  be  conveniently 
analyzed  by  considering  the  average  cost  per  repair  on  a  subsystem  basis. 

The  average  cost  per  repair  for  a  given  subsystem  is  determined  by  dividing 
the  yearly  logistic  support  cost  by  the  number  of  subsystem  failures  during 
the  yearly  period.  This  data  is  compiled  by  the  AFLC  Air  Materiel  Areas  and 
includes  Field  Maintenance  Cost,  Specialized  Repair  Activity  Cost,  Packing 
and  Shipping  Cost,  Condemnation  Cost,  and  Base  Material  Cost.  Once  the 
average  cost  per  repair  is  established  in  this  way,  logistic  support  cost 
projections  can  be  made  for  future  years  as  shown  below.  (For  subystems  not 
in  the  inventory,  estimates  must  be  made  based  on  similar  subsystems  of 
comparable  complexity.) 


39 


The  average  yearly  logistic  support  cost  for  the  1-th  subsystem,  LSC^, 
Is  given  by 


LSC1  =  (CR1)(R1)  ,  (III-l) 


where  CR^  is  the  average  cost  per  repair  for  the  1-th  subsystem,  and  R.  is 
the  expected  number  of  repairs  for  the  i-th  subsystem  during  the  year.  The 
expected  number  of  repairs  can  be  expressed  as 


R 


1 


(II 1-2 ) 


where  T^  is  the  yearly  operating  time  of  the  i-th  subsystem,  and  t*  is  the 
MTBF  of  the  i-th  subsystem. 

Therefore,  the  average  yearly  logistic  support  cost  for  the  i-th 
subsystem  is  given  by 

T, 

LSCi  =  (CRjK^)  ,  ( 1 1 1-3 ) 

and  the  logistic  support  cost  for  y  years  is  given  by 

i 

Tj 

LSCjy  =  (yXCR^X^i)  .  ( 1 11-4) 

The  total  system  logistic  support  cost  is  given  by 

Ns  T 

LSCy  =  y  S  (cri,(T7)  ■  (hi-5) 

i  ■  1 


where  N$  is  the  total  number  of  subsystems. 
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As  an  example,  Table  V  shows  MTBF  and  cost  per  repair  data  for  three 


hypothetical  subsystems. 


Table  V 

HYPOTHETICAL  SUBSYSTEM  MTBF  AND  AVERAGE  COST  PER  REPAIR  DATA 


SUBSYSTEM 

MTBF  (hours) 

AVERAGE  COST  PER  REPAIR  ($) 

A 

15 

100 

8 

300 

200 

C 

475 

1,000 

Assuming  that  the  total  10-year  operating  time  for  the  three  subsystems  is 
3  x  106  hours  (fleet  size  =  500,  average  monthly  operating  time  =  50  hours), 
the  relationship  between  10-year  logistic  support  cost,  MTBF,  and  average  cost 
per  repair  can  be  presented  as  shown  in  Figure  5.  (Any  other  value  of  operating 
time  would  result  in  merely  a  change  of  scale  for  the  ordinate  and  would  not 
change  the  conclusions. )  In  Figure  5,  it  is  immediately  apparent  that: 

(1)  Subsystem  A  needs  MTBF  improvement.  (Even  a  small  increase 
in  MTBF  will  result  in  significant  logistic  support  cost  savings.) 

(2)  Subsystem  B  appears  satisfactory  (relative  to  the  other 

subsystems). 

(3)  With  Subsystem  C  a  reduction  in  the  cost  per  repair  (rather 
than  MTBF  improvement)  could  lead  to  significant  savings. 

Figure  5  shows  the  ramifications  of  reliability  and  cost  per  repair 
on  logistic  support  costs.  Although  it  is  obvious  that  the  way  to  reduce 
logistic  support  cost  is  to  improve  reliability  and/or  reduce  the  cost  of 
repair;  when  this  type  of  analysis  is  applied  to  each  subsystem  it  systemati¬ 
cally  establishes  priorities,  indicates  realistic  goals,  and  allows  for  the 
proper  allocation  of  resources. 
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Figure  5.  Ten-Year  Logistic  Support  Cost  for  Three  Hypothetical 
Subsystems. 
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b.  Reliability  Optimization.  Reliability  optimization,  i.e.,  the 
trade-off  between  levels  of  reliability  and  lifetime  support  cost  to  decrease 
system  life  cycle  cost  is  a  very  important  adjunct  to  MCSP  considerations. 
Reliability  levels  can  be  established  during  development  by  analyzing  subsys¬ 
tem  options  and  their  associated  costs,  or  by  reliability  improvement  programs 
applied  after  the  system  becomes  operational. 

A  hypothetical  example  of  reliability  optimization  is  shown  in 
Figure  6.  In  Figure  6. a,  logistic  support  costs  for  a  given  subsystem  are 
seen  to  decrease  as  a  function  of  MTBF  as  shown  in  Figure  5.  Figure  6.b 
depicts  the  development  or  improvement  costs  associated  with  establishing 
various  levels  of  subsystem  reliability.  The  life  cycle  cost  of  a  subsystem 
is  the  sum  of  the  acquisition  costs,  logistic  support  costs,  and  operating 
costs.  A  cost  reduction  in  any  of  these  areas  leads  to  reduced  life  cycle 
costs.  Reliability  levels  can  be  established  either  during  subsystem 
development  or  through  reliability  improvement  programs  such  that  the  sum  of 
lifetime  logistic  support  costs  and  reliability  development/ improvement  costs 
can  be  minimized.  This  is  shown  in  Figure  6.c  where  lifetime  logistic  support 
costs  and  reliability  development/ improvement  costs  are  combined  as  a  function 
of  subsystem  MTBF.  Thus,  reliability  goals  can  be  selected  which  minimize  the 
sum  of  reliability  improvement  plus  logistic  support  costs.  This  same  procedure 
can  be  utilized  in  conjunction  with  MCSP  models  to  reduce  acquisition  costs. 

This  subject  is  addressed  in  Section  2-d. 

c.  Cost  of  Repair.  Another  approach  to  logistic  support  cost  reduction 
is  through  reducing  the  repair  costs  for  certain  subsystems.  Specific 
methodology  cannot  be  developed  for  systematically  reducing  repair  costs, 
and  the  problem  must  be  dealt  with  on  a  subsystem  by  subsystem  basis.  In 
general,  during  design  and  development  repair  considerations  should  be 
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emphasized  which  make  equipment  easily  accessible  for  inspection  or  removal. 
Also  designs  utilizing  standardized  components,  tools,  and  test  equipment* 
can  significantly  reduce  maintenance  costs.  During  the  B-l  Mockup  Review, 

297  Requests  for  Alteration  (RFA)  were  developed  and  processed.  Of  the  93  RFA 
involving  maintainability,  34  focused  on  accessibility.  Although  precise 
estimates  of  the  potential  savings  from  this  type  of  preliminary  maintain¬ 
ability  assessment  are  not  available,  it  is  clear  that  such  emphasis  on 
minimizing  repair  costs  during  the  early  stages  of  a  program  can  have  a 
significant  impact  on  lifetime  logistic  support  costs. 

After  a  subsystem  becomes  operational  and  experiences  very  high  repair 
costs,  corrective  action  can  sometimes  be  taken  through  Increased  Reliability 
of  Operational  Systems  (IROS)  programs.  IROS  programs  attempt  to  pinpoint 
causes  of  low  reliability  or  high  repair  costs  and  then  make  recommendations 
for  modifying  the  equipment  to  alleviate  these  problems.  An  excellent 
example  of  this  can  be  drawn  from  the  IROS  program  on  the  A-7D  Air  Data 
Computer. 

The  A-7D  Air  Data  Computer  has  experienced  excessive  logistic  support 
costs  due  to  a  water  ingestion  problem  associated  with  the  pitot  static 
system.  Generally,  when  water  gets  into  the  Air  Data  Computer  it  must  be 
returned  to  the  depot  for  overhaul,  and  this  is  the  major  contributor  to  the 
high  logistic  support  cost.  A  modification  program  is  currently  under  way 
to  correct  this  problem.  Figure  7  shows  the  estimated  savings  in  logistic 
support  costs  that  can  be  expected  after  the  A-7D  fleet  is  modified.  By 
solving  the  water  ingestion  problem,  the  average  cost  per  repair  for  the 
Air  Data  Computer  will  be  decreased  significantly.  As  shown  in  Figure  7,  the 
reduction  in  average  cost  per  repair  is  dramatically  more  cost  effective  than 
doubling  the  MTBF  of  the  unmodified  subsystem.  Along  with  the  cost  per  repair 
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Figure  7.  Ten-Year  Loqistic  Support  Costs  as  a  Function  of 
MTBF  for  the  A-7D  Air  Data  Computer. 
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reduction  a  10-20  percent  increase  in  MTBF  is  expected  after  modification. 

d.  MCSP  and  Logistic  Support  Costs.  The  previous  section  demonstrated 
how  logistic  support  cost  savings  can  be  made  after  a  system  becomes  operational. 
By  utilizing  MCSP  models  and  analyzing  logistic  support  costs,  savings  can  also 
be  effected  during  the  acquisition  phase  of  system  procurement.  Section  3.b 
demonstrated  that  logistic  support  cost  can  be  minimized  if  the  optimum 
subsystem  MTBF  can  be  realized.  However,  it  is  not  always  possible  to  design 
to  optimum  levels  of  reliability  because  of  various  constraints  such  as  time 
factors,  limited  funding,  and  technological  barriers.  Therefore,  additional 
guidance  is  required  in  order  to  establish  realistic  reliability  goals  for 
each  subsystem.  MCSP  models  provide  this  guidance. 

In  order  to  obtain  required  system  performance  for  the  least  costs, 
there  should  be  several  options  available  for  each  candidate  subsystem. 

Figure  8  shows  a  hypothetical  example  of  subsystem  reliability  options.  In 
Section  3.b  the  reliability  development/improvement  graph  was  shown  as  a 
continuous  curve.  Actually  such  graphs  would  consist  of  discrete  points 
since  reliability  levels  would  be  established  in  discrete  steps  rather  than 
continuously.  Figure  8  shows  three  options  which  may  represent  the  same 
subsystem  modified  in  two  cases  and  an  entirely  different  subsystem  performing 
the  same  function  in  the  third  case,  or  any  combination  thereof.  The  length 
of  the  lines  for  each  option  represent  the  lower  and  upper  limits  or  ranges  of 
the  expected  MTBF  of  the  subsystem.  With  subsystem  reliability  options 
available,  MCSP  models  and  logistic  support  cost  data  can  be  used  to  select 
the  most  appropriate  option  for  each  subsystem.  This  selection  will  not 
necessarily  be  the  optimum  as  shown  in  Figure  6.c.  For  example,  it  may  not 
be  possible  to  achieve  the  optimum  MTBF  for  a  given  subsystem  because  of  the 
constraints  associated  with  reliability  development/improvement  mentioned 
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Figure  8.  Options  for  Performance/Cost  Determinations. 
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previously.  Also,  if  the  subsystem  is  especially  critical  to  successful 
completion  of  the  mission,  it  may  be  necessary  to  spend  money  nonoptimally 
(right  hand  side  of  curve  in  Figure  6.c)  for  some  subsystems  in  order  to 
achieve  the  required  performance  (MCSP)  while  not  improving  the  reliability  of 
other  less  critical  subsystems  at  all.  The  net  result  is  the  achievement  of  a 
given  level  of  performance  for  the  least  cost. 

The  importance  of  proper  reliability  management  cannot  be  overemphasized. 
Figure  9  shows  the  consequences  of  not  considering  required  system  performance 
and  logistic  support  costs  in  establishing  reliability  goals.  Figure  9. a  shows 
the  mission  completion  success  probability  as  a  function  of  MTBF  for  the  A-7D 
Navigation  Weapon  Delivery  Computer  (NWDC).  All  other  subsystem  MTBFs  have 
been  held  constant  at  Category  II  values.  The  tick  mark  at  the  left  indicates 
the  MTBF  of  the  NWDC  achieved  during  Category  II  testing  (35  hours).  This 
value  has  been  improved  somewhat  since  the  A-7D  has  become  operational,  but 
it  is  still  well  below  the  mature  system  predicted  level  indicated  by  the 
second  tick  mark  (499  hours).  However,  an  examination  of  the  curve  shows 
that  as  far  as  probability  of  mission  completion  (Pmc)  is  concerned  there  is 
no  reason  to  improve  the  MTBF  beyond  about  150  hours.  The  only  other  reason 
for  high  MTBF  requirements  would  be  to  reduce  logistic  support  costs.  Figure  9.b 
shows  the  logistic  support  cost  (LSC)  as  a  function  of  MTBF  for  the  NWDC  where 
tick  marks  are  again  used  to  indicate  the  Category  II  and  mature  system  MTBFs. 

As  shown  on  the  curve,  a  point  of  diminishing  returns  in  LSC  savings  is  reached 
for  MTBFs  greater  than  about  200  hours.  Since  Category  II,  the  mature  system 
MTBF  prediction  for  the  NWDC  has  been  revised  to  250  hours.  This  is  a  much 
more  realistic  value.  Unfortunately,  reliability  development/ improvement  data 
is  not  available  for  the  NWDC.  Such  data  would  complete  the  analysis  of  the 
NWDC  from  the  reliability  management  standpoint.  Even  without  the  reliability 
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development/improvement  data  this  example  demonstrates  the  first  principle  of 
optimization  in  system  development,  i.e.,  don't  buy  or  strive  for  reliability 
levels  that  are  unreasonable  or  unrealistic. 

3.  SUMMARY 

This  section  has  described  a  procedure  that  ensures  obtaining  required 
system  performance  levels  for  minimum  costs.  This  methodology  is  depicted 
graphically  in  Figure  10.  If  the  information  displayed  in  Figure  10  is 
available  for  most  of  the  mission  critical  subsystems  (generally  options  will 
not  be  available  for  every  subsystem),  realistic  goals  can  be  established  and 
options  can  be  selected  such  that  the  required  performance  of  the  overall 
system  is  obtained  for  minimum  cost.  The  major  limitation  in  this  approach 
is  that  curves  such  as  those  displayed  in  Figure  10  must  be  examined  for  each 
subsystem  for  which  they  are  available,  and  it  is  difficult  and  cumbersome  to 
establish  priorities.  This  is  particularly  critical  if  funds  are  limited. 

The  next  section  discusses  a  procedure  that  systematically  and  in  a  step  by 
step  fashion  selects  the  options  that  offer  the  biggest  payoffs  in  terms  of 
higher  performance/lower  costs. 
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Figure  10. 


Graphical  representation  of  the  Designing 
Cost  Methodology. 


to  System  Performance/ 
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SECTION  IV 

DESIGNING  TO  SYSTEM  PERFORMANCE/COST  MODEL 

1 .  INTRODUCTION 

The  previous  sections  have  laid  the  groundwork  for  development  of  the 
Designing  to  System  Performance/Cost  (DSPC)  model.  In  this  section  the 
detailed  DSPC  methodology  Is  presented  along  with  numerical  examples.  A 
digital  computer  program  listing  for  the  model  is  presented  in  Appendix  B. 

2.  STATEMENT  OF  THE  PROBLEM 

Consider  a  system  consisting  of  a  certain  number  of  mission  critical 
subsystems.  For  some  subsystems,  there  are  options,  each  characterized  by  an 
acquisition  cost,  reliability  (MTBF),  and  average  cost  per  repair.  The  objective 
is  to  select  one  option  for  each  subsystem  such  that  a  maximum  value  of  MCSP  is 
achieved  at  a  cost  not  exceeding  some  prescribed  limit.  (Conversely,  the 
problem  can  be  formulated  in  terms  of  achieving  a  prescribed  value  of  MCSP 
for  the  least  cost.)  Cost  can  be  either  acquisition  cost  only  or  the  sum  of 
acquisition  costs  plus  logistic  support  costs  for  y  years,  e.g.,  y  =  10  or 
15  year  logistic  support  cost.  The  methodology  can  also  be  applied  to  existing 
systems  when  it  is  desired  to  optimize  reliability  improvement  programs. 

The  optimization  procedure  developed  in  this  chapter  is  very  simple  and 
easily  implemented.  The  procedure  will  yield  a  curve  (such  as  that  depicted 
in  Figure  11)  consisting  of  straight  line  segments  connecting  vertex  points. 

The  curve  has  the  following  properties: 

a.  Each  vertex  point  represents  the  maximum  MCSP  achievable  at  the 
associated  cost. 

b.  No  combination  of  subsystem  options  will  yield  a  point  above 
the  curve. 

The  vertex  points  represent  optimal  decision  alternatives;  for  each  such  point 
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the  combination  of  subsystem  options  Is  Identified  which  yields  that  optimal 
combination  of  MCSP  and  cost.  It  will  be  shown  that  going  from  one  vertex 
point  to  the  next  involves  a  change  in  only  one  subsystem  option.  Therefore, 
in  a  sense,  intermediate  points  on  the  straight  line  connecting  two  consecutive 
vertex  points  can  be  realized  (on  a  fleet  basis)  by  equipping  only  a  certain 
fraction  of  the  fleet  with  the  new  option. 

3.  REQUIRED  INPUTS 

The  following  notation  is  introduced  to  describe  the  inputs  required  to 
implement  the  DSPC  model  (a  complete  list  is  given  in  Appendix  B): 

N  =  total  number  of  systems. 

N$  =  number  of  mission  critical  subsystems. 

m  =  average  number  of  missions  per  month  per  system. 

y  =  number  of  years  to  be  considered  in  the  calculation  of 
logistic  support  costs. 

tj  n  operating  time  of  i-th  subsystem  (1  =  1,2,  ....  N$)  during 
one  mission,  i.e.,  duty  cycle  of  i-th  subsystem. 

ctj  =  ratio  of  total  operating  time  to  mission  operating  time. 

T.j  =  12ymc^tj  =  total  y-year  operating  time  of  subsystem  i. 

P  . 4  =  probability  that  a  failure  of  the  i-th  subsystem  during  the 
j-th  mission  phase  will  cause  an  abort  of  the  mission. 

n(i)  =  number  of  options  for  the  i-th  subsystem. 

C.  •  =  Cost  of  the  j-th  option  for  the  i-th  subsystem  (j  =  1,  2,  .... 

*  J 

n(i ) ;  i  =  1,  2,  ....  Ns). 

t.  .  =  lower  MTBf  for  the  j-th  option  for  the  i-th  subsystem. 

*  J 

7. .  =  upper  MTBF  for  the  j-th  option  for  the  i-th  subsystem. 

*  J 

CR. .  =  average  cost  per  repair  associated  with  the  j-th  option  for 

'  J 

the  i-th  subsystem. 
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As  shown  in  Section  II,  the  MCSP  is  a  function  of  the  duty  cycle  t . ,  the 
abort  probability  P  ..,  and  the  reliability  (MTBF)  of  each  of  the  Nc  mission 
critical  subsystems.  The  performance/cost  tradeoffs  arise  from  the  different 
options  available  for  the  subsystems,  i.e.,  for  a  unit  acquisition  cost  of 
C. .  dollars  for  the  j-th  option  of  subsystem  i,  the  subsystem  will  have  an 

■  J 

MTBF  of  at  least  hours  and  possibly  as  high  as  7..,  and  the  average  cost 
per  repair  will  be  CR  . ..  If  one  option  is  selected  for  each  subsystem  the 

*  J 

MCSP  is  determined,  and  the  total  y-year  cost  (excluding  operating  costs)  is 
the  sum  of  the  acquisition  costs  plus  the  y-year  logistic  support  costs  of 
the  subsystems.  The  y-year  logistic  support  cost  for  the  j-th  option  for 
subsystem  i  is 


NT. 

— lcr,,  .  (IV-1) 

Tij  J 

Therefore,  the  total  y-year  cost  (excluding  operating  costs)  of  the  j-th  option 
for  subsystem  i  is 


The  options  for  each  subsystem  can  always  be  ordered  in  terms  of  increasing 
MTBF  such  that  t.,  ,  >  t..,  i.e.,  the  reliability  of  the  (j  +  l)st  option 

'  J  *  *  s) 

is  equal  to  or  greater  than  that  of  the  j-th  option.  This  relation  is  assumed 
to  hold  for  each  subsystem.  It  also  should  be  mentioned  that  to  optimize  with 
respect  to  acquisition  cost  only,  the  value  of  y  should  be  set  equal  to  zero. 
4.  DESCRIPTION  OF  OPTIMIZATION  PROCEDURE 

For  clarity  it  is  desirable  to  change  notation  slightly  from  that  in 
Section  II  and  to  express  the  MCSP  function  in  slightly  different  form.  Let 
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t ..  denote  the  MTBF  of  the  i-th  subsystem,  e.g.,  for  i  =  1 ,  2,  . .. ,  Ns,  i .  is 
one  of  the  values  x..,  j  =  0,  1,  ....  n(i).  Once  the  MTBF  of  each  subsystem 

*  J 

is  specified  then  the  MCSP  denoted  by  Pmc  is  given  by 

Ns 

Pmc  ■  D  P,(x,)  .  (IV-3) 

1  ■=  1 


where  P.(x^)  denotes  the  probability  that  the  i-th  subsystem  does  not  have  an 

l 

abort  causing  failure.  Observe  that  if  the  value  of  t.  is  changed  to  x.  then 
the  resulting  MCSP  becomes 


(IV-4) 


Letting 


l 


( IV— 5 ) 


the  incremental  change  in  Pmc  resulting  from  the  MTBF  change  from  x.  to  x.  can 
be  written 


I 


( IV-6) 


Thus,  Pmc  needs  to  be  calculated  only  for  the  baseline  system  (x.  =  x^Q  for 
i  =  1,  2,  ....  N  ) ,  and  any  changes  in  P  resulting  from  the  selection  of  a 

j  IIIV^ 

new  option  can  be  calculated  easily  using  the  above  procedure. 

It  is  clear  that  the  optimization  problem  can  be  formulated  as  a  zero-one 
integer  linear  programming  problem,  i.e.,  letting  x,- .•  =  1  if  the  j-th  option 

*  J 

for  subsystem  i  is  selected  and  0  otherwise,  the  problem  is  (for  some 
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prescribed  cost  constraint  C)  to: 


maximize  log  Pm„ 
3  me 


N$  n(i) 

2  xu 109  pi<Tu) 

i  =  1  j  =  0 


( I V- 7a ) 


subject  to 


i  j 


=  1,  i  -  1,  2, 


( I V- 7b ) 


( IV-7c) 


Although  algorithms  exist  for  solving  such  zero-one  integer  problems,  they 
require  rather  complex  computer  programs.  A  much  simpler  and  straightforward 
optimization  procedure  will  be  developed  which  will  yield  an  optimal  curve 
such  as  that  shown  in  Figure  11. 

To  determine  the  starting  point,  it  is  first  necessary  to  calculate  the 
baseline  MCSP  and  cost: 


mco 


n  Pi(Tio} 
i  =  1 


( IV-8) 
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Figure  11.  Optimal  DSPC  Curve. 
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The  next  step  is  to  check  each  subsystem  for  the  possibility  of  reliability 
optimization  as  described  in  Section  III.  This  means  that  each  subsystem  must 
be  checked  to  determine  if  there  is  an  option  j  >  0  such  that  C".  j  <  r.Q.  If 
an  option  with  this  property  exists,  then  an  MCSP  greater  than  Pmco  can  he 
achieved  at  a  cost  less  than  the  baseline  cost  since  the  MTBF's  for  the  options 
are  ordered,  i.e.,  t..  ,  m-.  Thus,  if  for  some  subsystem  i  there  exists 

'  J  *  '  SJ 

an  option  j  >  0  such  that  1^0»  then  an  adjusted  baseline  (with  a  higher 
MCSP  at  lower  cost)  is  determined  as  follows: 

a.  For  each  subsystem  i  calculate 

min  jej  .  (IV-10) 

0<j<n(i)| 

b.  Let  m(i)  denote  the  maximum  (in  case  there  are  several  minima) 
value  of  j  for  which  C . .  is  a  minimum.  Let  option  m(i)  be  the  adjusted 

'  J 

baseline  for  subsystem  i. 

c.  Reject  the  options  j  *  0,  1,  ...,  m(1)  •  1  since  these  all  result 
in  a  lower  MCSP  at  higher  cost. 

The  adjusted  baseline  for  the  i-th  subsystem  has  an  MTBF  of  at  a  cost 

C i m ( i ) .  The  only  remaining  options  for  the  i-th  subsystem  are 

Tij*  Sj  |  *  m(i )<j<n(i )  .  (IV-11) 

For  the  adjusted  baseline  system  the  new  values  of  MCSP  and  cost  are  calculated 
in  a  manner  analogous  to  that  described  in  equations  (IV-8)  and  (IV-9),  or  by 
repeated  application  of  equation  (IV-4). 

With  the  adjusted  baseline  system  as  the  new  starting  point  and  each 
subsystem  having  the  options  defined  by  (IV-11),  the  optimization  procedure 
can  now  be  developed. 
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A  procedure  to  be  discussed  first  is  the  method  of  steepest  slope  which 
can  be  described  as  follows: 

At  any  vertex  point  the  next  vertex  point  is  determined  by 
selecting,  among  the  remaining  options,  that  option  which  max¬ 
imizes  the  improvement  in  MCSP  per  dollar,  i.e.,  maximizes 


A^mc 

ACost 


( I V- 1 2 ) 


Letting  k(i)  denote  the  option  for  the  i-th  subsystem  at  some  vertex  point 
then  the  next  vertex  point  is  determined  by  finding  the  subsystem  i  and 
option  j  from  the  maximization 


max  max 

l<j<N$  k(i)<j<n(i) 


MTik(i ) *  Tij)  "  1 
^i j  "  ^ik(i) 


(IV-13) 


Notice  that  the  value  of  Pmc  at  the  vertex  point  does  not  appear  in  (IV-13) 
since  it  enters  only  as  a  constant  factor  as  shown  by  equation  (IV-6).  Although 
the  method  of  steepest  slope  works  in  many  cases,  it  fails  to  select  the  optimal 
vertex  point  when  there  exists  a  combination  of  subsystems  and  options  whose 
total  incremental  cost  is  less  than  the  incremental  cost  of  the  selected 
system  and  whose  combined  aP^.  exceeds  that  of  the  selected  systems.  For 
example,  suppose  the  seiecttJ  system  has  A  cost  =  10  and 


(1.2  -  1) 


=  0.2P 


me 


It  follows  that 


aP_ 


me 


ACost 


-  02P. 


me 


(IV-14) 
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Suppose  there  are  options  for  10  other  subsystems  each  with  an  incremental 
cost  of  1.0  and  with  AP|nc  =  Pmc  (1.0199  -  1).  Thus,  for  each  of  the  10 
subsystems 

3St  ■  f-0'")  Pmc  •  <IV-,5> 

which  is  less  than  the  slope  given  by  ( I V- 1 4 ) .  However,  if  all  10  subsystems 
are  selected  their  combined  effect,  by  repeated  application  of  equation 
(IV-6),  is 

AP  =  Pm„  (1 . 01 99  -  1)  =  .218Pm„  ( IV-16) 

me  me  me 

with  a  total  incremental  cost  of  10.  This  shows  that  the  method  of  steepest 
slope  does  not  always  select  the  best  option  since  there  can  exist  a  combination 
of  subsystems  with  smaller  costs  yielding  a  better  result.  However,  if  all 
incremental  costs  were  equal  the  method  would  work.  This  suggests  a  modifica¬ 
tion  of  the  method  of  steepest  slope  which  will  be  described  in  the  next 
paragraph. 

Let  k(i)  denote  the  option  for  the  i-th  subsystem  at  some  vertex  point. 

For  i  =  1 ,  2,  . . . ,  Ns  and  j  =  k(i)  +  1 ,  . . . ,  n(i),  calculate: 

(a)  X,fx1|((i).  Tjj)  .  ( IV-1 7a ) 

(b)  =  l/NjVSkO')}  ■  (IM7b) 

■  (M'ikd)  •  Tij> }  J  •  <nM7c> 

The  next  subsystem  (to  be  replaced)  and  its  option  is  determined  by  selecting 
i  and  j  such  that  is  a  maximum.  The  above  calculations  then  have  to  be 
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repeated  only  for  that  subsystem  and  option  which  was  added  (all  the  other 
a..  values  remain  the  same),  and  the  process  is  continued  by  selecting  the 

’  J 

maximum  a.,  among  the  new  set. 

*  J 

The  procedure  described  above  is  equivalent  to  considering  a  number  aC.. 

'  J 

of  separate  subsystems  each  costing  one  unit  of  cost  and  yielding  a  relative 
change  in  MCSP  of 

- 1  • 

These  aC..  pseudo-subsystems,  each  of  1  unit  of  cost,  have  the  property  that 

'  J 

when  all  aC. •  are  selected  then  the  incremental  change  in  MCSP  is 

*  J 

(  aC.  • 

4Pmc  =  J  -  1  i 

=  Pmc  j  V'ikO)  •  Ti j*  "  1 j  ’  (IV-18) 

and  the  incremental  cost  is  aC^j.  In  other  words,  the  selection  of  all  AC^j 
of  these  pseudo-subsystems  is  equivalent  to  selecting  the  j-th  option  of 
subsystem  i.  It  remains  to  be  shown  that  this  selection  process  is  optimal 
in  the  sense  described  above  in  the  statement  of  the  problem. 

If  the  value  of  MCSP  at  a  vertex  point  is  Pmc  and  if  the  i-th  subsystem 
with  option  j  is  chosen  for  the  next  vertex  point,  then  the  value  of  MCSP  at 
that  vertex  point  is 

<V‘Cijpn,c  •  "V-IS) 

with  an  incremental  cost  aC^..  For  an  incremental  cost  AC  <_  aC^  the  pseudo¬ 
path  between  the  two  consecutive  vertex  points  has  the  value 

AC 

(xij)  Pmc  "  Mx1k(i)  ’  Tij^  Pmc 
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(IV-20) 


It  is  easily  shown  that  the  value  (IV-20)  lies  below  the  straight  line  connect¬ 
ing  the  two  vertex  points.  Furthermore,  the  pseudo-path  between  two  consecutive 
vertex  points  has  monotonically  increasing  slope  and  has  the  form  depicted  by 
the  dashed  curve  in  Figure  12. 

Suppose  the  selection  process  leads  to  n  ordered  values  of  the  A...  Let 

*  J 

these  n  values  of  A.,  and  the  corresponding  incremental  costs  be  denoted  by 

'  J 

A1  >  x2  >  '*•  *  An  *  ( 1 V-21 a ) 

ACj ,  AC2,  ....  ACn  .  ( I V-21 b ) 

For  the  purpose  of  proving  that  the  procedure  is  optimal,  the  assumption  of 
strict  inequality  in  ( I V-21 a )  is  justified.  Assume  that  for  some  cost  C  there 
exists  a  combination  of  subsystem  options  with  a  total  cost  C  and  with  an  MCSP 
above  the  curve  generated  by  the  procedure  (IV-17).  In  terms  of  the  aC^ 
defined  in  (I V-21)  the  cost  C  can  be  written 

C  =  aCj  +  AC2  +  ...  +  ACk  +  r  ,  (IV-22) 


where 


0  <r  <_  aC^  +  -j  and  k  _<  n 

In  other  words,  C  lies  between  the  costs  corresponding  to  the  k-th  and 
(k  +  l)-th  vertex.  In  Figure  13  the  dashed  pseudo-path  leading  from  the  k-th 
to  the  (k  +  l)-th  vertex  is  shown.  In  reaching  the  point  (on  the  pseudo-path) 
corresponding  to  cost  C,  the  greatest  C  values  of  A  were  selected  to  yield 
the  MCSP  of 

AC,  aC,  ACk  r 

pCo(V  (*2>  {Ak  +  l}  *  (IV-23) 
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Cost 

Figure  12.  Pseudo-Path  between  Two  Vertex  Points. 
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Figure  13.  Optimal  P  at  Cost  C. 

illL 
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where  PCq  denotes  the  MCSP  of  the  adjusted  baseline  system.  Any  other  combina¬ 
tion  of  subsystem  options  with  total  cost  C  requires  that  some  of  its  corres¬ 
ponding  a  values  be  different  from  those  in  (IV-23).  Consequently,  for  any 
other  combination  some  of  the  a  values  in  (IV-23)  would  be  replaced  by  smaller 
values  since  (IV-23)  contains  the  C  greatest  a  values.  Thus,  any  other  combina¬ 
tion  of  subsystem  options  with  cost  C  leads  to  an  MCSP  lying  not  only  below  the 
straight  line  segment  but  also  below  the  dashed  pseudo-path.  This  completes 
the  proof. 

This  optimization  procedure  has  been  developed  using  the  higher  confidence 
reliabilities  t..  rather  than  the  upper  limit  7...  Once  the  optimum  curve  has 

I  J  '  J 

been  obtained,  its  upper  limit  can  be  obtained  by  substituting  7..  for  t..  in 

*  J  *  J 

the  appropriate  equations.  However,  this  curve  is  not  necessarily  optimal 
since  the  ordering  and  selection  process  for  the  options  would  in  general  be 
different  when  the  optimization  procedure  is  with  respect  to  the  7. ..  If  a 

'  J 

high  risk  program  is  to  be  considered,  the  model  should  be  exercised  in  both 
ways,  i.e.,  determine  the  upper  limit  of  the  optimum  curve  based  on  high 
confidence  MTBFs  and  optimize  with  respect  to  the  upper  limit  MTBFs. 

5.  NUMERICAL  EXAMPLE 

To  illustrate  the  procedure  described  in  the  previous  section  it  is 

/ 

instructive  to  carry  through  the  calculations  in  detail  for  a  specific  example. 

Consider  one  system  (N  =  1)  consisting  of  3  subsystems  (N$  =  3),  each 
subsystem  having  3  options  (n(l)  =  n(2)  =  n(3)  =  3).  The  cost  to  be  considered 
will  be  the  sum  of  acquisition  costs  plus  15-year  logistic  support  costs 
(y  =  15).  An  average  of  10  missions  per  month  (m  =  10)  is  assumed.  The  duty 
cycle  t.- ,  abort  probability  P  . ,  and  15-year  operating  time  T.  for  each 
subsystem  is  presented  in  Table  VI. 
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Table  VI 


SUBSYSTEM  OPERATING  CHARACTERISTICS 


SUBSYSTEM  (i) 

ti  (hrs) 

Pa1 

BOB 

1 

3 

0.2 

5400 

2 

0.5 

1.0 

900 

3 

2 

0.8 

3600 

The  subsystem  options  are  defined  in  Table  VII. 

Table  VII 
SUBSYSTEM  OPTIONS 


OPTION  0 

(Baseline) 

OPTION  1 

OPTION  2 

i 

Ci0 

Ti0 

Ti0 

CRiO 

Cil 

Til 

~il 

Ci2 

Ti2 

Ti2 

CRi? 

1 

3 

10 

15 

O 

6 

16 

20 

.04 

15 

30 

36 

.05 

1 

5 

8 

.10 

2 

10 

15 

.20 

14 

25 

30 

.20 

2 

8 

10 

.07 

8 

12 

16 

.09 

20 

22 

30 

.15 

The  unit  of  cost  assumed  in  this  example  is  $10,000.  Using  the  costs  and  MTBF 
values  for  each  option  given  by  Table  VII  and  using  the  operating  characteristics 
given  in  Table  VI,  the  values  C-  ■  and  P. (x . - )  are  calculated  using  equations 

I J  •  '  J 

(IV-2)  and 


VTij> 


(IV-24) 
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Those  results  are  presented  in  Table  VI 11. 


Table  VIII 

COST  AND  MISSION  PERFORMANCE  FOR  EACH  SUBSYSTEM  OPTION 


OPTION 

0 

OPTION 

1 

OPTION  2 

i 

Pi^TiO) 

Pi (Til ) 

^11 

pi(Ti2)  ci2 

1 

.9418 

24.6 

.9632 

19.5 

.9802  24.0 

2 

.9048 

19.0 

.9512 

20.0 

.9802  21.2 

3 

.8187 

33.5 

.8752 

35.0 

.9299  44.6 

3 

n  Pi(Ti0)  =  .6976  .  (IV-25) 

1  *  1 

3 

C10  ■  77.1  (IV-26) 

i  =  1 

Checking  each  subsystem  for  reliability  optimization  shows  that  Option  1 
for  Subsystem  1  should  replace  Option  0  since  it  yields  a  higher  MCSP  at  lower 
cost.  In  other  words,  the  increase  in  acquisition  cost  in  going  from  Option  0 
to  Option  1  is  more  than  compensated  for  by  the  savings  in  logistic  support 
cost.  Thus,  the  adjusted  baseline  system  consists  of  Option  1  for  Subsystem  1 
and  Option  0  for  Subsystems  2  and  3.  This  combination  of  options  will  be 


For  the  baseline  system: 


mco 
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denoted  by  (1,  0,  0).  The  MCSP  and  cost  of  the  adjusted  baseline  system  is 


.6976 


.9632 

794T8 


Cost  =  72 


.7135  . 

( IV-27) 


With  the  adjusted  baseline  system  established,  the  optimization  procedure 
described  by  equations  (IV-17)  can  now  be  applied.  Since  Option  1  has  been 
selected  for  Subsystem  1,  the  value  of  k(l)  is  set  equal  to  1.  The  values 
of  k(2)  and  k (3)  are  0.  Table  IX  lists  the  values  of 

Ai  (k(i ))  =  max  /  A. .  I  ,  (IV-28) 

j  >  k(i) 

from  which  the  optimal  options  can  be  determined. 


Table  IX 

EVALUATION  OF  OPTIONS 


SUBSYSTEM  i 

o 

11 

k(i)  =  1 

Aj(0) 

NEXT  ELIGIBLE 
OPTION 

»,-(!) 

NEXT  ELIGIBLE 
OPTION 

1 

-- 

— 

1.0039 

2 

2 

1.0513 

1 

1.0253 

2 

3 

1.0455 

1 

1.0063 

2 

Starting  with  the  combination  of  options  (1,  0,  0)  the  next  vertex  point 
is  determined  from  Table  IX  by  finding  the  maximum  of  x^(k(i))  where 
k ( 1 )  =  1,  k(2)  =  k(3)  =  0.  This  maximum  is  1.0513  which  means  Option  1  for 
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Subsystem  2  should  be  added  to  yield  (1,  1,  0).  The  next  option  is  determined 


from  the  maximum  > .  (k( i ) )  for  k ( 1 )  =  k(2)  =  1 ,  k(3)  =  0.  This  gives  (1 ,  1 ,  1 ) . 
Proceeding  in  this  manner  yields  the  sequence 

(0,  0,  0)  *  (1,  0,  0)  -*  (1,  1,  0)  -  (1,  1,  1)  * 

(1,  2,  1)  -  (1,  2,  2)  -  (2,  2,  2)  .  ( IV-29) 

Using  Table  VIII  to  determine  the  corresponding  Pmc  and  cost  of  each  configuration 
yields  the  results  shown  in  Table  X. 

Table  X 

OPTIMAL  MCSP  AND  COSTS 


CONFIGURATION 

P 

me 

(lower) 

COST 

P 

me 

(upper) 

COST 

(0,  0,  0) 

.6976 

77.1 

.7691 

56.8 

(1,  0,  0) 

.7135 

72.0 

.7768 

56.2 

(1,  1,  o) 

.7501 

73.0 

.7998 

58.0 

(1,  1.  1) 

.8019 

74.5 

.8492 

59.0 

(1,  2,  1) 

.8263 

75.7 

.8635 

65.0 

(1,  2,  2) 

.8779 

85.3 

.9049 

74.8 

(2,  2,  2) 

.8934 

89.8 

.9171 

80.5 

Figure  14  shows  the  optimal  MCSP  vs  cost  curve.  The  curve  corresponding  to 
the  upper  values  of  MCSP  is  not  plotted.  For  this  simple  example,  there  are 

3 

3  =  27  possible  combinations  of  options,  and  for  purposes  of  illustration  all 
combinations  were  calculated  and  are  plotted  in  Figure  14. 


70 


To  demonstrate  the  necessity  of  an  efficient  algorithm  for  calculating  the 
optimal  MCSP  vs  cost  curve,  it  is  instructive  to  discuss  the  case  of  the  A-7D. 
For  this  aircraft  a  total  of  36  mission  critical  subsystems  were  identified. 
If,  for  example,  during  the  planning  phase  there  were  3  options  for  each 
subsystem  then  the  total  number  of  combinations  of  subsystem  options  would  be 

336  =  1.5  x  lo17  .  ( I V- 30 ) 

Even  if  a  computer  required  only  1  millisecond  to  calculate  the  MCSP  and  cost 
associated  with  each  combination,  a  total  computer  time  of  4.8  million  years 
would  be  required  to  compute  all  combinations. 

It  is  easily  shown  that  the  procedure  described  in  this  chapter  requires 
at  most  the  calculation  of 


£  Ml).  1.)  ( IV-31 ) 

i  =  1 

values  of  the  A^.  The  ordering  of  these  values  then  gives  the  optimal 
options.  For  the  above  mentioned  A-7D  example  of  36  subsystems  each  having 
3  options,  the  maximum  number  of  calculations  of  the  a^.  values  is 

36  =  108  .  ( IV-32) 

It  is  instructive  to  apply  the  optimization  procedure  to  the  above  example 
when  the  system  is  optimized  with  respect  to  acquisition  costs  only.  The 
results  of  the  optimization  procedure  lead  to  the  following  sequence  of 
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configurations : 


(0,  0,  0)  -(0,1,0)-  (0,  1,  1)  -  (1,  1,  1)  - 

(1 ,  1 ,  2)  -  (1 ,  2,  2)  -  (2,  2,  2)  .  ( I V- 33 ) 

The  sequence  of  system  configurations  (IV-33),  optimized  for  acquisition  costs, 
differs  significantly  from  the  sequence  (IV-29)  which  is  optimal  for  the  sum  of 
acquisition  and  15-year  logistic  support  costs.  Table  XI  shows  the  acquisition 
cost  and  MCSP  for  the  sequence  of  configurations  (IV-33). 


Table  XI 

OPTIMAL  MCSP  AND  ACQUISITION  COSTS 


CONFIGURATION 

^mc 

(lower) 

ACQUISITION 

COST 

(0,  0,  0) 

.6976 

6 

(0,  1,  0) 

.7334 

7 

(0,  1,  1) 

.7840 

13 

(1.1,1) 

.8019 

16 

(1,  1,  2) 

.8520 

28 

(1,  2,  2) 

.8779 

40 

(2,  2,  2) 

.8934 

49 

It  is  instructive  to  investigate  the  consequences  (in  terms  of  total 
15-year  costs)  of  designing  a  system  to  acquisition  cost.  For  this  purpose 
the  total  15-year  costs  were  calculated  for  the  configurations  identified  in 
Table  XI.  In  Figure  15  the  results  (encircled  points)  are  compared  with  the 
curve  which  resulted  from  optimizing  to  total  15-year  cost.  As  shown  in 
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Acquisition  Plus  15-Year  Logistic  Support  Costs 

15.  Comparison  of  Acquisition  Cost  Optimization  with 
Total  Cost  Optimization. 


the  figure,  most  of  the  encircled  points  lie  far  off  the  optimal  curve.  This 
large  discrepancy  between  the  results  for  the  two  cases  indicates  the  implications 
in  designing  to  acquisition  cost  rather  than  total  15-year  cost. 

To  demonstrate  how  the  optimization  procedure  treats  redundancy  options  it 
is  now  assumed  that  Subsystem  3  can  be  redundant.  The  operating  characteristics 
and  subsystem  options  defined  in  Tables  VI  and  VII  remain  the  same  but  Subsystem  3 
can  also  have  redundancy  with  Option  0,  Option  1,  or  Option  2.  Standby 
redundancy  with  two  identical  units  (primary  and  backup)  is  assumed.  The 
optimization  procedure  will  identify  when  (at  what  cost  level)  redundancy  should 
be  considered  and  also  identify  which  option  should  be  redundant. 

The  values  of  the  cost  and  mission  performance  parameters  listed  in  Table 
VIII  remain  the  same;  however,  the  corresponding  values  must  be  calculated  for 
redundancy  of  Option  0,  Option  1,  and  Option  2  with  Subsystem  3.  Using 
equations  (11-14),  (11-20),  and  (11-26)  these  values  are  calculated  and  listed 
in  Table  XII. 


Table  XII 

COST  AND  MISSION  PERFORMANCE  FOR  STANDBY 
REDUNDANCY  OPTIONS  FOR  SUBSYSTEM  3 


OPTION  0 
REDUNDANT 

OPTION  1 
REDUNDANT 

OPTION  2 
REDUNDANT 

P3R^Ti 0  ^ 

C3R0 

P3R  ^  Ti  1  ^ 

C3R1 

P3R  ^ T i 2  ^ 

C3R2 

.  .9825 

38.7 

.9919 

44.8 

.9975 

65.4 

Starting  with  the  baseline  for  Subsystem  3  (i.e.,  k(3)  =  0)  Table  IX  shows 
that  the  next  eligible  nonredundant  option  for  Subsystem  3  is  Option  1  with 
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A^O)  =  1.0455.  This  value  must  be  compared  with  corresponding  a  values  for 
the  redundancy  options.  Going  from  Option  0  to  Option  OR  (with  redundancy) 
yields  a  a  value  of 

1.0357.  ( I V- 34 ) 

Going  from  Option  0  to  Options  1  or  2  with  redundancy  yields  the  a  values 

1.0171 

and  ( I V- 35 ) 

1.0062  . 

Since  all  A  values  for  the  reliability  options  are  less  than  a3(0)  -  1.0455 
the  next  eligible  option  is  Option  1  without  redundancy.  After  Option  1  is 
selected  Table  IX  shows  that  the  next  eligible  nonredundant  option  is  Option  2 
with  AgO)  =  1.0063.  This  value  must  be  compared  with  the  A  values  associated 
with  going  from  Option  1  to  redundant  Option  0,  1,  and  2.  These  values  are 
1.0322,  1.0128,  and  1.0043.  Thus,  the  next  eligible  Option  is  Option  1  with 
redundancy.  For  the  redundancy  options  of  Subsystem  3  the  values  correspond¬ 
ing  to  those  of  Table  IX  are  given  in  Table  XIII. 

Table  XIII 

EVALUATION  OF  OPTIONS  FOR  STANDBY  REDUNDANCY 
OPTIONS  FOR  SUBSYSTEM  3 


SUBSYSTEM  i 

k(i)  =  0 

k(i)  =  1 

k(i)  =  2 

^(0) 

NEXT  ELIGIBLE 
OPTION 

NEXT  ELIGIBLE 
OPTION 

Ai  (2) 

NEXT  ELIGIBLE 
OPTION 

1 

-- 

- 

1.0039 

2 

-- 

- 

2 

1.0513 

1 

1.0253 

2 

-- 

i 

3 

1.0455 

1 

1.0322 

0  +  Redundancy 

1.0015 

1  +  Redundancy 
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The  sequence  for  selecting  options  is  then 


(0,  0,  0)  >  (1,  0,  0)  ,  (1,  1,  0)  ,  (1,  1,  1)  ,  (1,  1,  OR)  , 

(1,  2,  OR)  ,  (2,  2,  OR)  (2,  2,  1R)  ,  (2,  2,  2R)  .  (IV-36) 

Observe  that  redundancy  was  not  selected  until  late  in  the  sequence.  The 
optimal  MCSP  and  costs  are  presented  in  Table  XIV. 


Table  XIV 

OPTIMAL  MCSP  AND  COSTS  FOR  STANDBY 
REDUNDANCY  OPTIONS  FOR  SUBSYSTEM  3 


CONFIGURATION 

P 

me 

COST 

(0,  0,  0) 

.6976 

77.1 

(1,  0,  0) 

.7135 

72.0 

(1,1,0) 

.7501 

73.0 

0.  i.  D 

.8019 

74.5 

(1,1,  OR) 

.9002 

78.2 

(1,  2,  OR) 

.9277 

79.2 

(2,  2,  OR) 

.9440 

83.7 

(2,  2,  1R) 

.9530 

89.8 

(2,  2,  2R) 

.9584 

110.4 

These  results  illustrate  the  fact  that  even  if  a  subsystem  can  be  redundant 
it  does  not  follow  that  redundancy  is  the  optimal  decision. 

6.  SUMMARY 

The  DSPC  methodology  represents  a  new  and  innovative  approach  to  system 
acquisition,  and  preliminary  results  indicate  that  this  technique  will  provide 
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very  valuable  information  to  the  decision-maker.  The  DSPC  model  is  compatible 
with  designing  to  system  cost,  or  performance,  or  both.  Once  total  system 
reliability  specifications  are  established,  each  individual  subsystem  has  a 
corresponding  installed  reliability  and  cost  goal,  which  allows  realistic 
and  continuous  evaluation  and  adjustments  as  the  subsystem  is  developed  to 
maturity. 

It  should  be  pointed  out  that  although  the  model  has  been  formulated  in 
terms  of  optimizing  the  jerformance  of  the  total  system,  the  methodology  can 
also  be  profitably  applied  to  individual  subsystems.  For  this  case,  the 
subsystem  is  considered  as  the  total  system  and  its  components  are  considered 
as  the  subsystems.  Then  the  reliability  optimization  procedures  are  applied 
such  that  component  reliability  levels  are  established  such  that  the  desired 
subsystem  reliability  is  achieved. 

As  indicated  above,  the  DSPC  methodology  appears  to  have  great  potential 
in  the  system  acquisition  process.  However,  there  are  two  important  caveats. 
First,  if  the  required  data  are  not  available,  it  will  be  impossible  to  design 
to  required  levels  of  performance  at  minimum  cost.  Second,  assuming  the 
necessary  data  are  available,  if  DSPC  techniques  cannot  be  incorporated  into 
system  acquisition  contracts,  then  it  will  be  impossible  to  achieve  required 
levels  of  performance  at  minimum  cost  except  on  a  chance  basis. 

Preliminary  investigations  by  OAS  indicate  that  a  great  deal  of  data  are 
available  (especially  at  AFLC  Air  Materiel  Areas).  In  some  cases,  rough 
estimates  are  necessary,  but  these  can  be  refined  as  more  emphasis  is  placed 
establishing  and  maintaining  a  DSPC  data  bank.  The  means  of  implementing  DSPC 
techniques  in  contractual  requirements  are  well  beyond  the  scope  of  OAS  efforts 
in  life  cycle  cost  analysis,  but  these  means  must  be  found  if  the  full  potential 
of  the  methodology  is  to  be  realized. 
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5LCT10N  V 


DESIGNING  TO  SYSTEM  PERFORMANCE/COST/EFFECTIVENESS 

1 .  INTRODUCTION 

In  developing  the  DSPC  model  the  objective  function  was  mission  completion 
success  probability,  and  mission  effectiveness  was  not  considered.  It  was 
tacitly  assumed  that  mission  effectiveness  met  mission  requirements.  As 
shown  in  the  previous  section,  the  DSPC  approach  can  be  a  very  valuable 
management  tool;  but  for  a  more  complete  system  evaluation,  mission  effective¬ 
ness  must  also  be  considered.  By  combining  the  results  of  a  DSPC  model  with 
those  from  the  appropriate  mission  effectiveness  model,  information  can  be 
generated  allowing  the  decision-maker  to  more  effectively  evaluate  the  system. 
This  is  particularly  important  when  system  configuration  changes  are  being 
considered,  or  when  there  are  competing  subsystems  during  system  development. 
This  section  presents  a  hypothetical  example  illustrating  performance/cost/ 
effectiveness  interactions. 

2.  PERFORMANCE/COST/EFFECTIVENESS  INTERACTIONS 

Figure  16  provides  an  overall  viewpoint  of  the  methodology  developed  in 
this  study.  The  graphs  show  probability  of  mission  completion  and  mission 
effecti vensss  as  a  function  of  life  cycle  costs  for  three  configurations  -  a 
baseline  configuration  and  two  other  configurations  in  which  subsystems  have 
been  added  in  order  to  increase  mission  effectiveness.  (Effectiveness  is 
defined  as  some  measure  as  to  how  well  a  system  accomplishes  its  mission. 

For  example,  for  weapon  systems  it  is  usually  some  function  of  weapons 
delivery  accuracy  or  targets  killed,  while  for  transport  aircraft  it  would 
generally  depend  on  amount  of  cargo  delivered.) 


79 


Effectiveness 


Configuration,] _ 


Life  Cycle  Cost 
b.  Effectiveness 

Figure  16.  Performance/Cost/Effectiveness  Interactions. 
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For  a  given  system  configuration  a  major  way  of  improving  effectiveness 
is  by  improving  performance,  i.e.,  subsystem  reliabilities.  Assuming  that 
optimum  subsystem  reliability  levels  have  already  been  established  (reliability 
improvement  costs  =  logistic  support  cost  savings),  additional  reliability 
improvement  can  only  be  achieved  with  additional  reliability  improvement  cost 
and  hence,  increased  life  cycle  costs.  In  Figure  16,  the  increase  in  Pmc 
due  to  the  reliability  improvement  is  translated  into  increased  mission 
effectiveness  since  Pmc  is  one  of  the  principal  parameters  in  determining 
mission  effectiveness.  (Reliability  improvement  also  increases  system 
availability  which  is  another  principle  parameter  in  determining  mission 
effectiveness.)  Mission  effectiveness  can  also  be  improved  by  adding  on  other 
subsystems,  for  example,  adding  subsystems  which  improve  weapon  delivery  accuracy. 
As  shown  in  Figure  16-a,  additional  subsystems  increase  life  cycle  costs 
(additional  acquisition  costs  plus  increased  logistic  support  costs)  and 
decrease  system  performance  (overall  system  reliability  is  lowered).  However, 
these  detrimental  effects  may  be  offset  by  increases  in  mission  effectiveness. 

This  is  shown  on  the  Configuration  1  curves  in  Figure  16-b.  On  the  other 
hand,  if  the  potential  benefits  of  a  configuration  change  are  negated  by  a 
decrease  in  system  performance  (P  ),  then  the  modification  results  in 
mission  effectiveness  below  the  baseline  level  as  shown  in  the  Configuration  2 
curves  in  Figure  16-b. 

This  type  of  analysis  makes  it  possible  for  a  decision-maker  to  readily 
evaluate  his  options.  For  example,  in  Figure  16-b  if  available  funds  are 
less  then  C^ ,  then  the  baseline  configuration  is  the  only  option.  If 
additional  funds  are  available,  Configuration  1  is  the  preferred  option 
while  Configuration  2  is  never  in  contention. 

In  the  next  section  examples  of  two  measures  of  effectiveness  for  fighter 


aircraft  are  presented. 
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SECTION  VI 


MEASURES  OF  EFFECTIVENESS  FOR  FIGHTER  AIRCRAFT 

1.  INTRODUCTION 

Previous  sections  have  developed  performance/cost  relationships,  and  the 
last  section  presented  an  overview  of  analyzing  performance/cost/effectiveness 
interactions.  It  is  not  possible  to  develop  generalized  relationships  between 
performance/cost  and  effectiveness  since  effectiveness  depends  on  particular 
systems  and  particular  missions.  Because  of  this  dependence  on  particular 
systems  and  missions,  it  is  sometimes  quite  difficult  to  develop  valid 
measures  of  effectiveness  for  given  systems  and  missions.  Each  system-mission 
combination  must  be  examined,  and  the  system  and  mission  parameters  scrutinized 
to  see  if  meaningful  measures  of  effectiveness  can  be  developed.  In  this 
section  two  measures  of  effectiveness  for  fighter  aircraft  are  developed. 

These  measures  appear  to  have  great  potential  in  fighter  aircraft  evaluations. 

2.  AIR-TO-GROUND  FIGHTERS  (TARGETS  KILLED) 

a.  Characteristic  Effectiveness  Parameters.  The  utility  of  a  tactical 
interdiction  aircraft  is  dependent  upon  the  following: 

(1)  Availability. 

(2)  Probability  of  reaching  target  without  a  critical 
subsystem  aborting  the  mission. 

(3)  Kill  potential  (e.g.,  number  of  targets  destroyed 
per  successful  sortie). 

(4)  Probability  of  survival. 

The  availability  of  an  aircraft  depends  upon  the  frequency  of  repairs  and 
the  average  repair  time  (time  to  restore).  The  probability  of  no  abort  is 
dependent  upon  the  number,  complexity,  and  reliability  of  the  mission 
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critical  subsystems.  The  kill  potential  depends  upon  the  number  and 
type  of  weapon  carried,  acquisition  probability,  delivery  accuracy, 
and  target  type.  Survival  probability  is  dependent  upon  the  strength 
and  type  of  enemy  defenses  and  such  aircraft  characteristics  as  ECM, 
radar  cross  section,  IR  signature,  armor,  and  other  protective  measures. 

The  worth  of  an  aircraft  cannot  be  assessed  by  considering  any  one 
of  the  above  factors  individually.  All  of  the  factors  must  be  considered 
simultaneously  to  account  for  their  interaction.  In  this  section  measures 
of  effectiveness  of  an  aircraft  are  developed  which  quantitatively  account 
for  the  interaction  of  the  characteristic  effectiveness  parameters. 

For  most  mission  types,  an  aircraft  will  be  sent  on  repeated  sorties 
provided  it  survives;  thus,  any  valid  measure  of  effectiveness  must  account 
for  the  cumulative  effect  of  repeated  sorties.  It  is  also  clear,  and 
will  be  shown  quantitatively,  that  survivability  is  of  the  utmost  importance 
since  it  determines  the  average  number  of  sorties  an  aircraft  can  complete. 

If  a  particular  scenario  is  specified,  then  for  a  given  aircraft 
the  characteristic  effectiveness  parameters  serve  to  characterize  that 
aircraft  and  scenario.  Since  actual  scenarios  change  sortie  by  sortie, 
the  determination  of  aircraft  performance  over  repeated  sorties  requires 
that  the  characteristic  parameters  be  specified  for  each  sortie.  Such 
a  detailed  sped  i'i cation  would  introduce  a  high  level  of  arbitrariness 
leading  to  an  unsuitable  measure  of  a  system's  worth.  However,  to  obtain 
a  measure  (not  a  predictor)  of  the  effectiveness  of  an  aircraft  in  a 
given  scenario  it  seems  reasonable  to  keep  the  scenario  fixed  (fixed 
characteristic  parameters)  and  to  determine  the  cumulative  effectiveness 
if  the  aircraft  flies  repeated  sorties  (when  it  survives)  in  that  fixed 
scenario.  This  is  the  basic  idea  underlying  the  measures  of  effectiveness 
uevaloped  in  this  section. 
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The  first  measure  developed  is  the  expected  number  of  targets  destroyed 
if  the  aircraft  flies  up  to  S  sorties  in  a  fixed  scenario,  i.e.,  aircraft 
flies  repeated  sorties  (up  to  a  maximum  of  S)  if  it  survives.  The  next 
measure  is  the  lifetime  targets  destroyed,  i.e.,  aircraft  flies  repeated 
sorties  as  long  as  it  survives.  However,  the  expected  number  of  targets 
destroyed  during  the  lifetime  of  an  aircraft  may  not  be  the  prime  measure 
of  effectiveness  since  there  are  situations  in  which  it  is  more  important 
to  Know  the  effectiveness  of  an  aircraft  over,  for  instance,  a  10  or 
20  day  period.  The  final  measure  of  effectiveness  developed  is  the  expected 
number  of  targets  destroyed  as  a  function  of  time,  which  yields  targets 
destroyed  over  any  prescribed  time  period. 

Although  the  discussion  is  in  terms  of  tactical  interdiction  aircraft, 
the  kill  potential  can  be  redefined  (for  example,  in  terms  of  cargo  tonnage 
delivered  or  enemy  aircraft  destroyed)  to  account  for  airlift,  counterair, 
or  other  type  aircraft. 

b.  Lifetime  Destruction.  The  definitions  listed  below  will  facilitate 
the  mathematical  developments  contained  in  this  section. 

Psj  =  Probability  aircraft  survives  to  release  its  weapons  on  target. 

P$2  =  Probability  aircraft  survives  return  trip  after  weapons  are 
released. 

Pc  =  Probability  aircraft  reaches  target  and  releases  weapons 
without  an  abort  causing  failure  given  that  it  survives. 

F$a  =  Probability  aircraft  aborts  before  releasing  weapons  and 
survives  the  return  trip. 

c  =  "Kill  Potential"  =  expected  number  of  targets  destroyed 
after  aircraft  reaches  the  target  area. 

P$  =  Single  sortie  survival  probability. 


* 
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S  =  Number  of  'sorties  aircraft  flies  (if  it  survives). 

T(S)  =  Expected  number  of  targets  destroyed  after  S  sorties. 

=  Expected  number  of  targets  destroyed  during  the  "lifetime" 
of  the  aircraft,  i.e.,  S  «*. 

For  a  single  sortie,  the  expected  number  of  targets  destroyed  by  an 
aircraft  is 


TO) 


cPcPsl 


(VI-1) 


The  main  problem  in  this  section  is  to  determine  the  expected  number  of 
targets  destroyed  if  the  aircraft  flies  a  maximum  of  S  sorties  (if  it 
survives).  The  time  required  to  complete  S  sorties  is  treated  in  the  next 
section.  The  probability  Pi  that  the  aircraft  starts  its  i -th  sortie 
(i  <_  S)  is  equivalent  to  the  probability  it  survives  the  first  i-1  sorties. 
Therefore, 


H-l 


Pi  =  i  PslPcPs2  +  Psa | 


*  P)'\  (1  ■  1.2....S),  (VI-2) 


where  P$  denotes  the  single  sortie  survival  probability.  The  expected  damage 
from  the  i-ch  sortie  is 

PicPCPs!  =  Ps''”PCP5l  •  (VI-3> 


Therefore,  it  follows  that  the  expected  number  of  targets  destroyed  after  S 
sorties  is 


T(S)  « 


E  Wsi 


i  =  1 


■  =PcPsl 


E  pr’ 


oPcPsl 


i  =  1 


1  '  Ps 


( VI  -  4) 
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Letting  S  -  -  in  equation  (VI -4)  it  follows  that  the  expected  number  of 
targets  destroyed  during  the  lifetime  of  the  aircraft  is 

pPA, 

tk  ■  t4t7  •  («-5) 

Of  course,  if  for  any  reason  there  is^an  upper  limit  to  the  number  of  sorties 
the  aircraft  would  fly,  then  this  number  should  be  used  in  equation  (VI-4) 
to  determine  the  expected  damage  during  the  useful  lifetime  of  the  aircraft. 

The  expression  (VI-5)  for  lifetime  destruction  was  derived  under  the 
assumptions  that  the  aircraft  flies  repeated  sorties  as  long  as  it  survives 
and  that  the  scenario  remains  the  same  for  each  sortie.  It  is  important  to 
point  out  that  this  measure  of  effectiveness  has  another  interpretation. 
Suppose  N(N  =  1,  2,  3,  ...)  aircraft  each  fly  one  sortie  where  the  parameters 
p,  P  ,  P  ,,  and  P  are  the  same  for  each  aircraft.  The  expected  number  of 

t  a  I  j 

targets  destroyed  by  the  N  aircraft  is 

NpPcPs1  .  (VI -6) 


The  expected  number  of  aircraft  lost  is 

NO  -  Ps)  .  (VI-7) 

The  ratio  of  the  quantities  (VI-6)  and  (VI-7)  yields  a  measure  of  targets 
destroyed  per  aircraft  lost  (exchange  ratio)  equal  to 

ppcPsl 

T-fr:  •  ( vi  -8) 

which  is  independent  of  the  number  of  aircraft.  This  exchange  ratio  is 
identical  to  expression  (VI-5)  for  lifetime  targets  killed. 
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The  expected  number  of  sorties  completed  during  the  lifetime  of  an 


aircraft  is 


<s>=  ^  jpsjn  -  ps)  =  •  (VN9) 

j  =  i 


T.iis  Measure  is  further  discussed  in  the  examples  in  Section  VI-2d. 

Although  the  measures  (VI-4)  and  (VI-5)  are  useful  indicators  of  the 
effectiveness  of  an  aircraft,  they  do  not  reflect  the  time  rate  of  damage. 

This  is  the  subject  of  the  next  section. 

c.  Targets  Killed  as  Function  of  Time.  Equation  (VI-4)  gives  the  expected 
number  of  targets  destroyed  after  S  sorties.  However,  in  evaluating  the 
effectiveness  of  an  aircraft,  it  is  also  essential  to  determine  the  expected 
time  required  for  the  S  sorties.  This  time  depends,  vOf  course,  upon  the 
mission  time  T|t)  and  also  upon  the  time  required  to  make  repairs. 

If  the  aircraft  completes  S  sorties  then  the  expected  number  of  repairs  is 


where  ts  is  the  MTBF  of  the  total  aircraft  system. 
Therefore,  the  expected  total  repair  time  is 


(VI-10) 


ST  t 
m  r 


(VI-11) 


wnere  t  is  the  mean  time  to  restore.  If  At  denotes  the  average  time  for 
normal  service  actions,  e.g.,  refuel  and  reload,  then  the  expected  time  to 
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complete  S  sorties  is 


t(s)  =  S 


T  +  JLI  +  At 
m  t 


(VI-12) 


I 


It'  the  service  actions  can  be  performed  while  repairs  are  being  made, 
tuen  At  in  equation  (VI-12)  should  be  replaced  by 


min 


it  -  ^  ,  0 


(VI-13) 


The  Reliability  Engineering  Handbook,  (Reference  4),  defines  avail' 
ability  as 


A  = 


1 


1 


t 

T 


(VI-14) 


From  this  it  follows  that 


I'  1 


(VI-15) 


Therefore,  equation  (VI-12)  becomes 


t(s)  =  S 


(VI-16) 


Equations  (VI-4)  and  (VI-16)  provide  the  expected  number  of  targets  destroyed 
as  a  function  of  time. 

In  the  following  section,  examples  will  be  given  to  show  how  the  individual 
characteristic  parameters  associated  with  an  aircraft  interact  in  determining 
the  effectiveness  of  an  aircraft. 
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d .  Examples . 

a.  Lifetime  Sorties. 

Figure  17  shows  the  expected  number  of  sorties  completed  during 

the  lifetime  of  an  aircraft  as  a  function  of  survival  probability.  Since  the 

lifetime  targets  killed  T„  is  a  constant  factor  multiplied  by  life-time 

sorties,  the  curve  for  L,  has  the  same  shape  as  the  curve  in  Figure  17.  Of 

1\ 

course,  the  curve  cannot  be  extended  indefinitely  since  there  is  an  upper 
limit  based  upon  the  service  life  of  the  aircraft  or  other  such  factors. 

Several  conclusions  are  apparent: 

(1)  Conditions  resulting  in  survival  probabilities  below 
.95  are  probably  unacceptable  in  most  cases  since  lifetime  sorties  is 
less  than  19. 

(2)  Small  improvements  in  survival  probability  in  the  region 
P,  .98  result  in  a  small  increase  in  lifetime  sorties.  However,  in 

the  region  of  high  (e.g.,  P$  -  .98)  any  small  increase  in  results 
in  a  dramatic  increase  in  lifetime  sorties.  For  example,  the  small  increase 
in  P$  from  .99  to  .995  more  than  doubles  the  number  of  lifetime  sorties 
( from  99  to  199) . 

(3)  Survival  probability  can  be,  by  far,  the  most  dominant 
factor  in  determining  the  effectiveness  of  aircraft. 

To  appreciate  the  magnitude  of  the  numbers  involved,  it  is 
instructive  to  consider  a  historical  but  recent  engagement  in  a  severe 
environment  where  U.S.  aircraft  flew  1000  sorties  against  heavily  defended 
targets.  During  this  period,  26  U.S.  aircraft  were  lost.  The  survival 
probability  in  this  case  was  P$  =  0.974  which  is  on  the  low  part  of  the 
curve  in  Figure  17.  Under  such  conditions  the  average  number  of  sorties 
ner  aircraft  is  only  37.5. 
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Figure  17.  Lifetime  Sorties  as  a  Function  of  Survival  Probability. 
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One  obvious  me <ins  to  increase  survival  probability  is  to  reduce 


the  enemy's  defenses  (gain  air  superiority).  Survival  probability  can 
also  be  improved  by  designing  the  aircraft  to  reduce  the  probability  of 
hit  (e.g.,  ECM  or  reducing  radar  cross  section  and  IR  signature)  and  to 
reduce  the  probability  of  kill  given  the  aircraft  is  hit  (e.g.,  armor, 
foam  in  fuel  tanks).  As  shown  above,  any  improvements  in  aircraft  sur¬ 
vivability  significantly  enhances  mission  effectiveness, 
b.  Comparing  Aircraft. 

Table  XV  shows  the  effectiveness  parameters  associated  with  5 
hypothetical  aircraft  labeled  A,  B,  C,  D,  and  E.  Although  each  parameter 
is  important  in  the  evaluation  of  an  aircraft,  it  appears  impossible  to 
rank  the  5  aircraft  by  studying  the  table.  The  table  does  show  that 
aircraft  A  has  the  best  kill  potential,  B  has  the  highest  probability 
of  reaching  the  target  without  an  abort,  C  has  the  highest  survival  probability, 
and  D  has  the  highest  availability. 


Table  XV 


EFFECTIVENESS 

PARAMETERS 

FOR  FIVE 

AIRCRAFT 

EFFECTIVENESS  PARAMETERS 

A  I 

R  C  R 

A  F  T 

T 

Y  P  E 

A 

B 

C 

D 

E 

Kill  Potential  (e.g.,  targets 
killed  per  successful  sortie) 

2.5 

1.8 

.80 

2.0 

2.1 

P  (Probability  of  reaching 
target  without  an  abort) 

.90 

.93 

.90 

.82 

.85 

P$  (Survival  Probability) 

.970 

.990 

.999 

.980 

.995 

A  (Availability) 

.85 

.87 

.83 

.90 

.83 

Mission  Time  (hr) 

2 

2 

2 

2 

2 

Service  Time  (hr) 

.5 

.5 

.5 

.5 

.5 

9  2 


Using  the  parameters  listed  in  Table  XV  together  with  equations 
(VI -4)  and  (VI -16),  the  expected  number  of  targets  destroyed  as  a  function 
ot  time  can  be  calculated  for  each  aircraft.  The  results  in  Figure  IB  are 
based  upon  continuous  operation,  i.e.,  aircraft  is  launched  as  soon  as  it 
is  ready.  Although  aircraft  E  does  not  dominate  the  others  in  any  of  the 
effectiveness  parameters,  when  all  parameters  are  integrated  aircraft  E  is 
superior  to  the  others  (E  and  A  are  about  equal  in  the  beginning)  at  least 
for  time  less  than  50  hours  (about  25  missions).  The  lifetime  targets 
destroyed  by  each  aircraft  indicates  where  the  curves  finally  level  off. 

The  lifetime  targets  destroyed  (LTD)  by  each  aircraft  are: 

A:  LTD  =  73 
B:  LTD  =  166 
C:  LTD  =  719 
D:  LTD  =  80 
E:  LTD  =  355 

This  indicates  that  C  might  be  better  than  E  since  its  curve  will  eventually 
rise  above  the  targets  destroyed  curve  of  aircraft  E.  Figure  19  shows  targets 
destroyed  by  C,  E,  and  B  as  a  function  of  time  when  time  is  carried  out  to 
2500  hours  (about  859  missions).  Although  aircraft  C  and  E  are  the  only  two 
competitors,  aircraft  B  is  shown  merely  to  demonstrate  that  its  low  surviv¬ 
ability  causes  its  curve  *o  level  off  early  at  a  LTD  of  166.  Figure  19  shows 
that  E  is  substantially  better  than  C  for  times  less  than  1830  hours  (629 
missions).  For  times  greater  than  this  the  higher  survivability  of  C  more 
than  compensates  for  its  lower  kill  potential  and  C  is  better  than  E.  The 
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Targets 

Destroyed 


Figure  19.  Destruction  as  a  Function  of  Time  for  Aircraft  B,  C,  and  E. 
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The  analysis  shows  that  E  is  better  than  aircraft  A,  B,  and  D.  However,  the 
selection  between  C  (with  a  higher  LTD)  and  E  is  dependent  upon  the  preference 
of  the  decision-maker,  i.e.,  whether  short  term  or  long  term  performance  is 
of  prime  interest. 

3.  AIR-TO-AIR  FIGHTERS  (IMPORTANCE  OF  FIRST  SHOT) 

It  is  particularly  interesting  to  apply  some  of  the  ideas  of  the  previous 
section  to  air-to-air  engagements  between  fighter  aircraft.  It  is  intuitively 
clear  that  the  probability  of  maneuvering  into  position  to  fire  the  first  shot 
is  an  important  factor  in  determining  the  effectiveness  of  a  fighter  aircraft. 
This  is  due  to  the  fact  that  the  first  shot  probability  has  a  strong  influence 
on  both  the  kill  potential  and  the  survival  probability  of  the  fighter.  The 
tools  developed  in  the  previous  sections  provide  a  means  to  show  quantitatively 
the  influence  of  first  shot  probability  on  the  exchange  ratio  (i.e.,  Red 
fighters  destroyed  per  Blue  fighter  destroyed).  This  exchange  ratio  can 
also  be  interpreted  as  the  expected  number  of  enemy  fighters  destroyed  during 
the  lifetime  of  a  Blue  fighter.  The  first  air-to-air  scenario  is  described 
in  the  next  paragraph. 

In  an  air-to-air  engagement  between  a  Blue  and  a  Red  fighter,  the 
probability  that  the  Blue  fighter  fires  the  first  shot  is  denoted  by  P-j . 

This  first  shot  probability  is  a  function  of  acquisition  and  tracking 


capabilities,  speed,  maneuverability,  and  pilot  skills.  The  fighter  firing 
the  first  shot  releases  its  air-to-air  weapons  destroying  the  other  fighter 


with  a  certain  probability  (P^  for  Blue  weapons,  P^ 


for  Red  weapons). 


If 


the  attacked  fighter  is  destroyed  the  engagement  is  finished;  however,  for 


this  first  scenario  it  is  assumed  that  if  the  attacked  fighter  is  not  des¬ 


troyed  it  maneuvers  into  position  to  launch  its  weapons  against  the  other 
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aircraft  (this  assumption  will  be  modified  later).  The  engagement  is  then 
finished  with  each  fighter  getting  at  most  one  pass.  Although  multiple 
passes  could  easily  be  considered,  it  requires  additional  assumptions  and 
contributes  little  to  the  understanding  of  the  problem  (especially  if  both 
fighters  are  assumed  to  have  highly  effective  air-to-air  weapons). 

The  first  quantity  to  be  derived  is  the  probability  that  the  Blue 
fighter  destroys  the  Red  fighter  in  a  given  engagement.  This  is  the  fighter 
kill  potential;  it  is  equal  to  the  probability  that  the  Blue  fighter  fires 
the  first  shot  and  destroys  the  Red  fighter  plus  the  probability  that 
the  Red  fighter  fires  the  first  shot  and  misses  the  Blue  fighter  and  the 
Blue  fighter  then  destroys  the  Red  fighter.  Thus, 

P  ■  F,Pkb  +  (1  -  P,)0  -  Pkr)Pkb 

3  pkbj'  •  V  -pi>}  ■ 

where  denotes  the  first  shot  probability  of  the  Blue  fighter,  P^b 
is  the  kill  probability  of  the  weapons  of  the  Blue  fighter,  and  Pkr  is 
the  kill  probability  of  Red  weapons. 

The  next  expression  to  be  derived  is  the  single  engagement  survival 
probability  Pg  of  the  Blue  fighter.  The  Red  fighter  will  be  prevented  from 
launching  its  weapons  only  if  the  Blue  fighter  gets  the  first  shot  and 
destroys  the  Red  fighter.  Therefore,  the  probability  that  Red  attacks 
tne  Blue  fighter  is 


1  - 


PlPkb 


(VI-18) 
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Tlie  survival  probability  is  then 


\ 


ps  ■  '  -  pkr<’  -  plPkb>  • 


From  equations  (VI-17)  and  (VI-19)  It  follows  that  the  exchange  ratio  is 
given  by 


ER  = 


rrr. 


s 


rrr T7P 


1  kb 


(VI -20) 


Figure  20  is  presented  to  illustrate  the  strong  influence  of  first  shot 
probability  on  both  the  probability  of  survival  of  the  Blue  fighter  and  the 
probability  of  survival  of  the  Red  fighter.  In  this  example  the  effectiveness 
of  Red  and  Blue  weapons  is  assumed  to  be  equal,  i*e-'  Pkb  “  Pkr  *  0,9‘ 

Although  weapons  are  equally  effective,  the  first  shot  capability  of  Blue  can 
cause  the  survival  probability  to  vary  from  0.10  to  0.91  and  the  kill  prob¬ 
ability  against  the  Red  fighter  to  vary  from  0.09  to  0.90. 

Figure  21  incorporates  both  the  kill  potential  and  survival  probability 
to  show  the  dependence  of  the  exchange  ratio  upon  the  first  shot  capability. 
Two  cases  are  presented  corresponding  to  Red  weapon  effectiveness  of 
P^r  -  0.6  and  =  0.9.  For  each  case  the  exchange  ratio  is  plotted  for 
pkb  =  0,6  and  Pkb  =  Several  conclusions  are  apparent: 

(1)  Effective  Blue  weapons  and  a  high  first  shot  capability  are 
both  necessary  for  achievement  of  a  high  exchange  ratio  for  Blue. 

(2)  Even  when  P^  *  0.9  and  P^  =  0.6  a  first  shot  probability 
below  P,  *  0.22  results  in  an  exchange  ratio  below  1.0,  i.e.,  the  advantage 
of  a  superior  weapon  can  be  nullified  by  a  poor  first  shot  capability. 
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Figure  20.  The  Effect  of  First  Shot  Probability  on 

Survival  and  Kill  Probability  (Pb  =  Pkr  =  0.9), 
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Figure  21.  Exchange  Ratio  as  a  Function  of  First  Shot 
Probability. 
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(3)  The  disadvantage  of  a  poor  weapon  (e.g.,  =  0.6  and 

Pkr  =  0.9)  can  sometimes  be  more  than  compensated  for  by  a  high  first  shot 
capability. 

In  the  previous  scenario  it  was  assumed  that  whenever  the  fighter  firing 
the  first  shot  missed,  the  other  fighter  then  maneuvered  into  position  to 
fire  its  weapons.  However,  a  fighter  may  fire  the  first  shot  and  miss  but 
still  have  the  capability  to  outmaneuver  the  other  fighter  thereby  avoiding 
being  fired  upon.  To  account  for  this,  the  following  probabilities  are 
introduced: 

Pmb  =  Pr°bab^ty  B^ue  ffqhter  avoids  being  fired  upon  whenever  it 
fires  first  shot  and  misses. 

Pmr  *  Probability  Red  fighter  avoids  being  fired  upon  whenever  it 
fires  first  shot  and  misses. 

The  probability  that  the  Red  fighter  is  destroyed  becomes 

»*pkb{pl  +  <’ -pkr>(' -V}  • 

The  probability  that  the  Blue  fighter  is  destroyed  is 


'  -  Ps  3  pi(l  -  PhJO  -  P„JP„,  Ml  -  pi)P 


kb 


nt> 7  kr 


1 7  kr 


=  P 


kr  I1  '  pl(pkb  '  rmb  "  Wkb 


■r  P  j.  “  Pj.Pi, 


>}  • 


(VI-22) 


From  equations  (VI-21)  and  (VI-22)  it  follows  that  the  exchange  ratio  is 


ER 


.  Pkb{Plt(1-P!>(1-pkr><1-pnr>) 
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r  kb  mb  mb  kb' 


F 


(VI -23) 


For  =  pmr  =  equation  (VI -23)  reduces  to  equation  (VI -20). 
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The  most  favorable  case  for  Blue  is  when  =  1  and  Pmr  =  0;  the  most 
unfavorable  case  is  Pn|b  =  0  and  Pmr  =  1.  Using  these  extreme  cases,  the 
bounds  for  the  exchange  ratio  are  shown  in  Figure  22.  The  solid  curves  are 
identical  to  those  in  Figure  21 ,  i  .e.,  P^  =  Pmr  =  0.  As  seen  by  comparing 
the  solid  and  lower  curves  in  the  Figures,  If  Red  can  outmaneuver  Blue  after 
getting  first  shot  but  Blue  does  not  have  this  capability  (P^  =  0,  Pmr  =  1) 
this  has  little  effect  on  the  exchange  ratio  since  it  has  no  effect  on 
Blue's  survival  probability.  However,  if  P^  =  1  and  Pmr  =  0,  Blue  can 
improve  its  survival  probability,  and  hence  the  exchange  ratio  is  improved 
significantly  if  the  first  shot  probability  is  high;  furthermore,  the  lower 
the  value  of  P^  the  greater  the  importance  of  the  capability  of  Blue  being 
able  to  outmaneuver  Red  after  firing  the  first  shot. 

4.  CONCLUSIONS 

(1)  An  evaluation  of  the  effectiveness  of  an  aircraft  must 
account  for  the  interaction  of  availability,  abort  probability,  kill  potential, 
and  survivability.  Individually,  these  characteristic  parameters  do  not 
determine  the  worth  of  an  aircraft. 

(2)  Any  valid  measure  of  effectiveness  must  also  account  for 
the  cumulative  effect  of  repeated  sorties. 

(3)  The  measures  of  effectiveness  developed  here  provide  a 
simple  means  of  integrating  the  characteristic  effectiveness  parameters  to 
determine  the  cumulative  damage  accrued  by  repeated  sorties. 

(4)  Survival  probability  can  be  the  most  dominant  factor  in 

determining  the  lifetime  effectiveness  of  an  aircraft.  For  example,  a  5% 

increase  in  kill  potential  results  in  a  5%  increase  in  lifetime  damage; 

however,  a  5%  increase  in  survival  probability,  say  from  P„  =  .95  to 

s 

Ps  =  .9975,  results  in  a  2100%  increase  in  targets  destroyed  during  the 
lifetime  of  the  aircraft. 
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Figure  22.  Importance  of  Maneuverability  after 
First  Shot  (Pkr  =  0.6), 
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(5)  Since  survivability  is  of  such  great  importance  it.  warrants 
special  emphasis  during  design  and  testing.  Survival  probability  is  an 
extremely  important  factor  in  comparing  two  aircraft;  for  instance,  one 
aircraft  may  have  a  poorer  weapon  delivery  accuracy  and  yet  be  far  superior 
because  of  higher  survivability. 

(6)  In  the  case  of  air-to-air  fighter  aircraft  the  exchange 
ratio  (Red  aircraft  destroyed  per  Blue  aircraft  destroyed)  is  an  important 
measure  of  worth. 

(7)  The  exchange  ratio  for  air-to-air  fighter  aircraft  can  be 
expressed  as  a  function  of  three  fundamental  parameters:  weapon  effectiveness, 
first  shot  probability,  and  the  capability  to  maneuver  away  (avoid  being 
fired  upon)  after  firing  first  shot  and  missing. 

(8)  The  most  important  parameter  affecting  the  exchange  ratio  is 
the  first  shot  probability.  The  advantage  of  a  superior  weapon  can  be 
nullified  by  a  poor  first  shot  capability;  and,  conversely,  the  disadvantage 
of  an  inferior  weapon  can  sometimes  be  compensated  for  by  a  good  first  shot 
capability. 

5.  MULTI-ROLE,  MULTI-MISSION  CAPABILITY 

The  measures  of  effectiveness  discussed  above  provide  very  useful  tools 
for  evaluating  different  candidate  systems  to  be  used  in  performing  the  same 
missions.  However,  when  the  systems  under  consideration  have  a  multi-role 
or  mission  capability,  a  more  elaborate  means  of  evaluating  the  effectiveness 
of  each  system  must  be  used.  For  example,  when  evaluating  two  systems  which 
can  perform  both  the  air-to-air  and  air-to-ground  missions,  some  means  must  be 
devised  to  allow  the  tradeoff  to  be  made  by  evaluating  the  system's  effectiveness 
under  both  roles.  Unless  one  system  clearly  dominates  the  other  in  both  roles, 

/ 

some  measure  of  effectiveness  must  be  used  which  allows  their  effectiveness  to 
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be  evaluated  in  conjunction  with  their  ability  to  complement  the  existing  base 
force  and  its  capabilities  in  performing  the  air-to-air  and  air-to-ground 
missions.  That  is,  the  current  base  force  may  be  weak  in  the  air-to-air 
capability  and  thus  the  addition  of  the  superior  air-to-air  system  would  be 
preferred.  On  the  other  hand,  by  adding  a  better  air-to-ground  system  to  a 
force  which  lacks  this  capability,  it  may  be  possible  to  allow  better  utiliza¬ 
tion  of  an  existing  air-to-air  capability  which  was  previously  needed  in  the 
air-to-ground  role. 

6.  SUMMARY 

Two  measures  of  effectiveness  for  fighter  aircraft  have  been  presented. 

In  the  case  of  air-to-ground  fighters,  it  was  shown  that  an  evaluation  of  the 
effectiveness  must  account  for  the  interaction  of  availability,  abort  probability, 
kill  potential,  and  survivability;  and  survivability  is  often  the  most  dominant 
factor.  For  air-to-air  fighters  the  exchange  ratio  (Red  aircraft  destroyed 
per  Blue  aircraft  destroyed)  is  an  important  measure  of  worth,  and  it  can  be 
expressed  as  a  function  of  weapon  effectiveness,  maneuver  capability,  and 
first  shot  probability  with  first  shot  probability  being  the  most  important 
parameter. 
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SECTION  VII 


SUMMARY 

In  summary,  Figure  23  presents  a  progression  of  the  new  and  innovative 
approach  for  obtaining  higher  operational  reliability  levels  using  the  models  and 
methodology  developed  in  this  study. 

Each  block  represents  a  necessary  step  in  the  process,  and  continuous  feed¬ 
back  and  iteration  is  required  to  realize  the  full  potential  of  the  approach. 

The  feedback  loop  extends  from  any  one  block  to  any  preceding  block.  By  estab¬ 
lishing  this  sequence  with  the  appropriate  feedback  and  iteration,  requirements 
and  achievable  operational  levels  can  be  kept  compatible. 

Starting  with  the  requirements  and  proceeding  to  reliability  testing  must 
involve  a  great  deal  more  than  MIL-STD-781B  demonstrations.  If  Initial  Operational 
Test  and  Evaluation  (IOT&E)  results  are  not  available,  then  the  equipment  should 
be  stressed  in  the  laboratory  in  such  a  way  as  to  uncover  as  many  reliability 
deficiencies  as  possible.  If  test  results  Indicate  that  the  equipment  in  its 
original  configuration  will  not  meet  operational  requirements,  then  the  requirements 
can  either  be  adjusted  or  reliability  improvement  programs  can  be  undertaken.  The 
MCSP  model  is  then  used  to  evaluate  the  original  configuration  by  identifying  the 
critical  components  and  determining  the  effect  of  critical  component  improvement 
on  overall  system  reliability.  The  next  step  is  to  determine  realistic  funding 
levels  for  the  reliability  Improvement  program.  This  is  accomplished  with  a 
reliability  management  program  in  which  reliability  options  and  logistic  support 
costs  are  considered. 

With  this  data  the  DSPC  model  can  be  implemented.  This  methodology  systemati¬ 
cally  identifies  those  subsystem  options  which  provide  the  highest  system  perfor¬ 
mance  at  any  prescribed  level  of  cost.  Along  with  the  DSPC  methodology  appropriate 
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Figure  23.  Implementing  the  Methodology  to  Achieve 
Higher  Operational  Reliability  Levels. 
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measures  of  effectiveness  must  be  tailored  to  the  particular  mission  of  interest 
and  related  to  the  system  performance  parameters. 

Two  measures  of  effectiveness  for  fighter  aircraft  have  been  presented.  In 
the  case  of  air-to-ground  fighters,  it  was  shown  that  an  evaluation  of  the  effec¬ 
tiveness  must  account  for  the  interaction  of  availability,  abort  probability, 
kill  potential,  and  survivability;  and  survivability  Is  often  the  most  dominant 
factor.  For  air-to-air  fighters,  the  exchange  ratio  (Red  aircraft  destroyed  per 
Blue  aircraft  destroyed)  is  an  important  measure  of  worth,  and  it  can  be  expressed 
as  a  function  of  weapon  effectiveness,  maneuver  capability,  and  first  shot 
probability  with  first  shot  probability  being  the  most  important  parameter. 
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APPENDIX  A 


MISSION  COMPLETION  SUCCESS  PROBABILITY  (MCSP)  COMPUTER  PROGRAM 
1.  DEFINITIONS  OF  INPUTS 

NSSYS  =  Total  number  of  subsystems  (NSSYS  <_  50) .  Each  subsystem  is 

identified  by  number  and  name. 

NPHASES  =  Number  of  mission  phases  (NPHASES  <  20). 

T0(i,  j)  =  Operating  Time  of  i-th  subsystem  during  j-th  mission  phase 
(i  =  1,  2,  ....  NSSYS;  j  =  1,  2,  ....  NPHASES). 

PA(i,  j)  =  Conditional  probability  of  mission  abort  given  that  the  i-th 
subsystem  has  a  failure  during  the  j-th  mission  phase 
(i  =1,2 . NSSYS;  j  =  1,  2,  ....  NPHASES). 

INDl(i)  =  Indicator  equal  to  0  or  1 .  If  INDl(i)  =  0  then  the  baseline 
for  the  1-th  subsystem  is  nonredundant.  If  IND1 (1 )  =  1  then 
the  baseline  for  the  1-th  subsystem  is  redundant. 

F(i,  1)  =  When  IND1 (1 )  *  0,  this  parameter  is  required  to  denote  the 

mean  operating  time  between  failure  of  subsystem  i. 

NR(i,  1)  h  When  IND1 (1 )  =  1,  this  parameter  is  required  to  denote  the 
number  of  redundant  units  for  subsystem  i  (NR(i,  1)  ^  5). 

R(i,  1)  =  When  IN01 (1 )  =  1,  this  parameter  is  defined  as  0  or  1 .  If 

R(i,  1)  =  0  then  the  redundancy  for  subsystem  i  is  operative. 
If  R(i,  1)  =  1  then  the  redundancy  is  standby. 

FR(i,  1,  k)  =  When  IND1  ( 1 )  =  1,  this  parameter  denotes  the  mean  operating 
time  between  failure  of  the  k-th  redundant  subsystem  of 
subsystem  1  (k  =  1,  2,  ...»  NR(i,  1)). 
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INO  2  -  Indicator  equal  to  0  or  1 .  If  IND  2  =  1  then  the  subsystems 

will  be  ranked  according  to  their  probability  of  causing  an 
abort. 

IND  3  =  Indicator  equal  to  0  or  1 .  If  IND  3  =  1  then  a  sensitivity 

analysis  will  be  performed  for  improvement  of  a  selected 
non redundant  subsystem. 

IS  =  When  IND  3  =  1 ,  this  identifies  the  number  of  the  nonredundant 

subsystem  whose  MTBF  is  to  be  incremented. 

DELTA  =  When  IND  3  =  1 ,  this  denotes  the  size  of  the  MTBF  increment 
for  subsystem  S. 

XLIMIT  =  When  IND  3  *  1 ,  this  denotes  the  upper  limit  of  the  MTBF 
increment  for  subsystem  S. 

2.  OUTPUTS  OF  THE  MCSP  MODEL 

The  model  calculates  the  probability  of  mission  completion  without  an 
abort  causing  falljre  of  a  subsystem.  If  IND  2*1  then  the  subsystems  will 
be  ranked  according  to  the  probability  of  an  abort  causing  failure  of  each 
subsystem.  If  IND  3=1  then  the  mean  operating  time  between  failure  of  one 
nonredundant  subsystem  will  be  incremented  and  the  corresponding  values  of 
MCSP  will  be  calculated.  To  perform  a  sensitivity  analysis  on  a  redundant 
subsystem  the  model  could  be  exercised  repeatedly  making  the  appropriate 
changes  in  the  F(i,  1,  k)  for  each  case.  It  should  be  mentioned  that  in 
the  case  of  standby  redundancy  of  a  subsystem  the  corresponding  MTBFs 
(FR(i,  1,  k),  k  *  1,  2,  ....  NR(i ,  1))  must  be  input  as  either  all  equal  or 
all  unequal. 

3.  MCSP  CARD  INPUTS 

The  card  inputs  to  the  MCSP  model  are  identical  to  those  for  the  DSPC 

model  (described  in  Appendix  B)  with  a  few  exceptions. 
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There  are  fewer  arrays  to  be  input,  with  several  having  a  constant  rather 
than  variable  dimension.  The  same  DSPC  rules  for  array  input  hold  for  these 
arrays  with  the  appropriate  dimension  set  to  1. 

There  is  an  extra  "array"  to  be  input  to  this  program.  The  mnemonic  is 
INO  and  the  fields  are: 


FIELD 

1 

2 

3 

4 

5 


VALUE  OF 
IND  2 
IND  3 
IS 

DELTA 

XLIMIT 
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o  r>  r»  r>  nm 


4  -  MCSP  PROGRAM  LISTING 


PROGRAM  MCSP ( INPUT  *OUTPnT «TaPE60*INPUT ) 
PRr»GP«M  MrSP(lNPI»T.OUTPUT.TAPE60*TNPUT> 

INPUT  THF  APPAYS 
CALL  INPUT 

NOW  EXECUTE  THE  MAIN  800V 

CALL  CONTROL 

THATS  ALL  FOLKS.... 

FNO 


tnor>o<“  u  n  o  n  m  >*000  r»r»r> 


SUBROUTINE  ImPIJT 


SlieROUTINr  JMPIJT  r  1 

COMMON  NSSYS.NPHASFS* IN02.  TN03«TS»DELTA»X|_IMIT»TO(50»20),PA(50»2P)  R  2 

I.NAmFMSOI.INOI  (SO) .NR(50«1>,W(50.1>.F<SO,1>,FR(50,1.6>,TP<50).P(5  R  3 

20.1 > .PMC(SO)  R  4 

DIMENSION  CAPOS(A)  R  5 

JAROPTsO  R  6 

fl  7 

READ  THE  FIRST  CARO  R  B 

R  9 

PFA3  19*  CAPOS  R  1 0 

NSSYS=CARnS (?)  R  U 

MPHASPS=CAHDSn)  R  12 

PPINT  20.  NSSYS.NPHASFS  R  13 

PFAO  PI.  (NAMFSm  .1=1. NSSYS)  R  14 

PPINT  22.  (NAMES(I) .1*1. NSSYS)  R  15 

R  16 

PFOIN  RFAOINC,  THE  ARRAYS  R  17 

R  1 R 

PFAO  23.  TAN  R  19 

IF  (EOF (60) )  1R.2  R  ?0 

PRINT  24.  IAN  R  21 

IF  ( I AN.E0.4HT0  »  GO  TO  3  R  22 

IF  (IAN.F0.4HPA  )  GO  TO  5  R  23 

IF  (IAN.F0.4HIN01)  GO  TO  7  R  24 

IF  (IAN.F0.4HIN0  )  GO  TO  9  R  25 

IF  (IAN.E0.4HNR  )  GO  TO  IP  R  26 

IF  ( IAN.F0.4HR  )  GO  TO  12  R  27 

IF  (IAN.E0.4HFR  )  GO  TO  14  R  28 

IF  (IAN.F0.4HF  )  GO  TO  14  R  29 

FPPOR  ON  CAPO  R  30 

PRINT  25  R  31 

IAR0RT=1  R  32 

GO  TO  1  R  33 

R  34 

TO  OPERATING  TIME  R  35 

R  36 

no  4  1*1 .NSSYS  R  37 

no  4  Jsl.NPHASES.fl  R  38 

PRINT  27,  T.J  R  39 

CALL  oEAOCn  (CAROS)  R  40 

DO  4  L=1 , A  R  41 

K=L*J-1  R  42 

IF  (K.GT.NPHASES)  GO  TO  4  R  43 

TO(I«K)=CARDS(L>  R  44 

CONTINUE  R  45 

GO  TO  1  R  46 

R  47 

PA  PROAARILTTY  OF  MISSION  ABORT  R  48 

R  49 

no  4  T=1 .NSSYS  R  SO 
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subroutine  input 


00  ft  J=l,NPHASES*R 

P 

SI 

PRINT  2B«  ItJ 

B 

S2 

CALL  PEADCO  (CAROS) 

R 

S3 

00  6  1=1 

B 

S4 

K=L*J-1 

B 

S5 

IF  (K.GT.NPHASES)  GO  TO  6 

B 

S6 

PA ( I »K) *CAPDS (L ) 

B 

S7 

ft 

CONTINUE 

B 

SB 

GO  TO  1 

B 

S9 

c 

B 

ftO 

c 

TN01  RFOUNOANT  SUBSYSTEM  INDICATOR 

B 

ftl 

c 

B 

ft2 

7 

00  fl  I=1»NSSYS»8 

B 

ft3 

PRINT  29,  I 

B 

ft4 

CALL  PEAOCD  (CAROS) 

R 

ft5 

00  R  J=1*R 

B 

ft  6 

K= I ♦ J- 1 

B 

ft7 

IF  (K.GT.NSSYS)  GO  TO  8 

B 

ftB 

TN01(K)*CARDS(J) 

B 

ft9 

R 

CONTINUE 

B 

70 

GO  TO  l 

R 

71 

C 

B 

72 

C 

INO  INDICATORS  OF  PROGRAM  OPTIONS 

B 

73 

C 

B 

74 

9 

CALL  REA0C01  (CAPOS) 

B 

75 

IN02=CAR0S(1) 

B 

76 

IND1=CAR0S (?) 

B 

77 

TS*CAPOS (3) 

B 

78 

DELTA=CAR0S(4) 

R 

79 

XLIMIT=CARDS(S) 

B 

RO 

PRINT  2ft,  IN02,IN03,TS«OELTA,XLIMIT 

B 

Rl 

GO  TO  1 

R 

R2 

C 

B 

P3 

C 

NR  NUMBFR  OF  RFOUNOANT  SUBSYSTEM*;  /SUBSYSTEM 

B 

R4 

C 

B 

R5 

10 

00  11  I=1,NSSYS 

B 

R6 

PRINT  30,  I 

R 

R7 

CALL  PEAOCD  (CAROS) 

B 

R8 

NP(I,1)=CARDS(1) 

B 

R9 

11 

CONTINUE 

B 

90 

GO  TO  1 

B 

91 

C 

B 

92 

f 

P  OPERATIVE  OR  STANDBY  REDUNDANCY 

B 

93 

c 

B 

94 

1? 

00  13  1*1 *NSSYS 

B 

95 

PRINT  31,  I 

B 

96 

CALL  PEAOCD  (CAROS) 

B 

97 

R ( I ♦ 1 ) *CAODS ( 1 ) 

B 

98 

13 

CONTINUE 

B 

99 

GO  TO  1 

B 

100 
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SUBROUTINE  INPUT 


C  R  lnl 

C  FP  MTRF  FOR  RFOUNOANT  SUBSYSTEMS  R  lop 

C  B  103 

IB  00  15  1*1 .NSSVS  R  10B 

NR0=NP(I»1)  R  105 

PRINT  32,  I  R  106 

CALL  PEAOCO  (CAROS)  R  107 

00  15  J*1.NR0  R  108 

FRd.lt  J)=CARnS<J)  R  109 

15  CONTINUE  R  HO 

60  TO  1  R  111 

C  R  112 

C  F  MTRF  FOP  NON-REOUNDANT  SURSYStEMS  R  1 1 3 

C  R  114 

16  no  17  1*1  .NSSYS  R  US 

PRINT  33.  I  R  116 

CALL  PEAOCO  (CAROS)  fl  1 1 7 

F ( I , 1 ) *CARDS (1 )  R  118 

17  CONTINUE  B  119 

GO  TO  1  R  1?0 

18  IF  (IAR0RT.E0.1)  CALL  EXIT  B  1?1 

RETURN  R  122 

C  R  123 

C  R  124 

19  FORMAT  (RF10.0)  R  1?5 

20  FORMAT  (1H1.SX.23HNUMBER  OF  SUBSYSTEM*;  *  •  I5.5X.26HNUMBEP  OF  PHASE  R  126 

1 S/M  I SSI ON  *,IS)  R  127 

21  FORMAT  (8A10)  R  1?8 

2?  FORMAT  </.17H  SUBSYSTEM  NAMES/. 5(10(3X.Al0)/> )  B  129 

23  FORMAT  (A4)  B  130 

24  FORMAT  (/♦  17H  NFW  ARRAY,  10*  .A4)  B  131 

25  FORMAT  (/RRH  ERROR  ON  ARRAY  TYPE  CARO.  WILL  CONTINUE  READING  PARAM  R  132 

1FTFR  DECK  BEFORE  ABORTING  THE  JOB.)  B  133 

26  FORMAT  («M  IND2  =  .I2.9H  ,IND3  *  »I2»l4H  .SUBSYSTEM  *  .12, HH  .INC  R  134 

1PEMENT  *  .F6.0.12H  .MAXIMUM  *  ,F6.0>  R  135 

27  FORMAT  (4M  TO ( • 12. lH, . I2.4H)  *  )  R  136 

2R  FORMAT  (4R  PA ( . 12, lH. • I2.4H)  *  )  R  137 

29  FORMAT  <6H  IN0K.I2.4H)  *  )  B  138 

30  FORMAT  (4H  NP(,I2.6H,1)  *  )  B  139 

31  FORMAT  (3H  R(.I2.6H.l)  s  )  B  UO 

32  FORMAT  (4H  FR ( « I2.8H. 1.1)  *  )  R  141 

33  FORMAT  (3H  F(,I2.6H.l)  *  )  R  142 

FNP  B  163- 
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U  U  >£)  N  or  O' 


SURROUTINF  READCn  (CARDS) 


1 

? 

C 

3 


4 


5 


<iURWOUTINF  PPADCH  (CARDS) 

DIMENSION  CARDS(«).  TCAPD (R) 

RATA  (IRslOM  ) 

NOP=0 
GO  TO  1 
ENTRY  PFAnrOl 
NOP«l 

READ  6«  ICAPD 
IF  (FOF(<S0))  2.3 
PRINT  7 
CALI.  FXIT 

CHFCK  TO  SFF  WHICH  IS  THE  LAST  NON-BUNK  WORD 
DO  4  1=1 

IF  (ICARD(I).EO.TB)  GO  TO  4 
IT*  I 

CARDS ( I ) =0. 

NC*IT*10 

DFCOOE  (NP.8.ICAR0)  (CARDS ( I ) • 1*1 •  IT) 

IF  (NOP.FO.l)  GO  TO  5 
PRINT  9.  (CAPOS (T)tI*l*IT) 

CONTINUE 

RETURN 


FORMAT  (8A10) 

FORMAT  (//» 61H  FND-OF-FTLF  READ  INSTcAO  flF  PARAMETER  CARD.  JOB  AB 
10RTFD....) 

FORMAT  (HFIO.P) 

FORMAT  (lH*»20X*flF10.2) 

END 


C  1 
C  2 
C  3 
C  4 
C  5 
C  6 
C  7 
C  8 
C  9 
C  10 
C  11 
C  12 
C  13 
C  14 
C  15 
C  16 
C  17 
C  18 
C  19 
C  20 
C  PI 
C  22 
C  ?3 
C  24 
C  25 
C  26 
C  27 
C  28 
C  29 
C  30- 
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FUNCTION  FI  (1«JI 


FUNCTION  FI  ( I« J)  0  1 

COMMON  NSFVS*NPHASES.IN02.TN03»IS«0ELTA*X|  IMIT»TO<50»20)*PAr50*20)  0  2 

1 • NAMES (50) • INOl (SO) «NR (SO* 1 ) «R (50* 1) «f (5ft . 1) »►  R <50» 1 *6> *TP (FO) tP  US  0  3 

20*1 > *PMC(FO)  0  4 

F1*EXP<-TP(I)/F(T*J)»  n  5 

RFTUPH  0  6 

FNO  0  7- 
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FUNCTION  F?.  (l.J) 


FUNCTION  P?  <I.J)  F  1 

COMMON  NSSYS.NPHASES*  IN02.  IND3.  IS.DFLtA.XI,  IMIT  .TO (5C .20)  .PA  (50*20)  F  2 

1. NAMES  (SO).  INDUS')). NR  (SO#  1>  .P  (SO.  1  >  .r  (SO,  1  >  *»•» <50. 1  *6)  .TP  (SO)  »P (5  f  J 

20.1) .PMC (SO)  F  4 

TFMPal.  E  5 

NPOaNR(I.J)  F  6 

HO  I  L=1.NR0  F  7 

1  TEMPsTEMP*(l.-EXP(-TP(I)/FR(T,J,L)))  E  8 

F?*).-TFMP  F  9 

RFTURM  E  10 

FNO  F  II* 
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FUNCTION  F3  r  T • 


FUNCTION  F*1  (Ifj)  F  1 

COMMON  NSSVS»MPHASES*IN02*INn3.IS»DFLTA*XI IMIT»TO(SO«20>.PA(50*20)  F  2 

1. NAMES  (SO)  .IN01(SO)»NR<50*1)  . R (50 . 1) .r (5ft ♦ 1 > .FR (50 ♦ l «6> • TP (FO > *P  (5  F  3 

?0*|) #^MC(S0i  F  u 

TFMP*TP(r)/FP(I*J,l)  F  5 

T1=FXP(-TFMP)  F  6 

SUM=0 •  F  7 

NPO=NP(I*J)  F  8 

00  1  L=l*NRO  F  9 

I  SUM=SUM*(  (TEMP)##(L-1I)/IFAC(L-I)  F  10 

F3=T1*SUM  F  11 

PFTURN  F  12 

END  F  13> 
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FUNCTION  IFAr  <I> 


FUNCTTON  TF AC  U) 

C  COMPUTES  THE  FACTORIAL 

ITFMPsl 

!F  (I.EO.O)  1.2 

1  IFAC=1 
RETURN 

2  no  3  K«l,r 

3  !TFMP=ITFMP*K 
IFACsTTEMO 

RETURN 

FNO 


r,  l 

o  2 

g  3 
G  4 
G  5 
G  6 
G  7 
G  8 
G  9 
G  10 
G  11- 
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FUNCTION  65  (I.  |«M*XX> 


FUNCTION  65  <I*J*M«XX)  H 

COMMON  NS6YS«N°HASES« IND2* IND3* IS* DELTA* *1 IMIT*TO(50*20)«PAf50«?0)  H 
1  ♦  NAMES  (501*1 NO 1(50)«NP( 50 *l)*R(50*l)«r(  50  «l)*LR( 50*1*6)* TP (cO)  *PI5  H 
20*1) «PMC<50)  M 

If  (M.EQ.ll  60  TO  3  H 

$UM=0,  H 

DO  ?  L*1*M  H 

SUM 1*1*  H 

DO  i  K=1 .M  h 

I"  (K.EQ.L)  GO  TO  1  H 

SUM1=SUM1*(FP(T*J*L)-FR(I*J.K) )  H 

1  CONTINUF  H 

2  SUMs  SUM  *FR  (I#J*L)##(M»l)*t'XP(“XX/rN(T.J«L))  /SUM1  H 

G5*SUM  H 

RETURN  H 

3  GS*FXP(-XX/FP<I.J.l) )  H 

RETURN  H 

FND  H 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
18- 
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FUNCTION  FS  fI*J> 


V 


FUNCTION  FS  ( I « J)  !  1 

COMMON  NSSYS«NPHASE5*  IN02* TN03* JS*OELTA*X|  IMIT  » TO (50*20) * PA (50*20 ) 
l«NAMFe(SO)  *  INOl  (SO )  *NR  (SO*  l  >  «R (50*  l )  *F (SO*  1 )  *FR(50*  1  *6)  »TP (50)  *P  (5 
20*1 ) *PMC(S0) 

MsNW(I.J) 

XX*TP ( I ) 

F5*G5( I *J*M»XX) 

RETURN 
FNO 
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SUHPOUTINF  TSOPT  <A,N> 


SUPPO'JTINF  TSOPT  (A»N) 

DTMFNSION  A<1)#  TUIM*  10(16) 

f  *1 

J*N 

M*0 

IF  (J.LF.n  GO  TO  9 

IJ*(I*J>/2 

Ks| 

L»J 

IF  (A(I).LF.A(J))  GO  TO  3 
T*A ( J) 

A ( J) *A IT) 

A 1 1 ) *T 
T*A ( I J) 

IF  (A(I).LF.T)  GO  TO  4 
AIT.J)=AII) 

A | I ) *T 
TeA(lJ) 

GO  TO  5 

IF  (T.LF.A(J)l  GO  TO  5 
AdJ)cA(J) 

A(J)=T 

T*A(IJ) 

L*L-1 

IF  (T.LT.ML))  GO  TO  5 
TT=A(L> 

K*K*1 

IF  (A(K).lT.T)  GO  TO  6 
IF  (L.LT.K)  GO  TO  7 
A  (L )  *A  (K  ) 

A(K)*TT 
GO  TO  5 

M*M#1 

IF  (L-I.LP.J-K)  GO  TO  8 

IL<M>=I 

IU(M)=L 

I*K 

GO  TO  10 
IL(M)=K 
IU(M)rJ 
J=L 

GO  TO  10 

IF  (M.FO.O)  PETUPN 
I*IL<M) 

J*IU (M) 

MsM-1 

IF  (J-I.GF.13)  GO  TO  2 
IF  (I.EO.l)  GO  TO  1 
1*1*1  * 

IF  (J.LT.T)  GO  TO  9 


J  1 
I  2 
J  3 
J  4 
J  5 
J  6 
J  7 
J  8 
J  9 
.1  10 
1  11 
J  12 
J  13 
J  14 
J  15 
J  16 
J  17 
J  18 
J  19 
J  20 
J  21 
J  22 
J  23 
J  24 
J  25 
J  26 
J  27 
J  28 
J  29 
J  30 
J  31 
J  32 
I  33 
J  34 
J  35 
J  36 
J  37 
J  38 
J  39 
J  40 
J  41 
J  42 
J  43 
J  44 
J  45 
J  46 
J  47 
J  48 
J  49 
J  50 
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SUBROUTINE  TSORT  <A,N> 


T*A  <  T  >  J  SI 

IF  <ACI-l>.Le.T>  60  TO  11  J  S2 

K*!-l  J  S3 

12  a<K*1)*A(K>  J  S4 

K*K-1  J  R5 

IF  (T.LT.A(K))  60  TO  12  J  S6 

A(K*1)*T  J  S7 

60  TO  11  J  S8 

END  J  S9- 
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SUBROUTINE  CONTROL 


7 


SUBROUTINE  CONTROL  K 

COMMON  NSSYS.NPHASES* IND2*IND3*IS*DELtA«XLIMIT*T0 (50*20) *PA (50*20)  K 
1 *NAMES (50 )*IN01(50)*NR(50*1)*R(50*1)«P(50«1)*FR(50*1*6) «TP(50) » P (5  K 
?0*1)*PMC<50)  K 

DIMENSION  XTFMP(50) «  ITEMP(SO)  K 

K 

PART  I.  CALCULATE  P  K 

K 

00  5  1*1 *NSSYS  * 

TH!)*0.  K 

00  I  J*1«NPHASES  K 

TP( I ) *TP ( I ) *T0( I ♦ J) *PA( I  * J)  K 

IF  (INDl(T).EO.l)  60  TO  2  K 

K 

N0N-RE0UN0ANT  SUBSYSTEM  K 

K 

P(I«1)*F1 (1*1)  K 

GO  TO  5  K 

IF  (R(I.l).EQ.l.)  GO  TO  3  K 

K 

OPERATIVE  RE0UN0ANCY  K 

K 

P(I*1)*F 2(1*1)  K 

GO  TO  5  K 

K 

STANDBY  REOUNOANCY  K 

K 

IF  (FR(I*1*1).NE.FR(T*1«2))  GO  TO  4  K 

K 

EQUAL  MTBF  K 

K 

P(I ♦ 1 ) *F3 (1*1)  K 

GO  TO  5  K 

K 

UNEQUAL  MTBF  K 

K 

P(I*1)*F5(I*1)  K 

CONTINUE  K 

K 

PART  II.  MCSP  K 

K 

PMC ( 1 ) *1 .  K 

00  4  I*1.NSSYS  K 

PMC ( 1 ) *PMC (1)*P(I*1)  K 

K 

PART  III.  RANK  ABORT  CAUSING  SUBSYSTEMS  K 

K 

00  7  1*1 »NSSYS  K 

XTFMP(I)*P(I,1)  K 

CALL  TSORT  (XTEMP.NSSYS)  K 


1 

2 

3 

4 

5 

6 
7 

e 

9 

10 

11 

12 

13 

14 

15 

16 
17 
IB 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 
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SUBROUTINE  COMTROL 


1 


no  M  T*1 »N5SY5 

K 

91 

00  «  Jal.NSSYS 

K 

52 

IF  (XTEMP(I).NE.P(J.l))  60  TO  P 

K 

53 

TTFMP ( T )  =J 

K 

B4 

a 

CONTINUE 

K 

55 

00  9  I=1,NSSYS 

K 

56 

9 

P(I*1)=1-XTFMP(I) 

K 

57 

PRTNT  15*  PMCU) 

K 

58 

IF  (IN02.FQ.0I  GO  TO  11 

K 

59 

PRTNT  16 

K 

60 

00  10  1*1 »NSST S 

K 

61 

JsTTFNPm 

K 

62 

PRTNT  17,  IfNAtfEMJ)  *J*PU»1) 

K 

63 

10 

CONTINUE 

K 

64 

c 

K 

65 

c 

PART  TV.  SESTT1VITY  ANALYSTS  OF  NON-RF0UNDANT  SUBSYSTEM,  IS 

K 

66 

11 

IF  (IN03.F0.0)  GO  TO  14 

K 

67 

MAXsXLIMTT/DFLTAM. 

K 

68 

FS*F(TS,1) 

K 

69 

TPS*TD (IS) 

K 

70 

DO  12  1*2, MAX 

K 

71 

0ELT*(I-1)*DFLTA 

K 

72 

12 

PMC ( I )  =PMC <11  *FXP (-TPS* <  1 ./ (FS^DELT)-l  ./FO  > 

K 

73 

PRTNT  IB,  NAMFS(IS),FS 

K 

74 

00  13  1*1, MAX 

K 

75 

DELT»(I-1)*DFLTA 

K 

76 

PRINT  19,  DELT,PMC(T> 

K 

77 

13 

CONTINUE 

K 

78 

14 

CONTINUE 

K 

79 

RETURN 

K 

AO 

c 

K 

A 1 

c 

K 

A2 

15 

FORMAT  (1H1,17X*7HMCSP  *  ,F8.6> 

K 

P3 

16 

FORMAT  <//,12X,9MSUBSYSTEM,5X»9HSUBSYsTEM,5X*nHPR08ABILITY,/f3X*4 

K 

84 

lHRANK,7X»4HNAME*9Xf6HNUMBER,9X,BH0F  ABORT » 

K 

P5 

17 

FORMAT  (4X,I?,6X,A10«6X,I2,UX«Fe.6) 

K 

86 

la 

FORMAT  (///,2X,3?HNAME  OF  SUBSYSTEM  INCREMENTED  *  « A10t/,5X,22HINI 

K 

87 

1TIAL  MTBF (HOURS)  *  ,F10*2//, 2X,9HTNCRFMENT*6X,4MMCSP> 

K 

88 

19 

FORMAT  (2X,FB.1,5X,FA.6) 

K 

89 

FNO 

K 

90- 
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APPENDIX  B 


DESIGNING  TO  SYSTEM  PERFORMANCE/COST  (DSPC)  COMPUTER  PROGRAM 

DEFINITIONS  OF  INPUTS 

NSYS  =  Total  number  of  systems  (e.g.,  fleet  size). 

NSSYS  =  Total  number  of  subsystems  (NSSYS  <_  40) .  Each  subsystem  is 

identified  by  number  and  name. 

NPHASES  =  Number  of  mission  phases  (NPHASES  _<  20). 

LASTP  =  Phase  through  which  MCSP  is  to  be  calculated  (this  is  usually 

NPHASES  or  the  target  phase). 

NYEARS  =  Number  of  years  to  be  considered  in  calculating  logistic 
support  costs  (e.g.,  system  lifetime). 

NMPM  h  Average  number  of  missions  per  month  per  system. 

TR(i)  s  Ratio  of  total  operating  time  to  mission  operating  time  of 

the  i-th  subsystem. 

T0(1 ,  j)  =  Operating  time  of  the  i-th  subsystem  during  the  j-th  mission 
phase  (i  =  1,2,  ....  NSSYS;  j  =  1,  2,  ....  NPHASES). 

PA(i,  j)  =  Conditional  probability  of  mission  abort  given  that  the  i-th 
subsystem  has  a  failure  during  the  j-th  mission  phase. 

N(i)  =  Number  of  nonredundant  options  (other  than  the  baseline 

subsystem)  for  the  i-th  subsystem  (N(i)  ^5). 

R0(i)  =  Number  of  redundancy  options  (other  than  the  baseline)  for 

the  i-th  subsystem  ( RO ( i )  <_  5). 

IND1(1)  =  Indicator  equal  to  0  or  1.  If  INDl(i)  =  0  then  the  baseline 

for  the  i-th  subsystem  is  nonredundant.  If  INDl(i)  =  1  then 
the  baseline  for  the  1-th  subsystem  is  redundant. 
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NR(i,  j)  Number  of  redundant  units  for  the  j-th  redundancy  option  for 

the  1-th  subsystem  (j  f.  R0(  1 )  +  1 ). 

R(i,  j)  =  Indicator  equal  to  0  or  I.  If  R(1,  j)  =  0  then  the  j-th 

redundancy  option  for  subsystem  i  has  operative  redundancy. 

If  R(1,  j)  =  1  then  the  j-th  option  is  standby  redundant 
(j  <_  NR(i ,  j)  +  1). 

F(i,  j)  =  Mean  operating  time  between  failure  for  the  j-th  nonredundant 
option  for  subsystem  i  (j  <.  N(i)  +  1). 

UC(i,  j)  =  Unit  acquisition  cost  of  the  j-th  nonredundant  option  for 
subsystem  i  (j  <  N(i)  +  1). 

CR(i,  j)  =  Average  cost  per  repair  of  the  j-th  nonredundant  option  for 
subsystem  i  (j  <.  N(1)  +  1). 

FR(i,  j,  k)  =  Mean  operating  time  between  failure  of  the  k-th  redundant 
subsystem  of  the  j-th  redundancy  option  for  subsystem  i 
(j  <  R0(i)  +  lj  k  <  NR(i,  j)). 

UCR(i,  j,  k)  =  Unit  acquisition  cost  of  the  k-th  redundant  subsystem  of 
the  j-th  redundancy  option  for  subsystem  i  (j  £  R0(i)  +  1, 
k  <  NR(i,  j)). 

CRR(i,  j »  k)  =  Average  cost  per  repair  of  the  k-th  redundant  subsystem  of 
the  j-th  redundancy  option  for  subsystem  i  (j  <  R0(i)  +  1, 
k  <  NR(i,  j)). 

2.  OUTPUTS  OF  THE  DSPC  MODEL 

The  model  outputs  are  printed  in  two  tables.  The  Adjusted  Baseline  System 
is  printed  first  to  define  those  options  which  lead  to  a  higher  MCSP  at  lower 
cost  (this  results  from  the  reliability  management  procedure  described  in 
Section  III).  The  Adjusted  Baseline  System  defines  the  starting  point  for 
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the  optimization  procedure.  The  form  of  the  printout  for  the  Adjusted 
Baseline  System  is: 

ADJUSTED  BASELINE  SYSTEM 

Subsystem  Subsystem  Option  for  Adjusted 
Number  Name  Baseline _ 

1 

0 

2 

NSSYS 

The  meaning  of  the  first  two  columns  is  self-evident.  If,  for  subsystem  i, 
the  number  z  appears  in  column  3  this  means  that  nonredundant  subsystem  with 
F(1,  z),  UC(i,  Jt),  and  CR(i,  z)  should  replace  the  baseline  for  subsystem  i. 
If  HR  appears  In  column  3,  then  the  baseline  is  replaced  by  that  redundancy 
option  corresponding  to  FR(i,  zt  k),  UC(1,  l,  k),  and  CRR ( i ,  z,  k)  where 
k  =  1,  2,  ....  NR(i,  «,).  If  i  =  1  appears,  then  the  baseline  is  the  best 
starting  point. 

The  second  table  presents  the  Optimal  Subsystem  Options  in  the  following 
form: 


Configuration 

OPTIMAL  SUBSYSTEM  OPTIONS 

MCSP  Acquisition  Logistic  Total 

Subsystem 

Option 

Identification 

Cost 

Support  Cost 

Clanged 

Selected 

_ £cn _ 

Cost 

Baseline  Cl  =  1 


Adjusted 
Baseline  Cl  =  2 

Cl  =  3 
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The  configuration  identification  merely  numbers  the  sequence  of  optimiza¬ 
tion  steps.  For  the  adjusted  baseline  system  the  corresponding  options  were 
defined  and  printed.  For  each  configuration  Identification  (after  the  baseline) 
the  last  two  columns  define  the  subsystem  changed  and  the  option  selected  for 
that  subsystem. 

3.  OSPC  CARD  INPUTS 

A  description  of  the  Input  cards  for  the  DSPC  model  Is  presented  in  this 
section.  As  mentioned  In  Appendix  A,  the  form  of  the  card  Inputs  to  the  MCSP 
model  are  identical  to  those  for  the  DSPC  model  with  the  few  exceptions 
described  previously.  It  must  be  pointed  out  that  on  all  numeric  cards 
each  value  must  be  followed  by  a  decimal  point;  on  alphanumeric  cards  no 
decimal  point  Is  allowed, 
a.  First  Card. 


FIO.D 


INPUT 


NSSYS 

NPHASES 


LASTP 


b.  Second  Card(s). 


NYEARS 


FIELD 


INPUT 


Name  of  subsystem  1 
Name  of  subsystem  2 


Name  of  subsystem  8 


Input  as  many  cards  required  to  name  all  NSSYS  subsystems.  The  number 


of  cards  required  Is 


where  the  notation  [y]+  denotes  the  smallest  integer  greater  than  or  equal 
to  y. 

c.  Cards  for  the  One  Dimensional  Arrays. 

One  dimensional  arrays  are  required  for  the  inputs  TR(i),  N(i),  R0(i), 
and  INDl(i).  For  each  of  these  inputs  i  *  1,  2,  ...»  NSSYS.  The  first  card 
of  a  one  dimensional  array  contains  the  array  mnemonic  beginning  in  column  1, 
i.e.,  starting  in  column  1  one  of  the  mnemonics  TR,  N,  RO,  or  IND1  is  printed. 
The  next  cards  for  each  mnemonic  are  as  follows: 


FIELD 


8 


INPUT 

value  corresponding  to  subsystem  1 
value  corresponding  to  subsystem  2 

t 

value  corresponding  to  subsystem  8 


Continue  until  all  values  are  defined  for  each  mnemonic.  The  procedure  is 
repeated  for  each  of  the  four  one  dimensional  arrays, 
d.  Cards  for  the  Two  Dimensional  Arrays. 

Two  dimensional  arrays  are  required  for  the  inputs  T0(i,  j),  PA(i,  j), 
NR { i ,  j),  R(i,  j),  F(i,  j),  UC(i,  j),  and  CR(i,  j).  The  first  card  of  any  two 
dimensional  array  contains  the  array  mnemonic  (TO,  PA,  NR,  R,  F,  UC,  or  CR) 
beginning  in  column  1.  For  example,  after  the  mnemonic  TO  the  next  set  of 
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cards  (corresponding  to  subsystem  1)  is  as  follows: 


FIELD  INPUT 

1  T0(1,1) 

2  T0(1,  2) 

•  • 

•  • 

•  • 

8  T0(1.  8) 

Continue  to  input  T0(1,  j)  until  j  reaches  its  maximum  value  (for  the  i-th 
subsystem  the  maximum  value  of  j  for  TO  and  PA  is  j  =  NPHASES;  for  F,  UC,  and 
CR  the  maximum  value  of  j  is  N(i)  +  1  £  6;  for  NR  and  R  the  maximum  value  of 
j  is  R0(i)  +  1).  The  next  set  of  cards  for  subsystem  2  are: 

FIELD  INPUT 

1  T0(2,  1) 

2  T0(2,  2) 

•  • 

•  • 

•  • 

8  T0(2,  8) 

Continue  to  input  T0(2,  j)  until  j  reaches  its  maximum  value.  Continue  the 
process  until  the  values  of  T0(i,  j)  are  input  for  i  =  1,  2,  NSSYS. 

The  process  is  repeated  for  each  two  dimensional  array  corresponding 
to  the  mnemonics  TO,  PA,  NR,  R,  F,  UC,  and  CR. 
e.  Cards  for  the  Three  Dimensional  Arrays. 

Three  dimensional  arrays  are  required  for  FR(i,  j,  k),  UCR(i,  j,  k), 
and  CRR(i,  j,  k).  For  each  of  these,  i  runs  from  1  to  NSSYS,  j  from  1  to 
R(i  )  +  1  <.  6,  and  k  from  1  to  NR(i,  j)  5.  As  before,  the  first  card 
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contains  the  mnemonic  beginning  in  column  1.  For  example,  after  the  card 
containing  the  mnemonic  FR  the  cards  are  as  follows: 


FIELD  INPUT 

1  FR(1.  1,  1) 

2  FR(1 .1.2) 


NR(1,  1)  FR(1,  1,  NRO,  1)) 

1  FR(1,  2,1) 

2  FR(1,  2,  2) 

2nd  Card  •  • 


R(1 )  +  1  Card 


NR(1 ,  2)  FRO,  2,  NRO,  2)) 

1  FR(1,  R(1 )  +  1,  1) 

2  FR(1,  RO)  +  1,  1) 


NRO  ,  R(1 )  +  1)  FR(1,  RO)  +  1,  NR ( 1 ,  R(l)  +  1)) 


Repeat  the  same  procedure  for  each  subsystem  where  the  i-th  subsystem  consists 
of  R(i)  +  1  <.6  cards. 
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r>  n 


4  -  DSPC  PROGRAM  LISTING 


PPOG«»M  OSPC  < INPUT •OUTPUT « T APgftO* INPUT) 

INPUT  THF  APPAYS 
CALL  INPUT 
C  NOW  FXFCUTE  THE  MAIN  ROOY 

CALL  CONTOOL 

C  THATS  ALL  FOLKS . 

FNO 


A  1 
A  ? 

A  3 
A  4 
A  5 
A  6 
A  7 
A  8- 


inn  i  v  —  n  r>  -> 


SliRPOIJTTNC  INPUT  R  1 

rOMMON  NfrYS.NSSYS.NPHASrfr.LASTP.NYFAPS.NMPM.TkUO)  .T0(40.20)  »PA(4()  R  2 

1 .20)  .N<40)  .ROUO)  •  TNtni  (40)  , NP  (40.6)  »P (40.fr )  «F(40.fr)  «UC(40*6)  *CB(40  R  1 

2.6) ,F® (4». fr.fr )  ,ijr9(4<1.fr.6)  ,rpp  (40 .fr.fr)  «NAwES(40) «T (40) *TP<40)  »TNP(  H  4 

1<*0)  »P (4rt.fr)  .nCA(40.fr)  .UCSUO.fr)  .PMC (4*)  .CA(40)  .CS(40) .C<40>  .PR (40*  B  5 

44)  ,, ICAR  U0.fr)  .IlCfrO  (40.6)  .HC°P  U0.fr)  »X  (40.20)  »  IX  (40.20)  *Y  (40.20)  .2  <  P  fr 

b4 0*20)  .MY  (40)  .  TASI.  (40)  oKS(40)  .YS(40)  »7S(40)  *XL(40)  .LAMBDA  (40)  *  I D 1  (  «  7 


640) «S^C(40I .frnfrf4rt) .TC(40»fr)  R  8 

•MMENSION  rA»0S<8)  r  9 

R  TO 

r  Uf  A,)  THE  FIRST  CAPO  H  U 

C  P  )2 

RF AO  12.  CARns  r  ,3 

Mfr YS-FAPOc ( ] ) $NSSYS=C ARDS ( 2 )  SNPH A5ES=rAPDS (3)  R  14 

l  AfrTP=CARnfr(4)  *NYFAPS=CAPnS(S)SNMPM=r.y»OS(6>  R  )S 

PPTNT  31.  MSYS.NfrSYfr.NPHASFS.LASTP.NYrARS.NMPM  R  )6 

PfAD  35.  (NAMFS(T).T=1.NSSYS)  r  ,7 

PPTMT  16.  (NAMES  (T) « T®1 .NSSYS)  P  18 

B  19 

UFOIN  REAiJNC  THF  ARRAYS  R  30 

B  21 

PF  AO  17.  TAN  R  2? 

IF  (EOF (60) )  11.2  R  ?3 

PRINT  14.  I  AM  R  24 

IF  (IAN.P J.4HTR  )  00  TO  1  R  25 

IF  (IAN.FP.4HTO  )  (,n  TO  S  R  26 

IF  (IAN.E0.4HPA  )  r-n  TO  7  R  27 

IF  (TAM.F0.4MN  )  00  TO  9  B  28 

IF  (TAN.FO.4HRO  )  00  TO  U  R  29 

IF  (IAN.F0.4HINDT)  00  TO  11  R  30 

IF  (IAN.F0.4HNR  )  00  TO  IS  R  31 

IF  ( I  AN.FO.4HR  )  OO  TO  17  R  32 

IF  (TAN.FO.4HF  )  00  TO  19  R  33 

IF  (IAN.F0.4HUC  )  00  TO  21  R  34 

IF  (I'N.FO.AHCR  )  00  TO  21  R  15 

IF  ( I AN.F0.4HFR  )  fiO  TO  ?5  R  36 

IF  (TAN.F0.4HMCR  )  00  TO  27  R  37 

IF  ( IAN.F0.4HCRR  )  00  TO  29  R  3ft 

FPPOP  ON  CARP  R  39 

PRTNT  12  R  40 

00  TO  1  R  41 

R  42 

TR  MTSSTON  OPFPATING  TImE/TOTAL  pPEPaTING  TIME  «  43 

B  44 

3  00  4  T=1 .MSSYS.P  R  45 

PPTMT  19.  T  R  46 

CAI  !.  RFAOCO  (CAPOS)  R  47 

00  4  J=l.°  P  48 

K=T ♦ J-l  R  49 

IF  (K.OT.MSSVS)  00  TO  4  R  SO 
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nnn  i-  -oonn  »  ^iooo>  jiooo 


5USP0UT  1  '••fc  fM-’tlT 

TR(K)rCAPn$( J)  R  SI 

4  CONTINUE  B  S2 

GO  TO  1  H  S3 

R  R4 

TO  OPERATING  TTmE  8  S5 

R  56 

00  S  I«1,MSSYS  R  S7 

00  6  J=l«NPHAScS.fl  R  SB 

'PRINT  40*  I*J  R  SR 

TALL  PFAOro  (CAPOS)  R  SO 

00  6  |.*1  «R  .  R  Si 

K=L*J-1  R  62 

TF  (K.GT.NPHASES)  GO  TO  6  R  S3 

TO(I,N)=CAROS(L>  R  S4 

CONTINUE  R  S5 

GO  TO  1  R  66 

R  67 

PA  PROP  ABILITY  OF  MISSION  ABORT  R  SB 

R  69 

00  1  I=1*MSSYS  R  70 

00  H  J=l.NPHASES»R  R  71 

PRINT  41 «  I.J  R  72 

CALL  RFAOCO  (CAROS)  R  73 

00  S  L=1«R  R  74 

K»L ♦ J“ 1  R  75 

TF  (K.GT.MPHASFS)  GO  TO  8  R  76 

PA(I«K)=C*ROS(L>  R  77 

CONTIVIJF  A  78 

«0  TO  1  R  79 

R  RO 

M  NIJWRFR  OF  NON-REOIJNOANT  OPTIONS  R  Rl 

R  «2 

OO  10  1=1 «NSSY5*R  R  P3 

PRINT  42*  I  R  H4 

CALL  PFAOCD  (CAPOS)  R  R5 

00  10  J= 1 « R  R  R6 

K=T ♦ J-l  R  q 7 

IF  (K.GT.NSSYS)  GO  TO  10  R  R8 

M(K)sCAPOS(J)  R  R9 

CONTINUE  R  90 

GO  TO  1  R  91 

R  92 

PO  NUMRFR  OF  RFniJNDANT  OPTIONS  R  93 

R  94 

11  no  1?  I=l«NSSYS*R  R  95 

PPTNT  41,  I  R  96 

TAIL  PFAnro  (CAPOS)  R  97 

00  l?  J=1»H  R  98 

K=IO-l  R  99 

TF  (K.GT.MSSYS)  GO  TO  12  q  100 
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I 

SUBROUTINE  INPUT 


PO(*)=C AQOS(J)  R  101 

12  CONIJNUF  R  102 

<H)  TO  1  R  103 

T  R  104 

C  IND1  INDICATOR  FOP  BASELINE  RFDUNO*NCV  R  lo5 

C  R  106 

H  00  14  I=1.NSSYS»8  R  107 

PPTNT  44.  I  R  lo8 

cai  i.  PFAorn  (capos)  fl  ln9 

00  14  J=lf«  R  1)0 

K  =  T ♦ J- 1  R  111 

IF  (K.OT.MSSYS)  00  TO  14  R  H2 

?NnnK)=c‘'pn«;(j)  r  l i 3 

14  CONTINUE  H  114 

no  TO  1  R  115 

C  R  116 

C  MP  N'IMPFP  OF  RCOUNOANT  SUBSYSTEM*:  R  H7 

C  B  11B 

15  00  16  I=1,NSSYS  R  119 

NPO=P0 ( T ) ♦ 1  R  120 

PPTNT  4S.  I  R  121 

CALL  RFAOCO  (CAPOS)  R  1?2 

00  16  JsI.NkO  R  123 

NR(l*J)=CA«OS(.|)  R  124 

16  CONTINUE  R  125 

np  TO  1  R  126 

C  R  127 

C  o  OOEPATIVP  OP  STANOPV  PFOUNOANrY  R  12B 

C  R  129 

17  no  18  I=1.NSCYS  R  130 

NPOsRO ( T ) ♦ 1  R  131 

PPTNT  46.  I  R  132 

CAI L  PFAncO  (CAPnS)  R  133 

no  18  J=1  .TiPn  R  134 

P(1.J)=CAOOS( J)  R  135 

i*  continue  r  136 

no  TO  1  R  137 

0  R  138 

C  F  MTRF  FOP  NON-PEOUNOANT  SUBSYSTEMS  R  139 

C  R  140 

19  no  20  1*1 *NSSYS  R  141 

NN'POaV(I)  *1  R  142 

PPTNT  47.  T  3  143 

CALL  0EAOC 0  (CAPOS)  R  144 

no  80  J=l,NNPO  R  145 

C(T.J)=CAPnS(J)  R  146 

20  CONTINUE  R  147 

no  TO  1  R  148 

C  R  149 

C  KC  IIMTT  cnST  F0°  NON-PEDUNOANT  SmBSVCTFMS  R  ISO 


I 
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SUBROUTINE  Tmpot 


r. 

R 

1F1 

?\ 

on  ??  t=i.nssvs 

R 

152 

mnposk1  ( 1 1  ♦  i 

H 

1*3 

PPT'IT  4R.  I 

R 

154 

r/»u.  oFanro  fCARns) 

R 

1S5 

nO  ??.  Jst.MNPO 

B 

1S6 

iirn«j)=c«»n<;(  j» 

B 

1S7 

?> 

rOMTIMJF 

B 

158 

r.o  TO  l 

R 

1  SR 

r 

B 

ISO 

r 

f.p  AVfPflOF  TOST  OF 

PF0ATR 

FOP 

NftM-RFnuNOANT  SUBSYSTEMS 

B 

101 

r 

B 

10? 

no  P4  T*1.M$SYS 

B 

103 

nnpo=n ( n ♦ i 

B 

104 

PRTNT  49.  T 

B 

105 

ran.  PFAorn  tcapoS) 

B 

106 

no  >U  Jsl.NNPO 

B 

lo7 

cP(i.j)=rA»ns< ji 

B 

1  *8 

?•* 

TONTTMUF 

B 

1*9 

Oft  TO  1 

B 

170 

r 

B 

171 

C 

FP  MTPF  FOR  Rtft'INOANT  SUBSYSTEMS 

B 

17? 

C. 

B 

173 

?■» 

Oft  ?h  T=1*NSSYS 

B 

174 

NPPsPO ( I  > ♦ 1 

B 

175 

no  ?6  j*  l  «npo 

B 

176 

NRS=N°( T  .  J) 

B 

177 

PRTNT  SO.  T.J 

B 

178 

rail  RFAftfn  <cawns) 

H 

17Q 

no  ?6  K=1,NRS 

B 

mo 

FR(T.J.K) sCAPOS ( K ) 

B 

mi 

continue 

B 

m? 

OO  TO  1 

B 

m3 

C 

R 

JB4 

r 

HOP  IIVTT  COST  FftO  PFOUNftANT 

SUHSYcTFMc; 

R 

185 

C 

R 

1P6 

?7 

no  pm  T  =  1  .NSOYS 

B 

1R7 

MPO  =  PO  < I ) ♦  1 

B 

1R8 

no  ?8  »=i . npo 

R 

189 

NRS=N° (T.J) 

B 

mo 

PRTNT  SI.  T.j 

B 

mi 

rau.  PFAorn  tcaROS) 

B 

19? 

no  ?f\  K=i  .npo 

B 

193 

iirPU,  j,K)=rawoS(K) 

B 

104 

?•* 

roKTT'.UF 

R 

105 

oo  TO  1 

B 

106 

C 

B 

107 

r. 

TRu  AVFPAOP  POST  OF 

RFPMR 

FOR 

PFoUNOANT  subsystems 

fl 

198 

c 

B 

199 

no  TO  T  =  1  .MS^YS 

B 

?no 
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SUHOOUTINF  TMPUT 


MPOaOOdUl 

R 

201 

on  in 

J=l .NRO 

« 

202 

nP«;=no(T.  J) 

H 

203 

PRINT 

S?.  T.  1 

R 

204 

tail  PSAnni  (CA»nsi 

R 

2  ns 

no  30 

K=| .MRS 

R 

206 

cpp ( I « 

J.K)=CARDS(K> 

R 

207 

30 

font  IN 

OF 

R 

208 

on  to 

1 

R 

209 

r 

THAT5 

ALL  OF  THE  INPUT  CAROS 

H 

210 

31 

CONTINUE 

R 

211 

PFTURN 

R 

212 

c 

R 

213 

c 

R 

214 

3-* 

format 

(8FI0.01 

R 

215 

3  1 

FORMAT 

C 1 4H 1  FLFFT  SI7E  =.TS,5X#??HNI»mRER  OF  SURSYSTEMS  *»IS«5X*?6 

R 

216 

lhNHMrtFP  OF  PHASFS/MTSSION  S.T5/14H  LAST  PHASt  *» IS*5X* 1RHLTFE  SPA 

R 

217 

?N(YFACS)  a.TS#3X.33HMIJMflFP  OF  MiSSIONS/MOvTH/SYSTFM  *»I5) 

R 

218 

34 

FORMAT 

</«I7H  NFW  ARRAY  10  r  ,A4) 

R 

219 

33 

FORMAT 

(RA 1 0 1 

R 

2*0 

3h 

FORMAT 

1/.16H  SUBSYSTEM  NAMFS/.5 < 1 0 (3y ♦ A  1 0 ) /l ) 

R 

2*1 

37 

FORMAT 

( A4 ) 

R 

2*2 

3rt 

FORMAT 

(//•?6H  FRR0O  ON  ARRAY  TYPF  CARD.) 

R 

2?3 

3<» 

FORMAT 

(4H  TP(*f3*4M)  =  1 

R 

224 

40 

FORMAT 

(4M  TO(«I3»lH** I ?*4H )  *  1 

R 

2*5 

41 

FORMAT 

C4W  PAC,I3.lH.*I2»4H)  *  ) 

8 

2*6 

42 

FORMAT 

(3M  N1.T3.4H)  =  ) 

H 

227 

43 

FORMAT 

(4«  R01.f3.4M1  =  1 

R 

2*8 

44 

FORMAT 

16«  TMD1 |,T3*4H>  s  1 

R 

2*9 

43 

FORMAT 

14“  NP1.I3.7H.  1)  r  ) 

R 

230 

46 

FORMAT 

(3H  R  1 » T3.7H.  1)  =  > 

R 

231 

47 

FORMAT 

(3W  Fl.T3.7H.  1)  =  ) 

R 

232 

43 

FORMAT 

(4H  Ur  1 »  I 3.7H.  1)  =  1 

R 

233 

4*^ 

FO°MAT 

(4M  CP1.T3.7H.  1)  a  1 

R 

234 

so 

FORMAT 

<4H  FP<.n*lM.,I?,7H.  1)  =  ) 

R 

235 

SI 

FORMAT 

1SH  UrR(,T3«lH,*I2.7M,  ])  =  ) 

R 

236 

52 

FORMAT 

1SH  CPP(»T3«lw«t!2»7H*  1)  =  1 

R 

237 

FMO 

R 

238 
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C  M 


1 


3 


r 

r 

4 

K 

ft 

7 


SUPPO'iTT'iF  RFAOC M  (CARDS)  C 
DIMENSION  CAOns(P)  ,  lCARf)(M)  C 
DATA  (mrlOH  )  r 
PFAi)  4*  KARn  c 
?F  (FOF(ftn))  \,?  c 
PRINT  5  C 
CAI  ».  fXTT  C 
CHFCK  TO  SFF  whI C«  TS  THE  LAST  NOK'-BLaN*  wORO  C 
no  1  **!•«  C 
TF  (KAPIHn.FO.TB)  00  TO  1  C 
!T=T  c 

CAPOS(T)=o.  C 

wC=IT*10  n 
nFCODr  (MC,ft.tCAon)  (CAROS( 1) . T*l* IT)  C 
DPTNT  7.  (CAPI)S  (  f)  ♦  T  =  1 » IT  >  C 
OFTMPM  C 

C 

c 

FOP  4  AT  HAJO)  C 
FOB^AT  (//*ftn  TNO-OF-FILF  READ  iNSTrAO  OF  PAkAMFTER  CARD.  JOB  Ap  C 
10PTF0....)  C 
FORMAT  (SF10.P)  C 
FORMAT  (1H*»50X«BF10»?)  C 
FMO  C 


1 

2 

3 

4 

5 

6 
7 
A 
9 

10 

11 

12 

13 

1* 

15 

16 
17 
IB 
19 
*0 
?1 
?2 
?3 
?4- 
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function  n  n.j>  o  1 

COMMON  NSYS.MSSYS.NPHASFS.LASTP.NYFA9«;»NMPM.TWC40).TO<40.20).PA<40  0  2 

l.?0)  .N(40»  .P0(40)  •IMOl(40)  .NP(40.6).9f40.M»F  (40.6)  *UC <40*61  .CP (40  0  3 

2.6)  *FP <40.6.6 ) ,UCO (40.6.6) .CPR (40.6.6)  .NAMES  (40)  *T  (40)  *TP(40)  *TNP(  0  4 

140 )  *P (40.4)  *IICA  (40*6)  *UCS(40«6)  .PMC <4-i)  *CA  (40 )  *CS (40)  *C(40)  .PR (40*  0  5 

46) .l)C*P<40.6)  »UCS9 (40.6) «IJC»P  (40 .6 > .  X  (40*?0 >  »  I X  (40 .20 )  .  Y  <40.20 )  ♦  Z  <  0  6 

S40.?0).MY(40).IAPL(40).XS(40)«VS(40».?S(40).XL(40).LAMBOA(40)»ID1(  0  7 

640) .SSC<40) .POS(40) .TC(40.6)  0  0 

F1*FXP(-T0(I)/F(T.J))  0  9 

PFTIJRN  0  1 0 

FNO  0  )1- 


FIINOT  T0N  F?  H.J*  F  1 

COMMON  N»;v<;.MSSV«;.NPHASFS.l.fl<;TP.NYFAe«;*NMOM.TK<40>#T0<40*?0**PA(40  F  2 

(401  .Pn(^O)  «  TMO|  (401  ,MR(40«M  .Qf40*M  *F  (4Q*M  »UC<40.6>  ♦CRUO  F  3 

*«*)  ,f-c  <40.*.* )  , l irp  (40,4,4)  .COp (40.4.4*  t names  (401  .T(401.TP(40) .TNP(  F  4 

*40)  ,P(40.4|  ,  MCA (40.4)  *UC5  (40.4)  , PMC  (  4o  )  . C  A  ( 40  I  « CS  (  40  )  « C  ( 40  )  , PW  ( 40  *  F  5 
4M  .i)C./'P(<*0.4)  .uCSP  (40*6)  .IIC*P<40.4)  »> (40.P01*  IX <40.201  .Y 140. 201  .7(  F  6 

14O.P0)  .MV (40)  .TAPL  (40)  «<S(40)  «yS(40>  »7S(40) .XL <401  tLAMHOA  C40)  . 101  (  F  7 
*401  .S«C<44)  .c<)S<40)  «TC<40.4)  F  H 

TFMPs] ,  F  9 

MPO=ND  ( I  .  )}  F  10 

no  |  |  =l«MPO  F  ii 

1  TFMP=TFMP»<l,-FXO(-T°U)/FOOfJ«L>n  F  1? 

F?s) ,-TFPn  F  13 

PFTijP^1  F  14 

FNO  F  15- 
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FUNCTION  ri  (f,j»  f 

COMMON  NSYS«NS$YS.NPMASFS.I_ASTP.NYEAPc.NMPM.TM(40>  .T0(40.?0)  «PA(40  F 
l.?0) ,N(40),PO(4rt) (40) .MR (40.6) »R (40.6)  .F (40.6! .UCI40.6) «CR(40  F 
P.6)  .FP  (40 .6,4)  .OCR  (40 .6. 6 ,  .CRR<40.6,6)  .NAMES  (40)  ♦  T  (40 )  ♦  TP  (4f> )  .  TMP (  F 
340)  .P(40«6)  ,1/0  4(40.6)  ,UCS(40,6)  .PMC  (A'*)«C4(40).CS(40)  .C(40)  .PH  (40.  F 
4M  .(/CAR  (40. 6)  .i/CSR  (40.6)  .UCRP  ( 40 . 6 )  «  X  (40  .  ?0 )  .  i  X  (40 .20  )  .  Y  (40  .?C  )  ♦  Z  (  F 
S40.?0).MY(40) .TABL(40)«XS(40) .YSI40) .7S(40).XL(40) .LAMBDA (40 ). 101 (  F 
640 ) «  SSC (40 ) .SOS (40) *TC(40«6)  F 

TFMp=TP(n/FC  (!♦  J.l)  F 

Tl=f XO(-TFMP)  F 

St|M=0.  F 

MPOiNP(T.J)  F 

00  1  1=1. MOO  F 

SIJM=SUMM  (TFVP)**(L-l )  )/IFAO(L-l)  F 

F.1*T1*SUM  F 

RFTHP*  F 

FNO  F 
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ft  INC  T  TON  TFAC  (I» 
COMPETES  THF  FACTOWTAL 

TTFMP=1 

IF  (T.FO.ft)  I.? 

IFAOsl 
»FT0«N 
no  3  K=1 . 7 

TTFMP=tTFMP#K 

|FAC=7TFMO 

»FTU«*' 

CNf) 


FUtonTTO*  p*  ( !  *  J)  H 

rOMMu*'  N4V9,nsSV5.NPMASF5.L*STP«NYEAP<*NPP«<*TW<40)  *TO<40*20>  *PA<40  H 
!•?())  **'<40)  *00(40)  tTMOl  <40 )  *NP<40*6)  *P  r 40 *M  »F  <40  » 4 )  *UC<40«6)  *CM(4A  m 
2*6 1  *FO (40*4.4)  *Uro  <40*6*4 1  *CRP(4G.4,4t  ,NAmES (40 ) • T <40 )  *TP<4(')  *TNP(  H 
14  0) «P (40.4) *UCA (40*6) *UCS (40 *6) .PMC (40) *f A <40)*CS  <401 *C <401. P« <40*  H 
44)  *UCA9(4A*6)  *tiCSP<40*6)  «ljf*PP  (40*6)  *X( 40*901  •  IX  (40*20)  *V<40*20) *7<  H 
540*90)  *MY  (40)  *t  ARLUO)  ,XS(40)  « YS (40 ) « 75 ( 4C ) *  XL  <40  )  *LAMHOA ( 40 > • 1 0 )  (  H 


440)  *55040)  .404(40)  *TC<40*4)  M 

HTMPNSION  XY (S)  H 

5| |M?*0  •  H 

T  1  =  1  (  T  )  /FP(  I  •  I*  1 )  w 

T?rC«0<I*.l*l)  M 

T1=rNP(  I)/FP( T.J.l)  H 

WRO=NP(I,J)-J  H 

00  9  L*1,X|P0  H 

l  L=L  * 

M)Ml=o.  4 

no  1  K=1,LL  H 

T4=(T3*«<K-1))/TFAC<K-1)  H 

1  SUM1=S0M1*T4  H 

XY(L)=SUMl*FXP(-n)  M 

2  SUM=SIJM-AL0G(XY(L)  )«CWR(  I.J.LM)  H 

F4=(T?*.4*SUM)*T1  H 

PFTllR\'  H 

FNO  H 


1 

2 

3 

4 

5 

6 
7 
b 
9 

10 

11 

12 

13 

14 

15 

16 
17 
IP 
19 

90 

91 

92 

93 

94 
95- 
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FUNCTION  nt,  (I.J.M.XX)  I  1 

COMMON  NSvsf ^ORYO.NPHASFS.L ASTP.NYFARs.NMDM.TK (40)  *TO(4O.P0)»PA(40  !  2 

1 .?0) »M(40) .90(40) .IMPI (40)  .NR  (40*4)  »Pf40.4 )  ,F (40 .4) *UC (40 .6)  »CP(40  !  3 

2.4) .FR (40.4.4) ,Ur» (40.6.6) .CRR (40.6.41 .NAMES (40) *T(40) *TP(40) .TNP(  I  4 
140)  .p  (41.. 4)  .DC A  (40.4)  »UCS(40.6)  ,PMC(4o)  .CA(40)  «CS(40)  .C(40)  .PR  (40,  T  5 
44) ,DCAP(40.6) .DCRP (40,6)  .IJCWP ( 40 ,4) « X f 40.20) • IX (40.20 ) « V (40,20 > .7 (  T  6 
540.20) *MY(40),(ARL(40).XS(40) • YS (40 ) • 7S (40) . XL (40 )  .LAMBDA (40) ♦ 101 (  T  7 
440).S<;C(40).ROS(40),TC(40.4)  T  8 

TF  (M.FG.l)  Of)  TO  3  T9 

SIIM=0.  !  )0 

00  2  L*1.M  T  11 

5UM1=).  I  IP 

00  1  Kal.M  l  )  3 

IF  (K.FO.L)  00  TO  1  I  )4 

SUM) sOUMl * (FP(I.J.L) “FR (I.J.K))  I  15 

1  CONTINDF  T  14 

?  SUMaSDM4FO(!.J.L)**(M-l)*FXP(-XX/FR(I.J.L))/SUMl  T  17 

05=SUM  I  )B 

RETURN  T  19 

1  OSaEXP(-XX/FP(T,J.l) )  !  20 

RETURN  T  ?1 

FNO  I  ?2- 


FUNCTION  Fq  n.J*  J  1 

COMMON  N<?YS<NSSV«;«NPMASFS«l.ASTO«NyEAeq*NMPM»TM(40I«T0(40*?0>«PA(40  J  2 
1  t?0)  »PO<40)  flNni  (40) «NR(40.A) »Mf40'*)«F(40*A>  .UC140.A1  *CP(40  I  3 

2.  A)  .FP<49,A,A»,IICR(4n»6«6).CPP(40*4t(St  .NAMES (40)  »T<40>  .TP  (40  *TNP(  )  & 

140) .P(40*A) .1 1C,  4(40. A) »UCS (40.6) «PMC(4n)  .CM40)«CS(40)*C(40)«PR(40.  J  5 
4M  ♦UCAPI40.6)  *UCSR(40,6)  .UCPP (40* A)  « X  (40*?0)  ♦  IX (40.20)  •  Y <40 .20)  .Z  (  J  6 
**40.20) *MY(40) .TABL(AO)  ,XS(4rt) ,YS(40) .7S(40).XL(40) .LAMBDA (40) « 101  (  J  7 
A40 ) *SSC(40) *^OS(40) «TC(40«4)  J  8 

MsNRfl.J)  J  9 

XX=TP(T)  J  10 

F5=0S(I.J,M,VX)  J  11 

RETURN  J  12 

FNn  J  13- 
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FUNCTION  P4  ff.J)  K 

rOMMO*1  NSY5.NSSYS*NrwA£PS.LASTP.NYEAP<;.NMPM.TK(40>  *10(40.20)  .PA(40  K 
1*201  «f'(40)  *00(401  ,  H  01  (40)  , NO  (40*4)  .9(40*4). F(40«4)«UC (40*6)  *C«(40  * 

2*6) *FP (40.4.4)  *uro (40*6*4) .CPP (40 .4.4) .NAm£S  (40) *T (40) *TP (40 >  « TNP (  K 
140)  *P  (40.4)  .MCA  (40.4)  ,IJCS(40,6)  « PMC  (4o )  *C A  (40 >  «CS (40 >  *C  (40 >  *PH(40*  K 
44)  «UCAR (40*4)  *UCSP(40.6).i)CPP(40.6).K/40.?0)  ♦  IX (40.20) ♦ Y (40 .20)  *Z(  K 
540.20) *MY (40) *TAPL(40) *  XS ( 40 ) *YS(40) *75(40) *  XL (40 ) .LAMBDA (40) *101  (  K 


440)  .S«C(4f>)  .505(40)  *TC(40*4)  K 

SUM=0.  K 

NROsNPU.  I)  K 

XX*TNP(I>  K 

00  l  w=2.NR0  K 

SUM=SUM-A)  0G(05(T.J.M-1.XX)  )*CPP(!.J.m)/Fc>(I.J.M)  K 

F6=T ( T ) * ( CPR  ( I .  J. 1 ) /PR ( I . J* 1 ) ♦ ,5*SUM )  « 

PFTURN  * 

FNO  K 


1 

2 

3 

4 

5 
4 
7 
ft 
9 

10 

11 

12 

13 
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n  n  t"  o  o  o  oj  m  *-  oo"» 


5 
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SMHPOUTINF  CONTROL 
INTEGER  SOS 
DIMENSION  XTFMPU2) 

COMMDN  NSYS.NSSYS.NPMASFS.L A5TP  »NY£ARs»NMPM»TR  (40 )  .T0(40»?0)  »PA(40 
l.?0)  ,N (40 )  «RO(40)  .IN'M  (40) .NR (40,4)  »P  (40. f>)  »F  (40.6)  ,UC <40,61  .CR(40 
2.6) ,FP (40,6,6)  .1/09(40,6,6)  »CPR(40»6.6)  .NAvES(40),T(4O)  .TP (40)  tTNP( 
140)  »P (40 , 6 )  ,t  1C  A  (40 ,6)  ,UCS(40,6)  ,PMC(4n)  ,CM40)  ,CS(40)  .C(40)  .PR  (40* 
46)  »U0AR(41,6)  ,UCSR (40,6)  .l/CRR (40 .6) ,  X  f  40 .20) . IX  (40.20 )  *Y  (40.20)  ,7< 
540.  PO)  *MV(40)  .TARL(40)  ,XS(40)  , YS (40) , ?S (40 > , XL (40)  .LAMBDA  (40 ).  101  ( 
640) »SCC (40 ) » SOS (40) »  TC (40 .6 ) 

PART  t.  baseline  mcsp 

00  11  I=1.NSSYS 
Tl=12.*NMOM#MYEARS»TP(I) 

S(|M=0. 

00  1  J=1.MPHASFS 
S(IM=SUM.TO(I.J) 

T(I)=T1*S0M 

SUM=0. 

00  2  J=1.LASTP 
SUMsSUM.TO ( I » J) *PA ( I ,  J) 

TP(I)=SUM 

SUM=0. 

00  1  J=1.NPHAS£S 
SUMsSUM.TO(I,l)*PA(T.J) 

TNP(I)sSUM 

IF  (TND1 (T) .FO.l)  GO  TO  4 

NON-RSOIJNOANT  SURSYSTEM 

P(I.1)=F1 (1.1) 

(irA(I.l)aUC(T.l) 

UCS (I.l)sT(I) *CR ( I • 1 ) /F ( I  .  1  ) 

TC ( I  *  1 ) *UCA (T.l ) .DCS ( 1. 1 ) 

GO  TO  11 

IF  (P(I.l).FO.l.)  GO  TO  7 
OPFPATIVE  PEniJNOANCY 
P ( T • 1 ) =F? ( I ♦ 1 ) 

urA(i,n*o. 

MROaNP ( I . 1 ) 

00  5  l=I,MRO 

IJCA  (1,1)  =t)CA  (1*1)  +UCP  ( 1 . 1  . L » 

SUM  =  0. 

00  4  l=l,NRO 

SUM=SUM.T(I)»CRR(I.1»L)/FR( l.l.L)  w* 

UCS(I,l)=SUM 

TCU.l  )=t)rA(T.l).()CS(I»l) 


L  1 
L  2 
L  3 
L  4 
L  5 
L  6 
L  7 
L  B 
L  9 
L  10 
L  11 
L  12 
L  13 
L  14 
L  15 
L  16 
L  17 
L  18 
L  19 
L  20 
L  21 
L  22 
L  23 
L  24 
L  25 
L  26 
L  27 
L  28 
L  29 
L  10 
L  11 
L  12 
L  13 
L  14 
L  15 
L  16 
L  17 
L  18 
L  19 
L  40 
L  41 
L  42 
L  43 
L  44 
L  45 
L  46 
L  47 
L  48 
L  49 
L  SO 
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SIJ8P0UTINF  COmTPOI 


00  TO  11  L  51 

7  IF  (FP(T*l«n.NF.rR(!«l.?>)  60  TO  9  L  “1? 

C  L  53 

r  stanopy  opoumd ancy  -  all  <tbf  ape  foiial  l  54 

c  L  *>5 

P I T  •  1 ) *F1 ( I t l 1  L  56 

NDOsNPU.n  L  57 

51IM=I).  L  58 

no  8  L  =  l«NRO  |.  pq 

8  SUM3SUM*ijrPn.l«l.l  L  60 

UCA(I.1»*5IJM  L  61 

UCS ( I • 1 ) =r4  <  T  *  1 )  I  62 

Tcn.i  )*urA(T»i>»ucs<!.n  l  63 

no  TO  11  L  64 

C  L  65 

C  STANDBY  PFnijNOANCY  -  ALL  MTBF  DTFFEoFNT  L  66 

C  L  67 

9  P ( T *  1 ) *FS ( I  *  1 )  L  6fl 

NPO=NO(I*1)  L  69 

Ur.A(If  11*0.  L  70 

no  10  L=l.NRO  L  71 

if  hca(I*i  )sucAn«i)«urp<i*i*u  l  72 

UC5U»1)*56(T*1>  L  73 

Tr(T*l)s'inA(T*l)4UC5II»ll  L  74 

11  CONTINUE  L  75 

PMC (11=1.  L  76 

CA(l)aCS(l)=0.  L  77 

HO  12  I=1.NSSYS  L  78 

PMcm«PMcm*piT*n  l  79 

CA  M  1  =CA  ( 1 )  ♦t'CA  (T  *1 )  L  «0 

12  CS<1  >=CS(1>*UCS<T.1»  L  8l 

CA(1)=CA(1)*NSYS  L  82 

rsm=csm*K'SYS  l  83 

rm*CA(ll4f5|l>  L  84 

c  L  85 

C  PAPT  TI.  SUCCESS  PROBABILITY  FOP  E  ACu  OPTION  L  86 

C  L  87 

00  22  I=l.NSSYS  L  88 

IF  <  N  f  I »  .50. 0. AND. POTT). EQ. 0.1  60  TO  ??.  L  89 

IF  (N(I).50.01  60  TO  14  L  90 

N|*N ( I ) ♦  1  L  91 

00  13  Ja?*NT  L  92 

r  L  93 

C  MOM-P50UNOANT  OPTIONS  L  94 

C  L  95 

°(T  » J)=F1  f  T  *.l>  L  96 

UCA  ( I  »  J)  siJC  ( T  «  J)  L  97 

UCS  ( I » J)  =T  (  T  >  *CP (I t  J) /F  (I ♦  J)  L  98 

H  TCI  !  •  J)  *UCA  (  T •  J)  ♦UCS(  1 1  J)  l.  99 

IF  (POU).FO.O.)  60  TO  2?  L  100 
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SMRWOUTINE  COkiTPOi. 


C 

C  PFOtl^nANT  OPTIONS 

C 

Non=»n(i)*i. 
no  21  j=2.npo 

IF  (P(I*J) .FO.l.I  GO  TO  17 

c 

C  OPERATIVE  PFnuNOANCY 

r. 

PP< I » J) *F2 ( I « J) 

SI)M  =  0. 

NRI=NP(I*.J) 
no  IS  Lsl.NPT 
IS  SUMsSl'M«urPU*J«L) 

UfAH<I,J>=SUM 

StlM=0. 

ro  16  L=t *NPT 

lh  SUM=SUM^T(I)*CPP(I*  |.L)/FP(!.J»L) 

ursP(i*.n=suM 

IICPP  ( T  *  J>  =IJCAW  ( T  «  J)  ♦IICSP  ( I « J) 

GO  TO  21 
C 

r  STANOBY  PFOUNOANCY 

C 

17  IF  <FP(T«.M),NF.FP<!,J,?H  PO  TO  IP 

r 

C.  FOUAL  MTPF 

C 

PR(I*J)sF‘m.J» 

NPI=NP(I,J) 

SUM=0. 

no  1M  L=1«NPT 
IS  SUM=SUM*urR(T«J»L> 

HCA»(!« J)=SUP 
UCSP(T«.1)sF4IT«J) 

itrRp  ( i  «.n  rucAPu  •  j)  ♦••csp(  i  .j) 

00  TO  2i 
r 

C  UNEQUAL  UTRF 

C 

IP  PP 1 1  *  J)  sFS  ( T «  J) 

NPT  sPB ( I • J? 

S|IM=0. 

no  ?0  L=l»NRI 

20  SUM=S'JM>UCP(T*J*L> 

UCAR(T«  J)=S»JM 

UCS»<T*J1=F6II*J) 

'»cpp(i*j)*ucAP(i,j)*iicspn.j) 

21  CONT!*'lJF 

2  2  CONTINUE 


I 

L 

L 

L 

L 

L 

L 

l_ 

L 

L 

L 

l 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

I. 

t. 

L 

L 

l 

L 

L 

I. 

L 

L 

L 

L 

L 

L 

L 

L 

L 


nl 

02 

03 

04 

05 

06 

07 

08 

09 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
2  9 
20 

31 

32 

33 

34 

35 

36 

37 

38 

39 
AO 
A  1 

42 

43 
A4 
45 
A6 
47 
a8 
A  9 
SO 
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surroutinf  COMTROt 


c 

L  181 

c 

PART  III.  ADJU5TF0  BASELINE  SYSTEM 

L  152 

c 

L  153 

00  10  Tsl.NSSYS 

L  154 

IF  (N(  I ) .FO.O, ANO.PO  ( I )  .FO.O, )  GO  TO  7ft 

L  155 

TN=N<n*l 

L  156 

TP*R0(I)M 

L  157 

INF=1«IRE=2 

L  158 

c 

PUT  THE  ARRAYS  INTO  A  TEMPORARY  ARRAY 

L  159 

no  23  J=INF.IN 

L  160 

27 

XTFMPfJ)*P(I.n 

L  1 6 1 

no  24  JsIOF.TR 

L  162 

24 

XTFMP(IN*J-1>*PP(I*J» 

L  163 

JX*IN+ IP-1 

L  164 

TALL  TSORT  (XTFMP.J*) 

L  165 

C 

MOW  PIACP  INTO  X 

L  166 

00  28  Jsl.JX 

L  167 

X(T.J)=XTFHP(J) 

L  168 

00  26  L=INF.IN 

L  169 

IF  (P(I*L).fO.XTFMP<J»)  25. *6 

L  170 

?S 

IX ( I . J) *L*33° 

L  171 

TX ( I • J) *SUIFT (IX(I*J>.6) *558 

L  172 

Y (I ♦ J) *TC ( I »L1 

L  173 

7  (I »  J)  slJCA  ( I  »L) 

L  174 

26 

CONTIMUE 

L  175 

00  28  l*irf.tr 

1.  176 

IF  (PO(I.U.F0,XTEHP(J))  ?7«?8 

L  177 

V»7 

TXO.J)*L*330 

L  178 

IXCI*J)*SHIFT(IX(I*J>*6) ♦228 

L  179 

Y  ( I » J) =UCRR ( T«L ) 

L  180 

7  ( T» J1 »UCAR (I »L) 

L  1«1 

28 

CONTINUE 

L  182 

YMIN*1.EOOO 

L  183 

C 

nETFRMINF  THF  MINIMUM  OF  Y 

L  184 

no  29  J=1 . JX 

L  185 

if  (Y(I.j>.gt.ym;n>  no  to  ?r 

L  186 

TF  (XII. J) .LT.P(Ifl) »  GO  TO  29 

L  187 

YMJMsY  (  I  •  J ) 

L  1PB 

MY(I)=J 

L  189 

20 

CONTINUE 

L  190 

IAPLm=IX(I.MY(I>» 

L  191 

XSCI)sXIT.MY(T) » 

L  192 

YSU>=Y(T,MY(m 

L  193 

7sm»zci.MYiin 

L  194 

30 

CONTINUE 

L  105 

PMC(2)*I. 

L  196 

00  71  Ial.NSSYS 

L  197 

31 

PMC(2)=PMC<2>*XSm 

L  198 

OA (?) =C (21=0. 

L  199 

00  3?  T*1.NSSYS 

L  200 
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SUBROUTINE  COmTPO) 


r4(P)?r4(?i^sm  l  2ol 

12  r(?|rC(?)*vsm  L  2*2 

CA(P)=CA(2)*\iSYS  I  ?n3 

C(?)*NSYS*C(2)  L  ?04 

cs(2)=c(?)-ca(?)  L  2o5 

c  L  206 

C  PART  TV.  OPTIMAL  OPTIONS  L  207 

0  L  20fi 

00  15  T=1.NSSYS  L  2()9 

J*sN(I)^(I»M  L  210 

IF  (MYU).NF.JX)  00  TO  33  L  211 

*L(I)=0.  L  212 

GO  TO  IB  L  213 

31  YMAXr-l.E».100  L  214 

IF=MY(T>«TF1=IF*1  L  2l5 

00  34  J=IF1«JX  L  216 

nr=Y(T.J)-Y(T.IF)  L  217 

YMs(x(Ttj>/xr i«rn )**n./oc)  l  2ib 

IF  (YM.LT.YMax)  00  TO  34  I.  2)0 

L  AMMQA ( I ) “ J  L  2? 0 

YM/usYH  L  221 

3a  C0NTINUF  L  222 

PCsYU.LAMpnam  )-Y(T.IF)  L  223 

XL ( T )  =  ( X ( T «LAM50A (I))/X(l«IP))**(l. /Or)  L  224 

35  rONTIMIE  L  2 25 

K=?  L  226 

34  KsK^l  L  2*7 

00  17  Tsl.NSSYS  L  2 28 

IF  (XL(I).NF.O.)  00  TO  18  L  229 

17  CONTINUE  L  2l0 

MAX=K-1  L  221 

00  TO  42  L  212 

11  YMAX=-1.F*300  L  233 

00  29  T=1.NS5YS  L  214 

IF  (XL  (T)  .LT.VMAX)  00  TO  .19  L  235 

ini(*)=I  L  216 

YMAX=XL(T)  L  217 

39  CONTINUE  L  2 18 

IA=101 (K)  L  219 

IPsl.AVROAdA)  L  240 

ins^YdA)  L  241 

PMC(K)=PMr(K-l)*(X(TA*I8)/X(IA.T0>)  L  242 

CA(K)sCA(K-l)»NSYS*(7(IA«IR)-Z(IA*I0)»  L  243 

C(K)=C(K-l)*NSY5*(Y(TAf Ifl)-Y(IA.ID) )  L  244 

C S(K)=C(K)-CA(K)  L  245 

SSC(K)=IA  I.  246 

SOS (* ) =1 X ( T A . 15)  L  247 

MY(IA)=LAM*OA(!A)  L  248 

TJ=N(IA)*PO(TA)M  L  249 

IF  (MY(TA).NF.TJ)  GO  TO  40  L  250 
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SUBROUTINE  COnTROI 


Tt.(lA)*0.  L  251 

GO  TO  36  1.  252 

4ft  ymax=-1.EOOO  L  253 

Tr=MY(I4)5IFl=IF*l  L  254 

00  41  J-IF1.TJ  L  255 

PC=y<TA.J)-Y(TA.IF)  L  256 

YM= (X(IA,J)/X(IA»TF)) ** ( 1 . /OC>  L  257 

IF  (YM.LT.YMAK)  00  TO  41  L  258 

LAMRDA(TA)=J  L  259 

YMAX=YM  L  250 

41  CONTINUE  L  261 

PCrY(TA.LAMROA(IA))-Y(IA.IF)  L  262 

XL(IA)  =  (X(IA,LAMRnA(TAn/X(IA.IF)1**(l./0r>  L  263 

GO  TO  36  L  264 

C  L  265 

C  FINISHED.  NOW  PRINT  L  266 

C  L  267 

C  ADJUSTFO  PASFLINF  SYSTEM  L  268 

C  L  269 

42  PRINT  46  L  270 

00  43  I=1.NSSYS  L  271 

PRINT  47.  T .NAMES (II • I ARL ( I 1  L  272 

4.3  CONTINUE  L  2*3 

C  L  274 

C  OPTIONS  L  275 

C  L  276 

PRINT  4«  L  277 

PRINT  49,  PMC(1).CA(1).CS(1),C(1>  L  278 

PRINT  50.  PMC (2) *CA(?)*CS(?)»C(?)  L  2 79 

IF  (MAX.LT.31  GO  TO  45  L  2P0 

00  44  1*3, MAX  L  281 

print  5i.  i.PMcm.cAm.csm.cm.ssccii.sosm  l  2«2 

44  CONTINUE  L  2R3 

45  CONTINUE  L  2«4 

RETURN  L  285 

C  L  286 

C  L  287 

46  FORMAT  (lHl,QX,?4HAn.jUSTEn  RASELTNE  5ySTFm//,<»X.9H5URSYSTEM,5X.9HS  l  280 

1IIRSYSTFM.5X.1QH0PTI0N  FOR  ADJUSTED/, 3x.6HMUMBER.9X.4HNAMF.13X.8HRA  L  289 

2SFLINF)  L  290 

4 7  FORMAT  (IX.I8.6X.A10.14X.R?)  L  2«1 

4H  FORMAT  (1H1.43X.25H0PTIMAL  SUBSYSTEM  OPTIONS. //. 2X. 13HC0NFIGURATI0  L  2«2 

1N,17X.llHAC0tlI5lTT0N.7X.flHL0CiI5TlC«9X,5HTOTAL.7X.9HSUBSY5TEM.7X.6H  L  2«3 

20PTI0N/.PX.14HI0ENTIFICATT0N.6X.4HMC5O.10X.4HC0ST.RX.12HSUPP0RT  CO  L  204 

3ST.8X.4HC05T .8X.7HCHANGF0* 7X ,8HSf I.ECTrO)  L  295 

40  FORMAT  (/. 16H  RASELINE.  CI=1 .3X.F6.3,6X.rl0,2*8X.F10.2.4X.Fl0.2)  L  296 

50  FORMAT  (/.16H  A0JU5TE0.  01*2,/, 10H  PASF| INE»9X.F6.3.6X,Flft.2»RX,  L  297 

1F10.2.4X.F10,?)  L  208 

51  FORMAT  (/.4X,3HCT=«I4.8X.F6.3*6X.F10.?.RX,F10.2.4X.F10.2*7X.F5.0.9  L  209 

1X.R2)  L  300 
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SUBROUTINE  COmTPOI 


FNO 
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SUBSOMTINF  T909T  (A.N) 

M 

1 

HTMFNfTOM  A(ll.  TL  < 1  A) •  I*M  1 A) 

M 

? 

T  =  1 

3 

J=N 

M 

4 

Hrfi 

M 

5 

1 

if  u.lf.m  on  to  9 

M 

6 

2 

!J=(l*J>/? 

7 

*  =  T 

M 

8 

1=  J 

M 

9 

TF  (AMM.F.MJ))  GO  TO  3 

M 

10 

TsA(J) 

M 

U 

a ( j» -a  <  r » 

M 

12 

ACT) =T 

M 

13 

3 

T=A(IJ) 

* 

14 

TF  Hm.I.F.TI  GO  to  4 

M 

IS 

A(TJ)=A(T) 

16 

A(I)=T 

M 

17 

TsA(IJ) 

M 

18 

00  TO  5 

M 

19 

4 

IF  (T.LF.4(jn  GO  TO  5 

M 

70 

A(TJ)=A(Jl 

M 

71 

A(J»=T 

M 

?2 

T= A  ( I  i) 

M 

23 

5 

L*L-1 

M 

24 

TF  (T.LT.Ad  ) 1  GO  Tn  <5 

M 

?5 

TT=A(L) 

M 

26 

b 

«=*<♦  | 

M 

7l 

TF  (A(K).LT.T)  GO  TO  6 

M 

78 

IF  <L.LT.K)  on  to  7 

M 

29 

A(L)=*  CK ) 

M 

30 

A (K ) sTT 

M 

31 

GO  TO  5 

M 

32 

7 

MsM« ) 

M 

33 

TF  (l.-I.LF.J-O  GO  TO  8 

M 

34 

IL(^)=I 

M 

35 

TU  M )  =L 

M 

36 

T=K 

M 

37 

GO  TO  10 

M 

38 

P 

TLTM)=K 

M 

39 

TtMd=J 

M 

40 

Jsl 

M 

41 

GO  TO  10 

M 

42 

9 

TF  (M.FQ.r)  OBTUPN 

M 

43 

T  =  TI.(M) 

M 

44 

JsKMM) 

M 

45 

MrM-1 

M 

46 

Kl 

TF  T  J-T.GP.n)  GO  TO  2 

M 

47 

TF  (T.FO.ll  GO  TO  1 

M 

48 

11 

T  =  T* » 

M 

49 

TF  (  J.LT.T)  GO  TO  9 

M 

BO 
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SU»»0‘ITiWF  TSOOt  (A.N> 


t  *  A  (  f  )  M  s;  1 

if  (a(t-m.lp.T)  rn  to  u  •*  q? 

Ksr-1  VI  IT  3 

1?  A(K4l)<A(K)  vl  r,A 

KsK-1  y  c;5 

»F  (T.LT.A(K)I  00  TO  12  VI  c,b 

A(K*1)=T  m  q7 

00  TO  1 1  v|  Csij 

FMO  M  CQ- 
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LIST  OF  SYMBOLS 


Availability  of  an  aircraft. 

Ratio  of  total  operating  time  to  mission  operating  time. 

Unit  acquisition  cost  of  a  redundant  unit. 

Cost  of  the  j-th  option  for  the  i-th  subsystem  (j  =  1,  2,  ..., 
n(i);  i  =  1,  2,  ....  N$) 

Logistic  support  cost  of  a  single  redundant  unit. 

Average  cost  per  repair  of  the  i-th  subsystem. 

Average  cost  per  repair  associated  with  the  j-th  option  for  the 
i-th  subsystem. 

The  i-th  failure  mode  of  a  certain  subsystem. 

Average  yearly  logistic  support  cost  of  the  i-th  subsystem. 
Total  system  logistic  support  cost  during  y  years. 

Average  number  of  missions  per  month  per  system. 

Total  number  of  systems  (fleet  size). 

Number  of  mission  phases. 

Total  number  of  subsystems. 

Number  of  options  for  the  i-th  subsystem. 

Probability  of  abort  given  a  failure  of  a  certain  subsystem. 
Conditional  probability  of  mission  abort  given  that  the  i-th 
subsystem  fails  during  the  j-th  mission  phase. 

Conditional  probability  of  abort  due  to  safety  factors  given  a 
failure. 

Conditional  probability  of  abort  during  phase  j  given  no  abort 
before  phase  j. 

Probability  that  a  failure  of  the  i-th  redundant  unit  is  an 
abort  causing  failure. 
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LIST  OF  SYMBOLS  (continued) 


Probability  aircraft  reaches  target  and  releases  weapons  without 
an  abort  causing  failure  given  that  it  survives. 

Probability  that  a  subsystem  with  n  redundant  units  will  not 
cause  an  abort  during  operating  time  T. 

Probability  that  the  system  completes  the  2-th  mission  phase 
without  an  abort  causing  failure. 

Conditional  probability  of  reduced  effectiveness  given  a  failure. 
Probability  that  the  i-th  subsystem  completes  its  function  without 
an  abort  causing  failure. 

Probability  that  the  i-th  subsystem  completes  its  function  (i.e., 

N 

operates  for  time  t  , 

/  >  i  J 

j  =  I 

Probability  that  a  mission  is  completed  without  an  abort  causing 
failure. 

Single  sortie  survival  probability. 

Probability  aircraft  aborts  before  releasing  weapons  and  survives 
the  return  trip. 

Probability  aircraft  survives  to  release  its  weapons  on  target. 
Probability  aircraft  survives  return  trip  after  weapons  are 
released. 

Expected  number  of  repairs  of  the  i-th  subsystem  during  one  year. 
"Kill  Potential"  =  expected  number  of  targets  destroyed  after 
aircraft  reaches  the  target  area. 

Number  of  sorties  aircraft  flies  (if  it  survives). 

Operating  time  of  a  certain  subsystem. 


LIST  OF  SYMBOLS  (continued) 


Expected  number  of  targets  destroyed  after  S  sorties. 

Expected  number  of  targets  destroyed  during  the  "lifetime"  of 
the  aircraft,  i.e.,  S  -*•  «. 

Operating  time  of  the  i-th  subsystem  during  the  j-th  mission 
phase. 

Total  y-year  operating  time  of  subsystem  i. 

Operating  time  of  i-th  subsystem  (i  =  1,  2,  ...,  Ns)  during 
one  mission,  i.e.,  duty  cycle  of  i-th  subsystem. 

Aircraft  mission  time. 

Mean  time  to  restore. 

Mean  operating  time  between  failures  for  a  certain  subsystem. 

For  the  discussion  of  failure  modes  a  subscript  on  this  symbol 
would  unnecessarily  complicate  the  development. 

Mean  operating  time  between  abort  causing  failures. 

Mean  operating  time  between  abort  type  failures  of  the  i-th 
standby  redundant  unit. 

Mean  operating  time  between  abort  type  failures  of  a  redundant  unit. 
Mean  operating  time  between  failures  of  the  i-th  subsystem. 

Mean  operating  time  between  failures  of  the  i-th  redundant  unit. 

MTBF  of  the  total  aircraft  system. 

Lower  MTBF  for  the  j-th  option  for  the  i-th  subsystem. 

Upper  MTBF  for  the  j-th  option  for  the  i-th  subsystem. 

Number  of  years  to  be  considered  in  the  calculation  of  log?stic 
support  costs. 
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